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Free  Boundary  Problem  for 
the  Heat  Equation  with  Applications  to  Problems 
of  Change  of  Phase* 

I.  GENERAL  METHOD  OF  SOLUTION 

I.  I.  KOLODNER 

In  the  physics  of  change  of  phase  processes  (melting,  evaporation, 
recrystallization,  dissolution)  one  often  encounters  a  boundary  value  problem 
for  the  heat  equation  in  a  domain  whose  boundary  is  not  fully  known  but 
which  must  be  found.  It  is  customary  to  call  such  boundaries  free  bound¬ 
aries. 

1.  The  Problem 

Let  Dr  be  an  open  and— for  simplicity  of  the  argument— a  simply 
connected  domain  in  the  x,  <-plane.  Let  its  boundary,  as  in  Figure  1,  con¬ 
sist  of:  the  unknown  arc  AB  represented  by  the  equation  *  =  R(t),  the 
given  arc  AC,  the  segment  BC  of  the  characteristic  t  =  i. 


The  coordinates  were  so  chosen  that  the  unknown  arc  meets  the  given 
arc  at  a  point  (a,  0).  Without  loss  of  generality  we  may  assume  that  a  =  0, 
although  in  some  cases— as  in  examples  3-5  of  Section  6— it  is  convenient 
to  have  a  ^  0. 

The  problem  is  to  find  R(t),  0  ^  t  ^  i,  and  «(*,  t)  —  u(P),  P  e  D, 
satisfying1: 

•This  paper  represents  results  obtained  under  Contract  N6ori-20I,  T.  O.  No.  1,  sponsored 
by  the  Office  of  Naval  Research. 

'The  following  abbreviations  are  used  throughout: '  DE— differential  equation,  IC— 
initial  condition,  BC— boundary  condition,  FBC— free  boundary  condition,  JC — jump  con¬ 
dition,  qo  C— condition  for  large  |*|. 
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DE 

Lu  =  u„  —  ut  =  0, 

P'D*. 

BC 

Fu  —  0 

on  arc  AC, 

FBC1 

u(r,  t)  =  m 

for  0  ^  t  i, 

FBC2 

£ 

11 

OQ 

for  0  <  2  Sj  /. 

As  to  the  functions  /(/)  and  g(t),  we  shall  not  require  that  they  be  ex- 
plicitely  given  as  functions  of  t,  but  permit  them  to  be  functionals  of  R(t). 
In  all  the  examples  considered  in  Section  6,  they  are  functions  of  R(t), 
R(t)  etc.  We  shall  occasionally  use  the  notation  f[R],  g[J?].  We  assume 
that  if  q{t)  is  continuously  differentiable  in  the  half-open  interval  (0,  t]  and 

(1.2)  \Q\^At-^,  q(0)  =  a, 

then  f[q]  is,  as  a  function  of  t,  continuously  differentiable  in  the  half-open 
interval  (0,  f],  and  g[o]-Lipschitz  continuous,  and 

(1.3)  \f\^At-'+\  e  >  0,  |g|  ^At-'l*. 

The  nature  of  the  boundary  operator  T along  the  given  arc  AC  is  best 
described  in  terms  of  an  auxiliary  boundary  value  problem.  Let  D  be  the 
domain  of  points  P  lying  to  the  left  of  AC,  and  let  »  be  a  differentiable 
function  along  AC.  F  is  required  to  be  such  that  the  problem  for  w(x,  t)  = 


w(P),  P  e  D, 

DE 

Lw  =  0, 

PeD, 

(1.4)  BC 

r(w  +  v)  =  0 

on  arc  AC, 

IC 

w(x,  0)  =  0 

x  <  a, 

00C 

\w(x,  i)|  ^  A, 

|  Vt  wx(x,  /(  ^  A, 

has  a  unique  solution. 

In  the  following  we  shall  derive  a  functional  equation  which  the  un¬ 
known  boundary  function  R(t)  must  satisfy.  In  the  process  of  the  derivation 
we  shall  also  obtain  a  representation  of  u(x,  t)  involving  the  boundary 
function.  Thus  the  problem  will  be  solved  once  the  functional  equation  for 
the  boundary  is  solved. 

It  is  remarkable  that  in  a  problem  of  the  type  considered,  where  the 
notions  of  the  boundary  and  of  the  solutions  of  the  differential  equations 
are  interlocked,  it  is  eventually  possible  to  eliminate  all  reference  to  the 
latter.  Due  to  this  reduction,  the  onus  of  the  problem  is  transposed  to  the 
solving  of  the  equation  for  the  boundary.  While  this  need  not  be  a  trifle— 
this  equation  turns  out  to  be  a  non-linear  integro-differential  equation— we 
shall  here  adopt  the  attitude  that  an  essential  simplification  has  been 
achieved,  and  that  we  are  faced  with  a  problem  of  a  lesser  order  of  difficulty. 
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In  chapters  to  follow,  we  shall  study  and  solve  some  of  these  equations. 
Already  in  Section  6  of  this  chapter,  it  will  become  apparent  that  the 
equations  obtained  in  cases  of  physical  interest— complicated  as  they  are- 
are  all  of  Volterra  type  and  of  the  second  kind.  Hence  they  are  amenable  to 
numerical  treatment.  Further  analysis  will,  of  course,  be  required,  if  the 
solution  in  the  large  (or  bounds  for  it)  are  desired. 

We  have  here  formulated  the  simplest  free  boundary  value  problem; 
namely  the  problem  in  which  the  medium  is  described  by  a  single,  one¬ 
dimensional  heat  equation  and  which  is  bounded  by  a  free  boundary  on  one 
side.  From  the  context  of  Section  2,  it  will  become  apparent  that  various 
extensions  can  be  easily  achieved.  One  such  extension— to  two  media  each 
governed  by  a  different  heat  equation  and  separated  by  the  free  boundary— 
is  considered  in  full  detail  in  Section  7,  while  others,  also  briefly  discussed 
in  that  Section,  will  form  the  object  of  further  studies. 

Special  problems  of  type  (1.1)  have  been  considered  for  over  a  century2, 
although  the  problem  never  received  as  much  attention  as  the  corresponding 
free  boundary  problem  for  the  potential  equation.  Except  for  the 
very  special  cases  when  the  solution  can  be  found  explicitly,  the  only 
rigorous  solutions  can  be  found  in  [2,  13,  14,  3,  10,  11,  12].  The  last  three 
of  these  references  are  based  already  on  the  method  proposed  here.  The 
demand  for  solutions  resulted  in  the  development  of  various  approximation 
procedures  which  apply  to  selected  cases  with  varying  success.  For  samples 
of  these  consult  [4,  5,  8]. 

2.  Heuristic  Considerations 

Since  the  ordinary  boundary  value  problem  is  usually  considered  of 
lesser  difficulty  than  the  free  boundary  problem,  one  is  tempted  to 
give  up  one  of  the  free  boundary  conditions,  and  first  solve  the  boundary 
value  problem  so  posed  with  a  fixed  though  arbitrary  boundary,  then  to  use 
this  solution  for  the  determination  of  the  boundary  function  by  requiring 
that  the  second  free  boundary  condition  also  be  satisfied.  Let  g[t)  be  an 
arbitrary  function  satisfying  condition  (1.2),  and  denote  by  DQ  the  domain 
formed  analogously  to  DR  (see  Figure  1),  with  arc  A B  represented  by 
x  =  g(t).  The  first  part  of  the  proposed  program  is  to  find  the  functional 
«[?./]  (P).  p*dq  satisfying 

DE  Lu  =  0,  P  e  De, 

(2.1)  BC  Tu  =  0  on  arc  AC, 

BC  u(g,  t)  =  f(t)  for  0  ^  t  ^  t. 


•For  references  that  appeared  prior  to  1929  consult  the  expository  paper  of  Brillouin 
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The  second  step  is  to  select  among  all  functions  g(t)  those  for  which  u[g,  f ] 
satisfies  the  second  free  boundary  condition.  That  is,  to  solve  the  functional 
equation  for  g, 

(2-2)  ux[Q,  /](&<)  =g{t)- 


This  program  is  defeated  at  the  first  step,  since  it  is  impossible  to  deter¬ 
mine  the  functional  u[q,  /]  in  closed  form.  Consequently,  equation  (4.2) 
cannot  be  written  out.3 

A  successful  attack  stems  from  a  different  interpretation  of  the  solution. 
Denote  by  D%  the  domain  sought  (heretofore  denoted  by  DR),  and  let  DR 
be  the  complementary  domain,  P  e  DR[P:  x  <  R(t),  0  <  t  t].  Extend 
the  notion  of  the  solution  into  the  domain  D ^  by  defining 

(2.3)  u(x,  t)  =  0,  P  e  DR. 

With  this  definition,  u(x,  t )  has  the  following  properties: 


Lu  —  0,  P  e  DR{j  D%, 

ru  =  0  on  arc  AC, 

(2.4)  u(x,  0)  =  0  for  x  <  a, 

M  =  f(t), 

[«.]  =  g(*)>  o  <  t  ^  I. 


Here  [  ]  denotes  the  jump  across  the  line  x  —  R{t),  i.e., 


[m]  =  Iim  u(x,  t)  —  lim  u(x,  t),  t  fixed. 

P€D%  P€D~ 

This  suggests  that  we  seek  the  solution  among  all  functions  which  satisfy, 
in  addition  to  the  first  three  conditions  (2.4),  the  two  jump  conditions  across 
an  arbitrary  curve  x  =  g{t). 

Let  D~  and  D+  be  domains  formed  as  DR,  D%,  with  x  =  g(t)  replacing 
the  separation  line.  The  auxiliary  functional  u*(x,  t)  is  defined  by 


DE 

Lu  =  0, 

P  <  D~  u  D+, 

BC 

ru  =  o 

on  axe  AC, 

IC 

u9(x,  0)  =  0 

for  x  <  a, 

JC1 

M  =  /M 

for  0  ^  t  ^  i, 

JC2 

04]  =  gM 

for  0  <  t  ^  /. 

•An  approximation  procedure  based  on  this  program  is,  nevertheless,  frequently  used. 
In  this  procedure,  one  selects  the  zero-th  approximation  Rt  to  R  arbitrarily  and  solves  (2.1) 
with  q  =  Rc  to  get  the  first  approximation  «,  to  u.  Rt  is  now  obtained  by  solving  (2.2) 
with  g  replaced  by  Rv  and  higher  order  approximations  are  determined  by  reiterating  this 
process.  In  applications  one  takes  R ,  =  a.  Usually  then  «,  and  Rl  can  be  determined  explicitly, 
while  the  continuation  presents  already  all  the  difficulties  of  the  program.  The  sequences 
u„,  f?„,  were  shown  in  [3]  to  converge  to  the  unique  solution  u,  R  of  the  problem  in  a  special  case. 
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We  assume  that  q  is  subject  to  condition  (1.2),  and  that  q  has  a  finite  number 
of  extremum  points.  Since  g[g]  may  be  singular  at  t  =  0,  we  cannot  expect 
w®  to  have  a  bounded  x  derivative  in  D~  u  D+,  Some  restrictive  conditions 
are  nevertheless  required  if  the  auxiliary  functional  is  to  be  uniquely  de¬ 
termined. 

Let  25*  be  the  closure  of  D*  in  the  set 
£[—  oo  <  x  <  oo,  0  <t  f^i,  |  *  —  a  |  -f  |  |  >  0].  Let  D  =  25+  j  D~, 
and  let  D  be  the  interior  of  D.  It  is  observed  that  both  D  and  D  are  inde¬ 
pendent  of  g(t).  Define  uQ±  (P),  P  e  D*  by 

(2.6)  *®±(P)=lim*®«5),  QeDf. 

Q-*P 

We  impose  in  addition  to  the  requirement  of  existence  of  «®*  and  w|*  the 
following  conditions: 

|w®±|^A 

(2.7)  |-\4<t|  ^  A 

m«(-  co,  t)  =  0. 

It  will  be  shown  in  Section  4  that  the  problem  for  the  auxiliary  function¬ 
al  m8±  has  a  unique  solution.  This  solution  will  be  determined  explicitly 
provided  that  we  know  how  to  solve  a  certain  boundary  value  problem  in 
the  determined  domain  D.  Granted  the  uniqueness  and  existence  of  uQ± ,  we 
have  the  following 

theorem  1:  If  the  free  boundary  problem  has  a  solution  u(P),  R(t), 
P  e  D#  ,  then 

u  =  uR,  P  e  D+ 

and  R(t)  satisfies 

(2.8)  uR+(R,  t)  =  /[/?],  R( 0)  =  a, 

(2.9)  <+(*,  t )  =  g[R],  R( 0)  =  «. 

The  truth  of  this  statement  follows  from  our  reinterpretation  of  the  free 
boundary  value  problem,  and  from  uniqueness  of  the  auxiliary  functional. 

theorem  2  (uniqueness  theorem).  If  either  one  of  the  equations  (2.8), 
(2.9)  has  at  most  one  solution,  then  the  free  boundary  value  problem  has  at  most 
one  solution. 

Indeed,  for  every  R(t),  u{x,t)  solving  the  free  boundary  problem, 
R(t)  must  satisfy  both  equations. 

Although  we  have  thus  established  that  R(t)  must  satisfy  both  equation 

(2.8),  (2.9),  what  is  actually  desired  is  to  show  that  the  solution  of  either 
one  of  these  equations  will  also  satisfy  the  other;  for  it  would  then  be  possible 
to  prove  the  existence  theorem  for  the  solution  of  the  original  problem.  This 
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will  be  shown  in  Section  5.  Here  we  remark  that  if  u *  and  m®  were  bounded, 
we  would  have  for  R  satisfying  (2.8), 

uR~  =  0  on  the  boundary  of  DR, 

Lu*~  =  0,  P  «  Dr. 

Hence,  by  a  well  known  theorem,  see  e.g.,  [6],  Chapter  XXIX 


Consequently, 

and 


P'D~r. 


<+(R,  t)  =  g[R]  +  «?-(*.  t)  =  g[R], 

which  is  equation  (2.9).  Actually  m®  need  not  be  bounded,  but  as  we  shall 
see,  its  behavior  is  sufficiently  mild  to  enable  us  to  prove  that  uR~  =  0, 
P  e  Dr,  is  a  consequence  of  (2.8). 

The  construction  of  the  auxiliary  functional  is  based  on  properties  of 
integrals  analogous  to  the  single  and  double  layer  potentials  of  potential 
theory  which  we  now  proceed  to  study. 


3.  Integrals  Analogous  to  Simple  and  Double  Layer  Potentials 


In  this  and  in  part  of  the  next  section  we  refer  to  the  set  E,  the  interior 
of  the  set  E  defined  previously  and  to  sets  E*.  E~  is  the  same  as  D~,  while 
E+  is  the  set  of  points  lying  to  the  right  of  x  =  g{t).  Let  U(z,  f,  t,  r)  be  the 
fundamental  solution  of  the  heat  equation, 


(3.1) 


U{x,  S,  i,  x)  = 


, - exP  {-(*  -  f)74(<  ~  T)}> 

2Vn(t  -  t) 

0, 


The  simple  and  double  layer  potentials  are  defined  by 
SJLoOK*.  *)  =  fa  U(x,  q(t),  t,  r)A(r)  dr, 

K »[**](*, 0  =  J*  ut(x>  e(*)>  l>  t)A(t)  dr. 
If  a  —  0,  b  —  t,  we  use  the  simplified  notation 


(3.2) 


t  >  x 
t  <.  r. 


50><  =  S,  Dot  —  D. 

S9±,  De±  are  defined  in  E±  similarly  to  u9±. 

These  potentials  have  the  following  properties: 

1)  If  q  e  C[0,  f],  \h\  ^  A  in  the  half  open  interval  (0,  t],  and  h  is 
integrable,  then  5®  [A]  is  a  regular  solution  of  the  heat  equation  in  u 
and  is  continuous  in  E.  Furthermore,  5*[A](g,  t)  has  a  meaning  and 

S®±[A](e(/),  t)  =  S®[A](g,  t). 

2)  If  g  e  C(0,  t],  and  A  is  Lipschitz  continuous  in  the  half  open  interval 
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(0,  t]  and  is  also  integrable,  then  D*[h]  is  a  regular  solution  of  the  heat 
equation  in  E~  u  E+Q  but  has  a  jump  of  magnitude  h(t)  across  the  line 
x  —  q(().  Furthermore  the  De±  have  a  meaning  in  E± ,  DQ[h](Q(t),  t)  has 
a  meaning  for  t  >  0,  and 

0  =  t)  ±  **«.  t  >  o. 

3)  If  q  and  h  satisfy  the  conditions  of  2,  then 

sim  =  -  mhi 

4)  If  q  e  C( 0,  t],  and  h  exists  in  (0,  t],  and  is  integrable,  then 


Dim  =  -  ^  exp  {-  z*(x,  0)}  -  D*[he]  -  S°[h]. 


Here 


z(x,  t) 


»  -  eW. 

2  Vt  —  V 


z(x,  0)  = 


a?  -  g( 0) 

2  Vi 


Properties  of  simple  and  double  layers  were  established  by  Holmgren 
[7]  in  the  case  when  the  conditions  imposed  on  q  and  h  held  in  the  closed 
interval  [0,  i].  Our  extension  is  based  on  the  observation  that  for  any  d  >  0, 


—  E0  d  -f-  St  t , 

D  —  Dol  +  Dlt . 

Both  So  t  and  Dot  are  solutions  of  the  heat  equation  analytic  in  x  and  t 
for  t  >  d,  while  in  the  interval  [<5,  t]  Holmgren’s  conditions  are  satisfied. 

The  formula  for  S°  is  obtained  by  differentiation  under  the  integral 
sign.  To  get  the  formula  for  Dl ,  we  first  write 

D0[h]  =  DB[h]  -  SQ[hq]  +  Se[Ag]. 

Now 

D"[h]  -  S9[ke]  =  J  -  ^|=)  exP  {-  z2(x’  r)}h(x)  dx 

1  r* 

=  — =  zT(x,  r)  exp  {—  z2(x,  z)}h{z)  dz 

\  n  J  o 

= - -=  exp  {—  a *}  da - |  h(r)  (  f  exp  {-  a*}  da  I  dz, 

V  71  J  #<*,<)  VjlJo  U  1(1,1)  ' 

on  integrating  by  parts.  Here 


(3.4) 


*(*, 


•'>  =  { 


+  00 
—  oo 


if  x  >  g(t) 
if  x  <  g(*). 


8 


I.  I.  KOLODNER 


Differentiation  and  use  of  the  formula  for  S®  now  yield  the  result  in  the 
desired  form. 


4.  Construction  of  the  Auxiliary  Functional 

We  first  consider  a  special  case  when  D+  is  E+  ,  that  is  when  the  given 
arc  AC  consists  of  the  x-axis  for  x  >  a,  and  when/1  is  the  identity  operator. 
The  solution  in  this  case  will  be  referred  to  as  the  fundamental  part  if. 
It  is  required  to  satisfy  equations  (2.5)-(2.7)  with  the  exception  of  BC  in 
(2.6)  which  is  replaced  by 

if(x,  0)  =  0,  x  >  a, 

(4.1)  BC 

ii®(oo,  t)  =  0. 

We  will  now  show  that 

(4.2)  if  =  Ef[f\  -  5®[g  +  fg] 

is  the  unique  solution  of  this  problem. 

First  recall  the  conditions  imposed  on  g,  f  and  g,  namely: 

q  <  Cj(0,  /],  q  has  at  most  a  finite  number  of  extremum 
feC^O.i],  points, 

(4.3)  g  is  Lipschitz  continuous  in  (0,  /], 

|e|  ^  Atr*  g(0)  =  a, 

|/|  ^  At~1+e,  e  >  0, 

\g\  ^  At-^. 

Here,  and  thereafter,  A  denotes  a  positive  constant,  not  always  the  same. 

1.  v  satisfies  the  differential  equations  in  E~\jE+. 

This  follows  from  properties  1  and  2,  discussed  in  the  previous  section, 
since  the  conditions  on  g  and  h  =  f  in  the  case  of  the  double  layer,  and  on  g 
and  h  =  g  +  fe  in  the  case  of  the  single  layer,  are  satisfied. 

2.  if  satisfies  the  jump  conditions. 

Using  properties  1  and  2,  if ±  is  defined,  and 

t)  =  if(g,  t)±\f(t), 

m  =  m- 

Now,  using  the  differentiation  formulae  (properties  3  and  4)  we  have 

Vl=~2^kt  ^  {~Z'{X'  °)}  +  ^  ~ 
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Again,  the  conditions  on  g  and  h  =  g  in  the  case  of  the  double  layer,  and  on 
g  and  h  =  /  in  the  case  of  the  single  layer  are  satisfied;  by  properties  1  and  2, 
ifl*-  exists  and 

t)  =  »S(0.  0  ±  kw, 

[*4]  =  g(t)- 

3.  v*(a:,  0)  =  0  for  x  ^  a. 

Using  the  definitions  of  S9  and  D9,  equations  (3.2), 

Let  6  >  0,  0  ^  t  <  <5.  Then  |  o(t)  —  a\  ;£  2Ay/fr,  as  seen  from  conditions 
(4.3).  Assume  that 

I*  —  a|  ^  2A\/S  +  2e,  e  >  0. 

We  get,  for  finite  x, 

0  <  2e  ^  |a;  —  a\  —  2  A\/d  Si  J*  —  g(£)| 

==  |  (*  —  a)  —  (g(t)  —  a) |  \x  —  a\  +  2A\/d  ^  A, 

if  0  ^  t  <  S.  Therefore,  using  (4.3) 

I /(*)*(*,  t)  exP  {-  z2(*.  *)}|  ^  —  exp  {-(«/VPr)2}, 

V*  —  T 

l«W  +  /WeWI  exP  {“  **(*>  *)}  ^  exp  {— (e/V^r)*}- 

11,1 -A  +  eXP  {“  Wr=r'  )!}  * 

=  2  (-  +  Vtt)  A  f  exp  {—a2}  da. 

' e  '  J  e/Vi 

This  shows  that  |t>*|  0  as  t  ->  0  and  x  ^  a  since  5  and  e  are  arbitrary. 

4.  w®(±  oo,  t)  =  0. 

Let  \x  —  a\>  2A\/rt,  for  any  t.  Then 

\x  —  g(r)|  ^  jx  —  a\  —  2 A-\/i  =  k  for  0  r  ^  t 
and  k  ->■  oo  as  |xj  ->■  oo.  Now 

1  2  2  , 

1  ^  r.  **(*. T>- 


Hence 


t^r~k* 
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Also,  for  any  value  that  z  may  assume, 

1 

j  z  exp  {—  z*}  |  ^ 


Therefore, 


f(r) 


and 


t  —  x 
g(*)  +  /( T)Q(T) 

Vt  —  T 


V2e' 

/  3  \’/» 

|z3exp{- 22}|  ^  -  . 


2  exp  { —  Z2} 
exp  {  —  z2} 


^  ^  |2*  exp{— 22}|  ^  ^ 

-4r  l*exp{-z2}| 

V xk  V  T« 


oo. 


5. 

Observing  that 


l/2(*,  r)  q(x)  \ 


we  get  on  integrating  by  parts, 

1  c*  eW 


2~r  f  ”7=  exP  {—  2*(x- T)}  dx~t~r\  exP  {-  ff2l da 
WnJoVt  —  x  Vnj%{x,t) 

- -=  f  /(r)  (  f  exp  {—  <r2}  i<r)  dx. 

The  last  two  integrals  are  clearly  finite,  since  |/|  ^  At~1+e.  As  to  the  first 
integral,  we  have 

j  -p=  exp  {  -  z2(*,  r)}  ^  ■  r  A/} - > 

I  yt  —  x  vt V<  —  r 

so  that 

1  r*  g(x)  A  r*  dx  Ay/n 

Zy/n  Jo  vT— t  6XP  ~  2  *’  T  ~  J  0  \/rVT— t  _  2 

Hence,  is  finite. 

6.  IVfoJl  ^  A. 

Using  the  definitions  of  S9,  D9,  we  have,  from  the  formula  for  v* 

-  -  Bi exp{- ^ o))  exp{_ ** i) * 

+  f  r^~  z(x’  T)  exp  {“  2*(x’ T)}  dx‘ 

2vn  J  0  1  —  T 
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Denoting  the  three  terms  by  a,  /9,  y,  consecutively,  we  have 
|«|  ^  At~v' ; 

\P\£A  f  ~~j= - =  At-y'  f  da  ...  ^  At-*,  t^i. 

J  0  T-X+V<  —  T  J0°-1+eVl  —  O 

To  estimate  y  we  split  the  interval  of  integration  into  two  parts:  from  0  to 
i/2  and  from  i/2  to  i.  For  0  r  ^  i/2,, 


g(*) 
i  —  r 


z  exp  { —  z2} 


For  i/2  <  r  <  i, 


exp  {—  z2} 


2(i  -  r) 

+ 

Hence 


|g(r)|  ^  4i-* 
z  q(  t) 


<?(*) 


2(i  -  r)  2Vi^r 


iVr 


exp  {  —  z2} 


exp  {-  z2}  ^  |zt(x,t)|  exp{  —  z2(z,  r)}  + 


yr  Vt  —  t 


.  .  ^  A  r11*  dr  A  r* 

J.  v?+vJ,s 


dr 


r, 


-  ,  .  exp  {  —  z2  (as ,  r) }  |  zt  (*,  r )  |  <fr. 

H 2  V*Vi  —  r  V*  J  t/t 

Now,  since  g(i)  has  at  most  a  finite  number  of  extremum  points,  the  same 
holds  for  p(i),  and  consequently 


w  T)  -  wh=T  (k^>  ~  *w) 


may  have  at  most  a  finite  number  of  zeros,  r(,  as  r  runs  from  i/2  to  i.  There¬ 
fore, 


**(*. *)}  !*,(*>  *)\dr=* 


/*T1  B-l 

J  +  2 

J  </2  (-1 


+  f 

J  T, 


where  «  is  the  number  of  zeros.  We  now  get 

\  >\<  —  f  *  —  4.  j!_  f1  da  i  Ay/n{n  +  1)  <  A_ 

^  ~  V  t  J  o  Vo  Vi)y,Vo{\  —  a)  V't  ~~  Vt 

Adding  the  three  results, 

I  V'tvll  ^  V<(M  +  \P\  +  \y\ )  ^  a. 

This  completes  the  proof  of  the  existence  theorem.  We  now  show  that 
7.  v*  is  the  unique  solution  for  the  auxiliary  functional  in  E~  u  E*. 
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Suppose  that  to  the  contrary  there  are  two  such  solutions  v{  and  v|  and 
let 

U  =  v\  —  v%. 

v  satisfies  the  same  conditions  as  v°,  except  that  both  v  and  vx  are  continuous 
in  E. 

We  use  the  identity 

(4.4)  2  JJ  (p\dxdt  =  J  <p2dx  +  2cp(pxdt 

domain  boundary 

satisfied  by  regular  solutions  of  the  heat  equation,  see  e.g.  [6],  which  we  apply 
to  v  in  the  domain  indicated  on  Figure  2.  This  domain  is  obtained  from  the 
rectangle  |x|  ^  k,  large  k,  0  ^  t  ^  t0,  t0  <L  i  but  otherwise  arbitrary,  by 
cutting  out  a  strip  of  width  d  about  the  curve  x  —  g(t)  and  a  rectangle  of 
width  2A\/e  and  height  e  centered  at  x  =  o(0),  t  =  0.  We  now  let  d  ->  0 


and  k  ->  oo,  keeping  e  fixed.  Contributions  to  the  line  integral  along  AF 
and  A'F'  vanish,  while  those  along  DE  and  D'E'  cancel  each  other.  Further¬ 
more,  on  AB,  A'B' ,  v  =  0,  so  that  also  there  the  contribution  is  nil.  Along 
BC  and  B'C  the  contribution  is  bounded  by 


j: 


Vi 


-  dt  =  4A  Vs, 


while  along  CD  and  C'D'  we  get  at  most 


/•a+X-y  « 

j  AHx  =  2A3Vs, 


see  conditions  (2.7).  Hence  these  also  vanish  as  e  -*  0.  Passing  to  the  limit, 
we  now  get 

2  JJ  v\dxdt  +  J  F2(x,  t0)  dx  =  0, 


o  s«s«0 
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i.e.,  v(x,  t0)  =  0.  Consequently  v(x,  t)  =  0,  and  vj  ==  i%  =  if. 

We  now  consider  the  general  problem  for  the  auxiliary  functional. 
Write  m®  as  the  sum  of  the  fundamental  part  ve  and  a  complementary  part  ufl, 

(4.5)  w®  =  if  +  uf,  P  e  D~  u  Z)J. 

Substituting  in  (2.5)  and  (2.7)  and  using  the  properties  of  if  we  find  that 
uf  satisfies 


DE 

Lwe  =  0, 

Pc  A 

BC 

r(isf  4-  if)  =  0 

on  arc  AC,  Figure  1, 

IC 

uf(x,  0)  =  0 

for  x  <  e(0)  =  a, 

00C 

\uf{x,  t)J  <:  A, 

j \ftisf(x,  t)\  A 

ve  is  clearly  continuously  differentiable  on  AC,  and  therefore,  by  the  assump¬ 
tion  on  the  nature  of  the  boundary  operator  r,  see  Section  1,  equations 
(4.0)  have  a  unique  solution.  ue  is  also  unique;  for,  if  there  were  two  such 
solutions,  u\  and  u\,  then 

w\=  u\  —  if, 

=  «|  —  if 

would  both  satisfy  equations  (4.6)  and  consequently  we  would  have 

w\  =  Mlg  =3  Ilf,  u\  S3  M®  =  if. 

5.  Reduction  Theorems 

It  will  now  be  shown  that  the  free  boundary  problem  can  be  re¬ 
duced  to  solving  a  certain  functional  equation  for  the  boundary.  Consider 
the  following  equations 


(5.1) 

Jf 

M. 

il 

O 

e(0)  =  a, 

(5.2) 

<{&,  t)  =  0, 

e(0)  =  a, 

(5.3) 

e(i)«®-(e,  t)  +  2u?(q,  t)  =  o. 

e(0)  =  a. 

Explicitly,  these  equations  are 

(5.1)' 

(e.  t)  -  feW  +  “fee.  t)  =  £f[f]  (p 

+  /<?]({?.') -few 

+  “fee.  t)  = 

(5.2)' 

"fee.  0  -  feW  +  “fee.  0 

=  -  f-^  exp  {-  *»(*(<),  0)}  +  &>[£](<>,  t)  -  &WM  ~ 

+  “fee.  *)  —  °. 

(5.3)'  e(t)if(e,  t) + 2v*(q,  t)  -  fee  (<)/(<)  +  2 g(t))  +  e“fee,  t)  +  2afee.  *)  =  0. 
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Assume  that  these  equations  have  solutions  Rv  R2,  R3  respectively  which 
satisfy  the  conditions  imposed  on  Q(t),  (see  conditions  (4.3)). 

theorem  3.  If  Ri(t)  exists  for  0  ^  t  then 

m  = 

u(x,  t)  =  uR'(x,  t), 

form  a  solution  of  the  free  boundary  problem. 


P'Dt 


Proof:  Consider  uRl(x,  t),  P  e  ,  and  apply  the  identity  (4.4)  to  uRl 
in  the  domain  ABCDEFA,  Figure  2,  with  g(t)  replaced  by  Rx (t)  and  (5  =  0. 
Along  AB  the  contribution  is  zero,  since  ue(x,  0)  =  0  for  x  <  a.  Along  DE 
the  contribution  is  zero  in  view  of  the  definition  of  Rv  Contributions  along 
A  F  and  BCD  tend  to  zero  as  one  lets  k  oo,  e  -*•  0,  and  one  gets 

*!(<) 

2  J*  j*  (uRl  )a  dxdt  -f-  J  {uRl  )2  dx  =  0. 

D- 

Hence 

uRl  (x,  t)  =  0,  Pe  , 

u*(R»  t)  =  o. 

Since  all  conditions  imposed  on  u(x,  t)  are  automatically  satisfied  by  uQ(x,  t ) 
for  any  q,  with  the  exception  of  the  free  boundary  conditions,  we  have  to 
verify  the  latter.  Now 

*)  =  »RT(Rv  t)  +  f[R1]  =  f[R,l 
u?(Ri,t)  =  uF(Ri,  *)  +  g[R]  =  g[R], 

completing  the  proof. 

theorem  4.  If  Rt(t)  exists  for  0  ^  t  ^  l,  and  Rt(t)  ig  0,  then 
R(t)  =  R2(t), 

u{x,t)  =  uRt(x,t),  Pc  Dr, 

form  a  solution  of  the  free  boundary  problem. 

Proof:  Proceeding  as  in  Theorem  3,  we  find 

0^2  J  J  («*')*  dxdt  +  J  ( uR‘  )*  dx  =  J  («*'  (Rt,  t))*Rt(t)  dt  <i  0. 


\ 

* 
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Hence 

=  0,  P€DrR%. 

The  remainder  of  the  proof  is  identical  with  the  preceding  proof. 

theorem  5.  If  Rz(t)  exists  for  0  ^  t  ^i,  then 
R(t)  =  Rs(t), 

u(x,  t )  =  uR‘(x,  t),  Pe  D%t 

form  a  solution  of  the  free  boundary  -problem . 

Proof:  One  proceeds  exactly  as  in  the  proof  of  Theorem  3. 

Remark:  Clearly,  every  solution  of  (5.1)  is  also  a  solution  of  (5.2)  and 
(5.3),  and  every  solution  of  (5.3)  is  also  a  solution  of  (5.1)  and  (5.2).  As  to 
solutions  of  (5.2),  Theorem  4  in  conjunction  with  Theorem  1  show  that  only 
those  solutions  of  (5.2)  which  have  a  nonpositive  derivative  are  also  solutions 
of  (6.1)  and  (5.3).  Although  we  feel  that  in  the  latter  case  the  sign  of  the 
derivative  is  immaterial,  we  were  not  able  to  prove  Theorem  4  when  this 
sign  is  positive.  This  is,  unfortunately,  a  serious  drawback.  To  determine 
R  we  have  at  our  disposal  three  equations.  (Many  other  equations  can  be 
also  obtained.)  In  applications— see  Section  6— it  turns  out  that  equation 
(5.2)  is  usually  the  easiest  to  handle.  The  existence  of  a  solution  R  of  this 
equation  insures  the  existence  of  a  solution  to  the  free  boundary  pro¬ 
blem  only  when  R  0.  If  it  happens  that  R  >  0  for  some  t  then  the 
existence  theorem  remains  unproven,  and  one  has  to  consider  either  equation 
(5.1)  or  equation  (5.3). 


6.  Applications 


We  consider  a  number  of  examples  drawn  from  physics  of  change  of 
phase.  In  each  case  we  shall  derive  the  auxiliary  functional  and  the  equation 
for  the  free  boundary.  Some  of  these  equations  will  be  considered  in  detail 
in  the  following  papers. 

We  define 


(6.1) 

(6.2) 

and  observe  that 


v(a)  =  — =  exp  {-  a*}, 
\n 


z(x,  t)  = 


f  -  g(T) 
2  Vf  —  T 


0  ^  T  gt, 


z(x,  t )  =  lim 

T-*t 


«  -  e(*) 

2  y/t  —  t 


+  00 

if  *  >  Q(t) 

0 

if  *  =  e(t), 

—  00 

if  x  <  g(t). 

t  >  0 


(6.3) 
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1.  Motion  of  a  plane  liquid-vapor  interface. 

A  liquid  in  the  presence  of  an  undersaturated  mixture  containing  its 
own  vapor  will  evaporate,  while  if  the  mixture  is  supersaturated,  the  vapor 
will  condense.  We  assume  that  the  process  occurs  at  a  constant  temperature 
and  denote  by  g  the  saturation  vapor  density,  and  by  c0  the  initial  vapor 
density.  The  problem  is  to  find  the  interface,  *  =  R(t),  R( 0)  =  0,  and  the 
vapor  density  u(x,  t)  satisfying 


DE  Duxx  =  ut  for  x  >  R(t),  t  >  0 

BC  u(x,  0)  =  c0  for  x  ^  0, 

FBCl  u(R,  t)  =  g  for  t  ;>  0. 

To  these  we  adjoin  a  second  boundary  condition  on  the  free  boundary 
expressing  the  conservation  of  mass  across  the  boundary,  namely 

FBC2  Dux(R,  t)  =  (q  —  g)R(t)  for  t  >  0. 

Here  D  is  the  coefficient  of  diffusion,  and  q  the  density  of  the  liquid. 

On  introducing  dimensionless  variables  and  functions,  x,  l,  R(i)  and 
u(x,  i),  defined  by 

x  =  lx,  t  =  RD-H,  u(x,  i)  =  U{-X‘  R(t)  =  lR(i), 

g  — Co 

the  problem  reduces,  on  omitting  bars,  to  finding  u(x,  t ),  R(t)  satisfying 


DE 

**m  =  Uf 

x  >  R(t),  t  >  0, 

BC 

u{x,  0)  =  0, 

x  >  R{ 0)  =  0, 

FBCl 

u(R,  t )  =  1, 

t  ^  0, 

FBC2 

xut{R,  t)  =  R, 

t  >  0. 

The  specifications  are  completed  by  adding  a  condition  at  infinity  which  on 
physical  grounds  is  u(oo,  t)  =  0.  In  the  above,  /  is  an  arbitrary  length,  while 

g~c o 
e-g 


a  = 
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If  the  vapor  is  undersaturated,  c0  <  g  and  a  >  0.  In  the  case  of  supersatu¬ 
ration,  c0  >  g  and  a  <  0.  Since,  however,  g  >  c0,  it  follows  that  in  this  case 


a  = 


gp~g 

e  —  g 


> 


-  i. 


The  equality  a  =  —  1  occurs  only  when  c0  =  o,  that  is,  when  the  liquid 
cannot  be  distinguished  from  its  vapor. 

Thus,  this  problem  is  of  the  type  formulated  in  Section  1,  with  f(t)  =  1, 
g(t)  =  a_1.R  (t),  a  =  0,  while  the  given  part  of  the  boundary  is  the  positive 
a;-axis,  and  r  is  the  identity  operator.  Part  of  the  given  boundary  has  been 
removed  to  infinity,  and  the  problem  would  not  be  properly  posed  without 
some  growth  condition  as  x  ->  oo.  We  chose  here  the  conditions  u(o o,  0)  =  0 
but  the  result  would  be  the  same  if  we  had  only  required  that  |  u(x,  t)  \ 
5S  A  exp  {for2},  t  <  oo.  Conditions  (1.2),  (1.3)  on  /  and  g  are  obviously 
satisfied. 

This  problem  can  be  solved  explicitly.4  The  result  is 

u  —  a  j  exp  {—  a2}  da, 

Wt 

R(i)  =  -  2 by/t 

where 


and  b  is  the  root  of 


a  exP  {“  <rz}  daj 

a  =  2b  exp  {i2}  j*  exp  {—  a2}  da. 


This  equation  for  b  has  a  unique  real  root  provided  a  >  —  1,  and  b  —  oo 
as  a  —  1.  This  does  not,  however,  imply  the  uniqueness  of  the  solution 
of  the  problem  itself. 

It  is  observed  that  in  this  problem,  R(t)  =  0(t~l12).  The  same  singular¬ 
ity  will  occur  in  some  of  the  other  problems  considered  here;  this  is  at  the 
root  of  our  assumptions  on  the  behavior  of  f[t)  and  g(t),  see  conditions  (1.2), 
(1.3).  Had  we  imposed  more  stringent  conditions,  these  problems  would  be 
eliminated  from  the  general  analysis. 

Although  a  solution  is  known,  it  is  nevertheless  instructive  to  apply  the 
method  discussed  in  this  paper.  Since  in  this  example  D+  =  E+,  and 


x  >  0 
t  ^  0 


‘From  dimensional  analysis  one  can  conclude  that  u(x,  t)  is  a  function  of  only 
and  the  result  follows  at  once. 
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the  auxiliary  functional  ue  consists  of  the  fundamental  part  ifi  only.  Using 
formula  (4.2), 


u 


'*■ l)=i  So  -  ^r=r)  ^ t})  lT 

r,{x,t)  l  r*  oh) 

=  t]{a)da  -  —  |  - ri(z{x,  r))  dr. 

J  t(x,  0)  J  0  —  T 


«£(*.  0  =  -  ^  °))  +  ^  J  iT- 


The  equations  (5.1)-(5.3)  are  in  this  case: 
f°  1  /•*  q(t) 

a)  J?(<r)4<r  -  —  -—=  ij(z(e,  x))dx  -  £  —  0,  e(0)  =  0, 

J  t(e,  0)  J  o  V  f  —  t 

b)  +  fz~cz^’r)r,(z^’r))dr~^  =  0.  e(0)  =  0> 

c)  M  “^/.^,(*(e’T))*)  -  Vi 0)) 

+  -  f  2(e>  *)»?(2(e.  »))*  -  (l  +  -)  q(*)  =  o.  e(0)  =  o. 

a  J  o  t  ~  r  \2  a/ 


The  first  of  these  equations  is  of  the  first  kind  and  would  be  difficult  to 
handle.  The  second  one  is  clearly  of  the  form 

b')  e(<>“G[d,  e(o)  =  o, 

suitable  for  an  iterative  attack.  Its  solution  need  not,  however,  imply  the 
solution  of  the  original  problem,  unless  the  result  has  a  negative  derivative; 
see  Theorem  4.  However,  every  solution  of  the  problem  must  satisfy  this 
equation,  so  that  a  statement  concerning  the  uniqueness  of  solutions  can  be 
made  independently  of  the  sign  of  the  derivative.  We  know,  however,  that 
the  free  boundary  problem  has  a  solution,  while  the  uniqueness  question 
is  still  open.  This  can  now  be  resolved  by  showing  that  equation  b)  has  a 
unique  solution. 

It  is  now  possible  to  indicate  in  what  sense  the  auxiliary  functional  u 9 
differs  from  the  functional  u[q,  /]  defined  in  Section  2  by  equations  (4.1). 
Even  in  this  simple  case  w[g,  /]  cannot  be  determined  in  closed  form.  If, 
however,  we  further  restrict  the  class  of  functions  q  to  functions 

q  =  —  2  bVt, 


one  finds 


—  oo" <;  b  <  oo, 
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«[e,  /]  =  |  exp  {—  a2}  exp  {-  a2}  do j  ,  P  «  D+. 

For  «e,  we  get,  after  somewhat  involved  computations, 

1  /  2b  r,lx,t)  \  r‘(x,t) 

ue  =  — =(H - exp  {62}  exp{— cr2}i(r|  exp{—  o%}do,  PeD~  u  D+ 

V™\  2  J  6  /  J  xfl^t 

The  two  functionals  are  defined  as  functions  of  (a:,  t)  over  different  domains. 
They  are  both  defined  in  D+,  and 

ut+(x,  t)  =  u[q,  f](x,t) 

only  if 

(t  +  “  exp  {&*}  J  exp  {-  cr2}  doj  =  exp  {-.  <r2}  do\~\ 
or,  if  b  is  such  that 

2b  exp  {62}  J  exp  {—  or2}  do  =  «. 

For  this  value  of  b,  however,  uQ~  =  0.  Hence  the  only  function  which 
the  two  functionals  have  in  common  is  the  solution  of  the  free  boundary 
problem. 

2.  Motion  of  plane  liquid-vapor  interface  in  a  confined  medium. 

If  we  complete  the  picture  of  example  1  with  the  requirement  that  the 


vapor  remains  confined  in  the  domain  *  gi  M  .while  at  M  there  is  an  imper¬ 
vious  boundary,  see  Figure  4,  the  problem  is  identical  with  the  preceeding 
except  that: 

u(x,  t)  is  defined  for  R(t)  <  x  <  M ,  t  >  0,  ux(M,  t)  =  0. 

Here  /[g]  =  1,  g[g]  =  -,  Z)+[g(<)  <x  <  M],  a  =  0,  and 

a 


ru  = 


I  u(x,  0) 

t) 


} 


=  0 


for  0  <  *  <  M 
for  t  ^  0. 
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The  fundamental  part  is  the  same  as  in  example  1.  The  complementary 
part  can  be  found  by  the  image  method,  and  one  gets 


uft  =  vt(2M  —  x,  t), 


x  <,  M. 


Thus, 

/•“  1  r*  nlr) 

«e=  V(a)da  +  rj(a)da  —  —  —  [r)(z{x,  t)) 

J  »(*, 0)  J  i(2M-x,0)  J  0  V»  —  T 

+  t)(z(2M  —  x,  <))] 

Here  again,  equation  (5.2)  for  the  determination  of  R(t)  is  easier  to  handle 
than  the  other  two  equations.  This  equation  is 


<-(e.  t)  =  0,  e(o)  =  o, 


or 

Q{t)  =  — 0))-r,(z(2M-e,m+  f  [*(e»  *)) 

Vt  J 0 t  —  T 

—  z(2M  —  q,  t)tj(z(2M  —  q,  r))]  dr, 

6(  0)  =  0. 

Again,  the  solution  R(t)  of  this  equation  will  be  the  desired  free  bound¬ 
ary  only  if  R  (t)  turns  out  to  be  negative  and  this  is  the  case  when  a  >  0. 
When  a  <  0,  R  >  0  and  the  theorem  may  be  established  by  using  equation 
(5.3),  while  (5.2)  can  be  used  to  determine  the  solution.  The  same  remark 
will  apply  to  almost  all  the  examples  considered  here. 

The  problem  of  solving  the  equation  for  R(t)  will  be  considered  later. 
Here  we  note  that  a  good  approximation  is  obtained  by  setting  g  ==  0  in  the 
right  hand  side  of  the  equation.  Then 


R - 7=  (1  ~  exp  {-(M/Vf)2}), 

ynt 

R  ~  —  2  -^=  (Vt(l  —  exp  {—(M/'^/t)*}  +  2 M  f  exp  {—  <r2}  da]. 
V7t  J  Ml  Vi 


3.  Decay  by  evaporation  (or  growth  by  condensation)  of  a  liquid  drop. 

This  is  the  three-dimensional  analogue  of  example  1.  Assuming  that  in 
the  process  of  evaporation  (or  condensation)  the  drop  will  remain  spherical 
due  to  some  extraneous  mechanism  like  surface  tension,  and  that  the 
saturation  density  of  the  vapor  is  independent  of  the  radius,  the  radius  of  the 
drop  R(t),  anc 
mined  by 


e  density  of  the  surrounding  vapor 

u(r,  t)  will  be  deter- 

DE 

Au  =  w, 

for  r  >  R(t),  t  >  0, 

BC 

u(r,  0)  =  0 

for  r  >  R( 0)  =  1, 

FBC1 

u(r,  t)  =  1 

for  t  ^  0, 

FBC2 

ctur(R,  t)  =  R 

for  t  >  0. 

FREE  BOUNDARY  PROBLEM 


21 


The  derivation  of  these  equations  is  similar  to  that  of  example  1,  and  they 
have  been  written  in  dimensionless  form,  using  for  the  unit  of  length  the 
initial  radius  of  the  drop,  /. 

As  is  well  known,  c  =  ru(r,  t)  satisfies  the  one-dimensional  heat  equa¬ 
tion.  Using  c,  the  problem  is  now  formulated  as  follows: 


DE 

<Vr  =  Ct 

for  r  >  R(t),  t  >  0, 

BC 

c(r,  0)  =  0 

for  r  >  R( 0)  =  1, 

FBCl 

c(R,  t)  =  R 

for  t  ^  0, 

FBC2 

xur(R,  t )  =  RR  +  a 

for  t  >  0. 

In  case  of  condensation  (  —  1  <  a  <  0)  the  drop  will  grow  ad  infinitum. 
In  case  a  >  0,  it  is  expected  that  the  drop  will  disappear  in  a  finite  time  tQ, 
i.e.,  R(t0)  =  0.  For  t  >  t0  there  will  of  course  be  no  boundary  to  determine, 
and  the  free  boundary  conditions  will  be  replaced  by 

ur( 0,  t)  =  0,  for  t  >  t0, 
lim  r~xc{r,  t)  =  finite,  for  t  >  t0. 

r->0 


Here,  /[g]  =  g,  g[g]  = - f-  1,  D*  =  E+ ,  a  =  1,  and 

a 


r(u)  =  {“5*,0)  1  =  o 

|m(oo,  t)  J 


X  >  1 

t  ^  0. 

The  complementary  part  is  zero,  while  the  fundamental  part  is 
'  Q(i)  ,  ,  *_1e(r)?(T)  +  1  +  &{r)g{r)\ 


u  = 


i  j;  (igL  ,)  -  W+I  +  ^  r))  ir 

1  r*  o(t)oIt) 

—  v(zix-  °))  +  *  rj(a)da  -  —  - -  ij(z(x,r))dr. 

J  «(*,  0)  ^JoVt  —  T 


For  the  equation  for  the  boundary  one  now  obtains 
^(e-0)  -  e*  g(r)g(r) 


(A)  gg  =  —  2a J  T){o)do  —  -^T](z(e, 0))  +  J  ^~~-z{q,x)r]{z(Q,x))dx, 

f?(  0)  =  I- 


When  studying  a  free  boundary  problem  one  generally  feels  that 
an  extension  of  the  solution  into  regions  where  the  final  result  is  of  no  interest 
will  only  complicate  the  problem.®  In  this  example,  stemming  from  a  three- 
dimensional  problem,  one  almost  fears  extending  our  considerations  into 
the  region  of  negative  r.  Let  us  then  see  what  our  procedure  would  lead  to 
had  we  not  made  this  extension. 


•To  our  knowledge,  the  only  previous  exception  to  this  attitude  is  found  in  [15],  where 
analytic  continuation  is  used  to  produce  solutions  to  the  Helmholtz  problem. 
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We  restrict  the  domain  under  consideration  to  the  point  set  t  >  0, 
r  >  0.  The  domain  D~  is  therefore  the  point  set  (t  >  0,  0  <  x  <  Q(t)), 
and  at  *  =  0  we  impose  the  boundary  condition, 

«»(  0,  t)  =  0. 

The  auxiliary  functional  ue  is  easily  found  by  using  the  method  of  images. 
If  ^(x,  t)  is  the  fundamental  part— the  same  as  the  one  obtained  previous¬ 
ly-. 

w®(a;,  t)  —  o®( x ,  t)  —  r®(—  x,  t),  x  ^  0. 

Since  g(t)  >  0  for  t  <t0,ve(—  x,  t)  is  a  regular  solution  of  the  heat  equation 
in  D. 

The  equation  for  the  boundary  is  again 

t)  =  o. 

Now,  in  view  of  the  jump  conditions, 

*4“(e(0»  *)  =  ®J(e(<).  *)  -  +  i). 

while 

wjr(—  eW.O  =  »J(-  eW.  0- 

The  equation  for  the  boundary  is  therefore 

(B)  »(<)  -  2«K(e(f),  /)  -  *  +  »•(-  e(f),  /)].  e(o)  =  1. 

Equation  (B)  differs  from  the  equation  (A)  by  an  additional  term, 
2xv°x(—  g(t),  t)  and  thus  appears  more  complicated.  We  will  show  that 
every  solution  R(t)  of  (A)  also  solves  (B).  Indeed,  if  R(t)  is  a  solution  of 
(A),  then  vR~(x,  t)  =  0,  and  consequently  uR~(—  R(t),  t)  =  vR(—  R(t),  t) 
=  0.  Hence  for  n(t)  =  R(t)— and  in  this  case  only— the  left  hand  sides  of 
both  equations  coincide.  Furthermore,  since  for  x  >  R(t),  —  x  <  -m 
<  R(t),  vR+(—  x,  t)  ==  0  in  Hence, 

uR+  =  uR+,  P  « 

although  in  general,  «®  ^  w®.  It  would  be  quite  difficult  to  verify  the 
converse  had  we  not  known  it  in  advance. 

4.  A  theory  of  cloud  behavior. 

Let  the  cloud  consist  of  a  collection  (finite  or  denumerably  infinite)  of 
disjoint  spherical  drops  growing  or  decaying  according  to  the  mechanism 
assumed  in  example  3.  Let  §(  be  the  center  of  the  t'-th  drop,  R((t)  its  radius, 
R{( 0)  —  l{.  Denote  by  D^t)  the  point  set  \x  —  f4|  >  R^t),  and  by  St  its 
boundary  |x  —  §(|  =  R^t).  The  problem  is  to  determine  the  density 
function  u(x,  t)  and  the  radia  R((t)  satisfying: 
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DE 

Au  =  ut 

for  x  e9-D<W>  t  >  0. 

BC 

u  —  0 

for  x  £><(0), 

FBC1, 

u  =  1 

for  x  e  Sf,  t  0,  i  =  1,  2,..., 

FBC2j 

■£-  f  f  UriX>  *)  do>  =  RiV) 

for  t  >  0,  i  ==  1,  2,... . 

1  J  J 
xes, 

In  this  formulation,  the  problem  of  cloud  behavior  is  of  much  greater 
complexity  than  the  preceeding  problems  in  that  the  concentration  satisfies 
a  three-dimensional  heat  equation,  and  in  that  there  are  many  free  bound¬ 
aries  to  determine.  While  the  method  devised  in  this  Part  1  is  inadequate 
to  treat  this  problem,  it  will  be  dealt  with  later  when  extensions  of  our 
method  are  considered. 

For  some  purposes  the  behavior  of  the  cloud  is  satisfactorily  described 
by  assuming  that  all  drops  are  initially  of  the  same  radius  and  that  their 
centers  form  a  cubic  lattice.  The  problem  is  then  identical  with  that  of  a 
single  spherical  drop  located  at  the  center  of  a  cube  whose  boundary  is 
impervious  to  vapor.  If  the  side  of  the  cube  is  large  compared  to  the  radius 
of  the  drop,  a  further  satisfactory  approximation  is  obtained  on  assuming 
that  the  cube  is  replaced  by  a  concentric  sphere  of  the  same  volume.  The 
density  function  u  then  depends  only  on  the  distance  from  the  center  of  the 
drop,  and  the  problem  becomes  identical  with  Problem  3,  except  that 

u(x,  t)  is  defined  for  R(t)  <  r  <  M,  t  >  0, 

ur(M,  t)  =  0. 

Here  M  is  the  dimensionless  radius  of  the  sphere.  Using  again  the  function 
c  =  ru(r,  t)  the  problem  is  formulated  as  follows: 


DE 

<Vr  = 

for  R(t)  <r  <  M,  t  >  0, 

BC 

f  c(r,  0)  =  0 

for  1  =  R  (0)  <  r  <  M, 

[Mcr(M,  t)  —  c(M,  t)  =  0 

for  t  S:  0, 

FBCl 

c(R,  t)  —  R 

for  t  Si  0, 

FBC2 

<xcr(R,  t)  =  RR  +  x 

for  t  >  0. 

As  in  example  3,  /[g]  =  g,  g[g]  =  —  +  1,  but  Z)+[g(f)  <  x  <  M, 

« 

f  >  0],  a  =  1,  and 


r«  = 


u(x,  0) 
Mux{M , 


1=0 

,  t)  -  u(M,  t)  j 


for  1  <  x  <  M 
for  t  ^  0. 


The  fundamental  part  v  is  the  same  as  in  Problem  3,  and  the  complementary 
w  part  satisfies, 
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DE 

Lufi  =  0, 

P  e  D, 

IC 

ufi(x,  0)  =  0, 

x  ?S,M  , 

BC 

Mwl(M,t)  -w°{M,t)  =  -  {Mvl(M,  t)  - 

-ve{M,t)),  t^0 

ooC 

|w*(*.  Wtw%{x,t)\^A. 

One  finds, 

w°(x,  t)  =  y/tr){z{2M  -  a,  0))  +  f  v(a)  da  -  ~  f 

J  z{2M  —  x,0)  **JoVt T 

{i  o  M x llr® 

M* - M - ll  V  v(°)do 

2  —  T  2M  —  x  —  pit)  1  /  r®  \ 

■‘«sJoeWeWe,IP{ir - jfW  v^;  ’<*>*)*■ 

—  — jjj — +  *(2Jf  —  x,  r) 

The  simplest  equation  that  the  boundary  function  must  satisfy  appears  to 
be  the  equation 

Mu%-(q,  t)  -  0  =  0, 

or,  explicitly, 

*w  =  -  *(i  -  f ) 

-  ^  (i  -  (n(*(e,  o))  -  v(*m  -  e.  o))) 


+ 


-  (z(2M  -  q,  r)  +  —  ^  Tj  v(z(2 M  -  Q,  t))J  dr, 

e(  o)  =  i. 

See  [12]. 

5.  Dissolution  of  gas  bubble  in  liquid .* 

The  radius  R(t)  of  a  gas  bubble  submerged  in  a  liquid,  and  the  concen¬ 
tration  u(r,  t)  of  the  gas  dissolved  in  the  liquid  satisfy  the  following  equations 
with  c(r,  l)  =  ru(r,  t): 


DE  c„  =  ce 

BC  c(r,  0)  =  0 

FBCl  c(R,  i)  =  (1  -  l)R  +  l 


for  r  >  R(t),  t>  0, 
for  r  >  R( 0)  =  1, 
for  t  ^  0, 


FBC2  oxr(R,  t)  =  (R  +  m)R  -f  a  ^1  —  l  +  for  t  >  0. 


•A  simplified  formulation  according  to  J.  B.  Keller. 


% 
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Here,  l,  m,  x  are  constants,  and  l  ^  0,  m  ^  0,  a  >  —  1.  In  this  formulation 
the  effects  of  surface  tension  are  included.  When  these  effects  are  neglected, 
l  =  m  =  0,  and  the  problem  becomes  identical  with  Problem  3. 


Here,  f[e]  =  (1  —  l)e  +  l,  g[g]  =  (?-+  ”°g  +  1  -  l  +  L,  D+  = 

a  Q 

a  =  1,  and 


r«  -If- 01  l  =  o 

[m(co,  t)\ 


for  x  >  1 
for  t  >  0. 


This  problem  is  similar  to  Problem  3.  The  complementary  part  is  zero  and 
one  finds  for  the  fundamental  part  of  the  auxiliary  functional 

=  U*  =  -(1  -  0))  +  [(1  -  l)x  +  l]  J 

ri{o)da 
J  *{x,0) 

1  r‘  (eM  +  m)Q(r)  +  a/g_1(T) 


Vt~- 


rj(z(x,  r))  it. 


Equation  (5.2)  for  the  boundary  is  then 

»*(e.o 

(g(t)  +  nt)g(t)  —  —  2«(1  —  l)  rj(a)  da 

J  —oo 


Vt 


V(2(Q>  0))  -  xLg~l{t) 


+ 


8(0)  =  1- 


*  (eW  +  m)e( r)  +  odQ-'ir) 


t  —  T 


z{q,  t)t}(z(q,  r))dr, 


6.  Change  of  phase  without  appreciable  change  of  density  (e.g.,  recrystalli¬ 
zation). 

When  heat  is  supplied  to  a  solid  at  a  recrystallization  temperature,  a 
partial  recrystallization  will  occur.  Part  of  the  heat  will  be  used  to  supply 
the  latent  heat  of  recrystallization  while  the  remainder  is  used  to  increase  the 


R(t)\ 

recrystolized 
S,  part 
\  u  20 

t 

original 
crystols 
u  ;  0 

\ 

ux.k(t) 

X 

Figure  6 

temperature  of  the  recrystallized  solid.  The  moving  front,  R(t)  and  the  tem¬ 
perature  u(x,  t)  of  the  recrystallized  part  are  determined  in  appropriate 
dimensionless  variables  by  solving 


\  j 


i 


\ 


i 


4 


I 


1 

I  * 
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DE 

BC 

FBCl 

FBC2 


«.( 0,  t)  =  k(t), 

R(  o)  =  o, 
u(R,  t)  =  0 
M*(°»  t)  —  —  R 


for  0  <  x  <  E(t),  t  >  0, 


for  t  ^  0, 

for  *  >  0. 


Here  £(£)  S  0  is  the  heat  input  function  (in  appropriate  units),  and  it  has 
been  assumed  that  the  recrystallization  temperature  is  normalized  to  zero. 
Also  /[e]  =  0,  gj>]  =  -  q,  D+[g[t)  <  x  <  0],  a  =  0,  and 

ru  =  ux(o,  t)  —  k(t)  =  o,  t  ^  o. 


The  fundamental  part  is 

»*  =  i  f  -S===  y(2(x- T)) dx- 

J  o  V*  —  r 

The  complementary  part  is  obtained  by  using  the  method  of  images.  One 
finds 


W  =  j  f  -  e=L=^  ri(x(—  x,  t))  dx  4 
J  o  V  t  —  T 

The  equation  (5.2)  for  the  boundary  is 


^  (i^=?)  *  +  L  [Z(S'  1)1 

—  z(—  0,  x)rj(z(—  Q,  x))]dx, 


e(0)  =  0. 


7.  Generalizations.  Problems  of  Freezing  of  a  Lake  and  of 
Solidification  of  the  Terrestrial  Crust 

So  far  we  considered  the  simplest  type  of  free  boundary  problem 
for  the  heat  equation.  Two  kinds  of  complications  could  be  conceived.  First 
there  may  be  more  than  one  free  boundary.  Secondly,  the  differential  equa¬ 
tion  may  be  different  from  the  one  considered.  To  see  under  what  conditions 
our  method  can  be  applied  to  more  complicated  cases,  let  us  briefly  review 
the  main  factors  which  contributed  to  the  establishment  of  the  reduction 
theorems  of  Section  5. 

We  associated  with  the  free  boundary  an  auxiliary  functional  having 
appropriate  jumps  across  the  candidates  for  the  free  boundary.  The  con¬ 
struction  of  the  fundamental  part  of  the  auxiliary  functional  depended 
essentially  on  the  existence  of  a  fundamental  solution  of  the  equation  in  the 
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large.  Finally,  the  reduction  theorem  could  be  established  mainly  because 
of  the  uniqueness  theorem  for  the  first  boundary  value  problem  for  the  heat 
equation. 

Disregarding  then  all  technical  difficulties,  an  auxiliary  functional  could 
be  constructed  whenever  the  differential  equation  has  a  fundamental  solu¬ 
tion  in  the  large.  If  there  are  more  than  one  free  boundary,  this  functional 
will  contain  in  addition  to  the  complementary  part  the  sum  of  as  many 
fundamental  parts  as  there  are  free  boundaries.  Each  of  these  will  have 
appropriate  jumps  across  the  corresponding  boundary,  being  regular  else¬ 
where.  To  establish  the  reduction  theorem  we  would  need  again  the  unique¬ 
ness  property  of  solutions  for  the  first  boundary  value  problem. 

Clearly  then  our  procedure  can  be  extended  to  problems  involving  the 
heat  equation  in  more  than  two  independent  variables,  and  to  free  boundary 
problems  for  the  potential  equation.  The  latter  lead  to  certain  complications 
which  are  now  being  studied. 

Here  we  shall  consider  still  another  kind  of  extension.  Namely,  we 
shall  consider  problems  involving  two  media,  each  governed  by  a  different 
one-dimensional  heat  equation  and  separated  by  a  free  boundary.  Since  no 
additional  preparation  is  required  we  proceed  immediately  to  application 
to  two  still  unsolved  problems  of  great  historical  importance. 

1.  Freezing  of  a  lake  of  finite  depth. 

We  assume  that  the  lake  has  a  constant  depth  h,  and  that  its  bottom 
and  its  lateral  boundaries  are  nonconducting.  Let 


ci,2  specific  heat 
d12  density 

klt  coefficient  of  conductivity 
L  latent  heat  of  fusion 


Subscript  1  will  refer  to  conditions  in  the  layer  of  ice,  subscript  2  to  con¬ 
ditions  in  water.  We  shall  assume  that  in  the  process  of  freezing  there  is  no 
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heat  transfer  by  convection,  and— what  is  more  disturbing— that  all  the 
parameters  remain  fixed,  and  that  in  particular  =  d2  —  d.  Then  the 
temperature  u(z,t)  and  the  freezing  front  R(t)  satisfy: 


DE  aiui,xx  —ui,t’ 

BC  %(<U)  =  TX  , 


0  <  x  <  R(t)  1 
t  ^  Oj 


in  ice, 


DE 

BC 

FBC1 

FBC2 


a2Ui,x*  —  U2,t  > 

u2(x,  0)  —T2  , 
^2,  xi^i  t)~  ®  > 


R(t)  <  x  <  h 
0  =  R( 0)  <  x  <  h 
t  >  0 


in  water, 


t)  =  u2(R(t),  t)  =  0,  >  1  at  the 

kjU iiX(R(t),  t)  —  k2u2x(R(t),  t )  —LdR(t)  J  freezing  front. 


Let  us  first  consider  the  problems  in  ice  and  water  as  two  separate 
problems.  Let 

%,*(#(*),  *)  =  &(*). 

%,*(£(*).  t)  =  &(<)• 


The  auxiliary  functional  u\  for  the  problem  in  ice  is 


(*.  t)  =  j  j  -~=L[y(zx(z,  t))  -  x>  *))]  fa  +  2:riJ  ^  VW)  fa. 

2o,V« 

*  -  e(r) 


Zi(x,  t)  = 


2axV 1  —  t 


P  €  Dj.,  Q  U  ^1,  Q  > 


where  DXe  is  the  point  set  (t  >  0,  0  <  x  <  g(t))  and  D+e  is  the  point  set 
(t  >  0,  z  >  g(t)).  For  the  problem  in  water  we  get 

ul(x,  t)  =  —  ~ f  -^=L[t)(z2{x,  t))  +t)(z2(2h-x,  t))  rfr+rsf  “  ’  ’  _i?(<r)*r, 

ZJ0Vf — T  J  (x-2h)!2a,Vt 

x  -  e(r) 


■?*(*,  t)  = 


2^^  £  —  T 


P  e  Pg,  C  O  P  J  5  < 


where  Z)^  is  the  point  set  (t  >  0,  x  <  g{t)) ,  and  D^t  is  the  point  set  (t  >  0, 
g(t)  <x  <h).  The  equations  that  the  boundary  curve  must  satisfy  are: 
for  the  problem  in  ice,  t)  =  0,  or  gx(t)  =  2u\tX{g,  t),  g(0)  =  0, 

for  the  problem  in  water,  u^x(g,  t)  =  0,  or  g2(t)  =  2u2fX(g,  t),  g(0)  =  0. 

We  now  recall  that  these  are  not  independent  problems,  but  that  gx{t) 
and  g2(t)  are  related  by  FBC2.  Thus  we  get  three  equations  for  the  three 
unknowns,  gx(t),  g2(t),  and  o(t).  These  equations  are,  explicitly, 
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gi(t)  =  -  f  Th(«i(e.  *))  +*(— e.  *)v(h(-Q>  »))]* 

2T 

&(*)  =  ~r —  We.  *M**(e» T))  -  **(2A  -  e.  -  q,  t))] 

J  or  —  T 

+  [»/(^(e.  0))  -  v(zt(2‘ h-  Q,  0))], 

eW  =  4&W-TF&(<).  0(0)  =  0. 


2.  Solidification  of  the  terrestrial  crust. 

This  is  the  three-dimensional  analogue  of  the  preceeding  problem. 
Subscript  1  refers  to  the  crust,  subscript  2  to  the  liquid  interior  of  the  earth, 
all  symbols  having  a  meaning  similar  to  that  in  the  preceeding  problem. 
The  temperature  u(r,  t)  and  the  solidification  front  R(t)  are  determined  by 


DE 

a\Aux  =  ult 

R(t)  <r  <  a 

BC 

f-T 

II 

**<* 

H 

s 

t  ^  0 

DE 

s* 

II 

s” 

* 

0  <  r  <  R(t)' 

BC 

J  u2(r ,  0)  —  T2 

,  0  <  r  <  R( 0)  —  a 

o 

ii 

cT 

K 

t  ^  o 

FBCl 

ux(R,  t )  =  u% 

(*,*)  =  o  .  1 

FBC2 

l) 

ktu2  JR,  t)  =  LdR(t)  j 

in  the  crust 
a  =  radius  of  the  earth, 

in  the  liquid, 
at  the 

solidification  front. 


As  in  examples  3 — 5  of  Section  0  it  is  again  convenient  to  introduce  c  = 
ru(r,  t).  c(r,  t )  then  satisfies: 


DE 

BC 

DE 

BC 


aiCl,rr  —  ct  ’ 

Ci  (a,  T)  =  7\a, 

GyPi.rr  ~  C2,t  i 
f  c2(r,  0)  =  7>, 
}c2(0,  t)  =  0, 


R(t)  <  r  <  a 
t  ^  0 
0  <r<  R(t) 
0  <  r  <  R(0)  =  a 
t  ^  0 


in  the  crust, 


in  the  liquid. 


FBC1  cx(R(t),  t)  =  c2(R(t),  t)  =  0, 

FBC2  hc1%r{R(t),  l)  -  kjCitr{R{t),  t)  =  LdR(t)R(t). 


Proceeding  as  in  the  previous  example,  we  find  for  the  auxiliary 
functionals: 


«?(*,*)  =  -^J  [qfafa  t))  —  i/(zi(2«— *,  r)) ]  +«^ij  y(<j)do, 

P  e  t  u  D\,t , 


,{a—x)l2alV  t 
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a.  r*  s„(t)  »i»toi/a>|Vi 

“i(*.  0  =  -7 =====  !>(**(*.  »))  -  *?(*•( —  *.  *))  ]  *+  *Vb  |  »?(or)rfff 

z  J  0  v  t  —  x  J  , — wo„v; 


(*+o)/2o,V< 


^  T  /a?  —  a\  /x  +  a\“ 

+w,hy+,yj' 


vV 


Here  the  various  point  sets  are  defined  as  follows: 


(x— ai/Zd,'^  ( 

-f*  e  -^2,e  U  D£e  • 


d1q  [t  >  0,  —  00  <  x  <  $(/)], 

Di,e  P  >  °.  eP)  <»<«], 

^2,o  [^  >  0,  0  <  a;  <  g(f)], 

D;£e  P  >  eP)  <  *  <  +  00]. 


To  determine  the  position  of  the  interface,  one  has  to  solve  a  system  of  three 
equations  for  gx((),  g2{t)  and  Q(t),  namely, 


&P)  =  J  0i(e>  T)»?(zi(e,  r))  +  Zi(2a  —  p,  r)»?(z1(2a  -  q,  r))]dx 

-£ztv{Zl(e’0))' 

ftp)  =  -  f  7^-  [**(«»  *)v(zi(e.  +  *2(—  e,  ^)^(^*( —  e> »))]  * 

Jof_r 

.i,(e+2«,0)  T  a 

+  2r2  ^(<7)^0  -  — j  0))  +  *?(z2(e  +  2a,  0))] , 

J  *,(e,0)  “svt 

eeP)  =  7T  &P)  ~  tt  &P)  >  e(°)  =  «■ 


In  connection  with  the  last  two  examples  we  note  that  Rubinstein  [14], 
obtained  in  a  similar  problem  the  same  equations  for  the  free  boundary  that 
we  would  get  by  applying  the  method  of  this  paper.  The  problem  considered 
by  Rubinstein  is  a  generalization  of  problems  of  Stefan  (see  [1]),  namely  to 
find  x ,  /),  u%(x,  t),  and  R(t)  satisfying: 


DE 

s' 

II 

8 

t>  0,  x  <  R(t), 

IC 

«i(x,  0)  =  <px{x), 

x  <  0, 

DE 

4H,xx  =  H.U 

t  >  0,  x  >  R(t), 

IC 

H(x>  °)  =  9’*(x). 

x  >  0, 

FBCl 

FBC2 

*i(flp)>  t)  =  ut(R(t),  t)  =  0, 
*i«],«{^P).  0  ~  k»H,x(R(t)>  t) 

=  kR(t)- 

(When  fi(x)  and  9 ?a(as)  are  both  constants,  a  solution  can  be  obtained  in 
closed  form.) 
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It  was  first  observed  that  if  a  solution  to  the  problem  exists,  then 
and  u2  have  the  following  representations  derived  from  Green’s  formula: 

t  o 

«i(*.  *)  =  ^  J*  t))dx  +  J  «Pi(^)  exp  {-[(*-  £)/2  a^t}*}  d£, 

0  *  <  R(t), 

t  00 

ua(x,t)=~^j-^^=v(zt(x,r))dr+^^j  <p2(£)  exp{— [(*-f)/2 

0  0  *  >  R(<), 

where  gt(<)  =  ul>x(R(t),  t),  g2(t)  =  u2lx(R(t),  t).  The  equations  which  the 
boundary  function  must  satisfy  were  then  derived  by  requiring  that 

lim  ul  x(z,  t)  =  gl(t), 

lim  u2tX(x,  t)  =  g2{t), 

R(t)<x-+R(t) 

and  adjoining  the  FBC2.  We  note  that  this  attack  would  prove  unsuccessful 
if  the  domain  of  definition  of  or  u2  were  finite. 
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On  a  Class  of  Solutions  of  Maxwell’s 
Electromagnetic  Equations 
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§1.  In  discussions  of  electromagnetic  propagation  in  a  homogeneous 
isotropic  conducting  medium  [1],  it  is  frequently  found  convenient  to  regard 
the  solutions  as  built  up  from  the  real  parts  of  component  plane  wave  solu¬ 
tions  having  the  ’’damped  monochromatic”  form 

E  =  A  exp  {i(a)t  —  kn  •  x)  —  a„f  —  £0n  •  x}, 

H  =  B  exp  { i(cot  —  kn  •  x  —  y>)  —  —  /?0n  •  x}. 

Here  E  and  H  are  the  electric  and  magnetic  vectors  respectively,  A  and  B 
are  thus  real  constant  vector  amplitudes,  and  n  is  a  unit  vector  in  the  direc¬ 
tion  of  propagation;  k  and  are  the  propagation  and  absorption  constants, 
co  and  a0  are  the  frequency  and  decay  constants,  and  y>  is  a  constant  phase 
difference.  A  and  B  may  be  shown  to  be  orthogonal. 

The  form  of  such  solutions  suggests  that  one  might  obtain  more  general 
ones  by  substituting  a  real  general  phase  function  <p(x,  y,  z)  for  the  linear 
function  —  &n  •  x  and  real  general  vector  amplitude  functions  E0(a;,  y,  z) 
and  H0(a;,  y,  z)  for  the  factors  A  exp  {— /?0n  •  x}  and  B  exp  {— /?0n  •  x}, 
respectively.  The  quantity  yi,  too,  might  be  allowed  to  depend  upon  position. 
If  such  substitutions  are  made,  the  resulting  expressions  represent  "waves” 
which  are  in  general  no  longer  plane,  and  do  not  have  constant  amplitudes. 
It  is  clear  however  that  Maxwell’s  electromagnetic  equations  will  impose 
certain  restrictions  on  <p,  ys,  E0  and  H0  ,  which  it  is  the  purpose  of  this  paper 
to  investigate. 

§2.  Maxwell’s  equations  for  a  homogeneous  isotropic  conducting 
medium  are 

„  „  3H 

VxE  +  fi'di  =  0’ 


(1) 


VxH  — £-=  ctE, 
dt 

V  •  E  =  V  •  H  =  0, 


♦Work  performed  chiefly  while  Professor  Hudson  was  in  residence  at  the  Navy  Elec¬ 
tronics  Laboratory,  summer  of  1953. 
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fi,  e,  and  a  being  the  permeability,  dielectric  constant,  and  conductivity, 
respectively,  of  the  medium.  We  look  for  electric  and  magnetic  vectors  E 
and  H  satisfying  (1)  of  the  form 

E  =  E0  exp  {-  clJ  +  i(<p  +  cot)}, 

H  =  H0  exp  {—  ccgt  -f  i(<p  —  y>  +  cot)}, 

where  E0,  H0,  <p  and  f  are  real  functions  of  position  and  a 0  and  co  are 
real  constants. 


§3.  Certain  preliminary  results  will  be  needed. 

lemma  1.  [2]  A  necessary  and  sufficient  condition  that  a  vector  field  V 
be  of  the  form  V  =  /VA,  where  f  and  h  are  scalar  functions  of  position,  is  that 
V  •  V  X  V  =  0. 

LEMMA  2.  If  (i)  |V./|2  =  1, 

(ii)  v/  •  vjf  =  v<?  •  vs  =  vse  •  vs  =  o, 

(ih)  \v<tr\  =  \v*r\g(S), 

where  <?,  Sf ,  S  are  independent  scalar  functions  of  position,  and  g  is  a  non¬ 
vanishing  function  of  S  only, 


then  V*/  =  V2.?f  =  0. 

VS  X  VS? 


Proof:  We  first  note  that  VS  =  ± 
Then  we  consider  SVS.  We  have 


\VS\\VST\ 


v  x  (SVS)  =  vs  x  vs  =  ±  gvsr. 

Thus  V  •  (gVST)  =  g'VS  •  VST  +  gV2ST 

=  gVW 

=  ±  V  •  V  X  (SVS) 

=  0. 


Hence  V*Jf  =  0.  Similarly  V2/  ==  0. 

lemma  3.  Necessary  and  sufficient  conditions  that  there  exist  functions 
S,  S? ,  S  and  g(S),  satisfying  the  conditions  of  Lemma  2  are  that  either : 

(i)  g(S)  =  constant,  and 

(ii)  the  line  element  in  Euclidean  space  must  be  expressible  in  the  form 

(3)  ds *  =  p*(S,  #)  ( mS  +  n)*(df*  +  dJt?*)  +  dS* 


where 

(4) 


la/* 

r  djrv 

p  +  m2p2  =  0 
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and  m  and  n  are  constants,  or.  (i)'  g  or  g"1  is  a  linear  function  of  J ,  and  (ii)' 
the  line  element  is  expressible  either  in  the  form 


(3)' 

ds 2  =  d£2  +  J2dj£2  +  dJ2 

or 

(3)” 

ds2  =  S2d£2  +  dtf2  +  dJ2. 

Proof:  By  the  conditions  of  Lemma  2.,  we  may  take  the  quantities  £, 
JC  and  J  as  coordinates  in  a  triply  orthogonal  coordinate  system.  If  we  denote 
the  corresponding  scale  factors  by  h^,  h2  and  then,  since  |V«?|  =  hf1, 
jVJfj  =  Ag1,  and  |Vc/|  =  kg1. 
we  have 


(5)  h2  —  ghy ,  h3  —  1. 

Now  the  scale  factors  must  satisfy  the  conditions  of  Lame  [3],  (These 
are  the  conditions  for  the  vanishing  of  the  curvature  tensor  of  an  orthogonal 
line  element  in  three  dimensions).  For  brevity,  let  a  subscript  1,  2,  or  3 
following  a  comma  denote  differentiation  with  respect  to  £,  ££ ,  or  J , 
respectively.  Then  the  Lame  conditions  for  our  case  are: 


(6) 


h  hj 

K.taK  +  K.nK  —  K.iK.x  +  *i,8*a,a  ~  > 


^2^1,23  —  ^l,2^a, 3  >  ^1,33  — 

^1^2,13  =  ^2,1^1, 3  >  ^2,33  =  0- 


We  see  immediately  that  \  and  h2  have  the  forms 


(7) 


hx  =  x./  -f-  /? 

h2  —  yS  4*  <5 


where  a,  p,  y,  and  d  are  functions  of  £  and  only.  Furthermore  we  must 
have  from  the  last  two  Lame  conditions  on  the  left 


(8)  a<5fl  =  0ytl  <5a  x  =  yPt, 
since  they  must  be  true  for  all  J.  Now  we  have 

(9)  (yS  +  <3)  =  g{*S  +  P), 

and  differentiation  of  this  with  respect  to  £  or  yields  equations  from 
which  g  can  be  eliminated  using  (9)  (since  g  ^  0).  Then  identification  of 
coefficients  of  powers  of  J  yields  linear  homogeneous  equations  in  the  deri¬ 
vatives  of  a,  p,  y,  and  6.  If  we  consider  these  in  conjunction  with  (8)  we  find 
that  either  (a)  cad  =  py  or  (b)  a,  /?,  y  and  6  must  all  be  constants. 

In  the  latter  case,  simple  substitution  into  the  first  Lame  condition 


J 
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shows  that  either  a  =  0,  in  which  case  g  is  a  linear  function  of  J ,  or  y  —  0, 
so  that  g-1  is  linear  in  J ,  or  both  a  and  y  vanish.  If  a  =  y  =  0,  then  g  =  d/p 
=  constant,  which  will  be  considered  under  case  (a).  In  no  event  can  <5  and 
P  vanish  simultaneously.  Now,  if  g  or  g_1  is  a  linear  function  of  J ,  it  is  clear 
that  a  redefinition  of  J  in  terms  of  this  linear  function,  and  a  renormaliza¬ 
tion  of  £  and  will  bring  the  Euclidean  metric  into  one  of  the  forms  (3)' 
or  (3)".  Note  that  the  surfaces  J  =  constant  for  these  metrics,  are  circular 
cylinders  along  whose  generators  either  £  or  Stf  varies. 

In  case  (a),  we  quickly  see  that 

yK  =  *h2, 

<5/q  =  ph2 

so  that  by  (5),  g  =  constant.  A  rescaling  of  the  variable^5  allows  us  to  take 
this  constant  equal  to  unity.  Hence  also 

y  —  *, 

«*  =  /?. 

Equations  (8)  now  imply  adp  —  p den ;  this  is  satisfied  if  and  only  if 

a  ==  mp\£ ,  Jf), 
p  =  np(£,Jf) 

where  m  and  n  are  constants.  Simple  substitution  into  the  first  Lame  con¬ 
dition  shows  that  p  must  satisfy  the  equation  (4). 

From  (7)  we  have 

hx  =  h2  —  {mS  +  n)p{£,  3^) 

so  that  the  line  element  is  indeed  of  the  form  (3).  The  situation  for  which 
a  and  y  both  vanish  and  P  =  d  —  constant  is  obviously  a  special  case  for 
which  m  —  0  and  p  —  constant. 

lemma  4.  [4]  The  surfaces  J  =  constant  in  a  triply  orthogonal  system  of 
surfaces  £  =  £(x,  y,  z),3f  =  Jf(x,  y,  z),  and  J  =  J{x,  y,  z),  chosen  so  that 
the  Euclidean  line  element  can  he  written  in  the  form 

ds 2  =  S(S) &{£,*?)  d£*  +  g,{£,Jf)  djf*}  +  dS* 

are  either  spheres  or  planes. 

§4.  It  is  convenient  to  introduce  the  notation 

+  il^  =  fi(w  -f  in 0)e~ir, 

(10)  k2  +  ila  =  [e(<u  +  i«o)  —  iaY'r, 

E0  —  |E0|  ,  H0  =  j  H0|. 

Note  that  4  =  dkjdtp,  la  —  —  dkjdip.  Then  we  may  state  and  prove: 


* 
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theorem  1.  Necessary  and  sufficient  conditions  that  electromagnetic 
vectors  of  the  form  (2)  he  solutions  of  the  system  (1)  are  that 


(11) 


(12) 


(13) 


(i)  f  is  a  function  of  y, 

....  „  E0  X  H0  ,  k,k 2 

(u)  Vy  =  <? - — — ,  where  =  - 

E0H ©  1  —  y 

(iii)  E0  and  H0  are  orthogonal  vectors  for  which 


E 2  = 


W') 


Hi, 


(iv)  E0  and  H0  have  the  forms 

E0  =  exp  I  —  f  ~-d<p  V«?,  H0  =  exp  j  —  f  (1  —  ip')dq> 

l  J  <P0  "1  l  j  9>0  *2 


where  y0  is  an  arbitrary  lower  limit  and  £  and  are  scalar  function  of  position. 

Proof:  If  we  substitute  the  expressions  (2)  into  (1)  and  separate  the 
resulting  equations  into  real  and  imaginary  parts,  the  system  becomes 


(14) 


V  X  E„  ^H0  —  0, 

V  x  H0  -f-  /2E0  =  0, 
Vy  X  E0  +  ^H0=  0, 

(Vy  —  Vy)  X  H0  —  AjE0=  0, 


o-0- 


V-E, 

V  ■  H0  =  0, 
V9J  •  E0  =  0, 
(Vy  -  Vy)  ■  H0  =  0. 


The  equations  in  the  left  hand  column  will  be  termed  "vector”  equations, 
and  those  on  the  right  will  be  called  "scalar”  equations. 

From  these  equations  we  infer  the  additional  orthogonality  relations 

Eo  ‘  Ho  =  Vy  •  H0  =  Vy  •  H0  =  Vy  •  E0  =  0. 

Hence  the  vectors  Vy,  Vy,  E0  X  H0  are  parallel.  Thus  y  is  a  function  of 
<p  and  we  may  write 

dw 

Vy  =  — —  Vy  =  y'Vy. 
dip 

Also  Vy  has  the  form  (11).  For,  j  Vy|  E0  —  |ftj|  H0  and  |  Vy|  H0  = 
|&,/(1  —  y')  |  £0  from  the  last  two  vector  equations  in  (14).  Hence  |Vy|*  = 
l^i^s/(l  —  y')|,  from  which  (11)  follows  immediately.  Moreover  we  see 
that  E0fH0  =  |  (^1/^s)(l  —  y')|*  which  is  just  (12). 

We  note  that  E0  •  V  x  E0  =  Ho  •  V  x  Ho  =  0.  Application  of  Lemma  1 
permits  us  to  write 


(15) 


E0  =  eaV£, 
H0  =  e'VjP. 


5 


I 


l 
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Substitution  of  these  expressions  into  the  vector  equations  in  (14)  yields 

V93  x  V<f  =  —  kyVJfC e},~a , 

VPxV^=  — V3ea~\ 

1  —  ip 

Va  X  Vc f  =  l1V^ei~i,> 

Vi  x  VJf”  =  —  l^3ea~\ 

Clearly  a  is  a  function  of  9)  and  3  only,  while  i  is  a  function  of  9?  and  JP 
only.  Moreover  these  equations  show  that 

da  4 
3p  =  ~kl’ 


db 


4 


d<p  &2  ^ 

which  are  functions  of  <p  only.  Hence 

4 

b 


v>') 


p*  1 

a  —  ”  ~r  d<p  -f- 


r-  (1  — ?')<*!» +  &o 


where  99,,  is  any  arbitrarily  chosen  lower  limit.  By  examining  the  resulting 
expressions  for  E0  ,  H0  we  find  that  the  factors  e“»,  eto  may  be  absorbed  into 

V3,  Vjti?  (by  replacing  3,  by  jea°d3,  jeb*dJ(?).  Thus  we  finally  have 

p9  l 

a  =  a(q>)  =  —  j-d<p, 

(16)  J  y 

b  =  b(<p)  =  -  f  ~  (1  —  ip')  dip. 

Kk* 

We  now  note  that  (13)  is  a  consequence  of  (15)  and  (16),  which  completes 
the  proof  of  the  necessity  of  Theorem  1. 

As  an  immediate  consequence  of  the  above  we  see  that  if  J  is  defined 
by 

r*  ll  —  in'  |H 

(17)  J  = 


=  r  Lul! 
J  <pt 


dip 


where  <px  is  arbitrary,  the  functions  J ,  3,  and  all  satisfy  the  conditions 
of  Lemma  2.  For,  g  =  |  V3\  j\  VJfj  is  defined  by 

(18)  s  =■  exp{  J’  [i  -  £  (1  -  ✓)]  ^  (1  -  ¥) 


* 

i 


SOLUTIONS  OF  MAXWELL’S  EQUATIONS 


39 


and  this  is  a  function  of  J  through  (17).  In  proving  the  sufficiency  of 
Theorem  1  by  direct  differentiation  and  substitution,  use  must  be  made  of 
Lemma  2.  We  formalize  this  step  in  the  next  result. 

theorem  2.  The  "potential”  functions  8  and  3^  of  Theorem  1  satisfy  the 
equations 

(19)  Va«?  =  V2^f  =  0. 

The  surfaces  of  constant  phase,  <p  =  constant,  are  also  surfaces  on  which 
J  —  constant.  Hence,  in  view  of  Lemmas  3  and  4,  we  have  the  additional 


theorem  3.  The  surfaces  of  constant  phase  <p,  for  electromagnetic  waves 
of  the  type  (2),  are  either  concentric  circular  cylinders,  concentric  spheres,  or 
parallel  planes. 

In  proving  the  above  theorem,  it  will  be  seen  that  the  function  g  in  (18)  is 
not  unrestricted  but  is  either  a  linear  function  of  J ,  the  reciprocal  of  one,  or 
a  constant.  By  suitably  choosing  q>0  and  <px  we  may  infer  in  any  case  that 


^(1_V') 

«« 


are  the  corresponding  three  rather  complicated-appearing  integro-differential 
equations  for  y>.  The  exponent  on  the  left  is  +  1  if  g  is  a  linear  function  of 
J ,  so  that  8  varies  along  the  cylinder  generators;  it  is  —  1  if  g-1  is  a  linear 
function  of  J ,  so  that  varies  along  the  cylinder  generators;  it  is  0  in  the 
case  of  spherical  or  plane  waves.  Simpler  forms  result  if  we  recall  that 


(20a) 


dJ 

dip 


1  -w' 


Then  the  foregoing  equations  (20)  can  be  rewritten,  after  a  logarithmic 
differentiation  in  which  it  must  be  remembered  that  k[  =  l^xp'  and  k'2  = 
—  lty)',  in  the  form 


(20b) 


d 

dip 


These  results  lead  to 


theorem  4.  The  phase  difference  function  y>(<p)  in  Theorem  1  must  satisfy 
the  integro-differential  equations  (20),  or  the  differential  equations  (20a,  b). 

The  application  of  Lemma  3  in  the  proof  of  Theorem  3  may  also  be 
formalized  by 

theorem  5.  A  necessary  and  sufficient  condition  that  the  functions  8, 
Jf,  q>  of  Theorem  1  exist  is  that  the  line  element  in  physical  ( Euclidean )  space 
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is  expressible  either  in  one  of  the  two  forms  f3)' 

(i)  ds*  =  d£*  +  S'dJP*  +  dS' 

{for  cylindrical  waves  with  the  electric  vector  parallel  to  the  generators )  and  (3)" 

(ii)  ds*  =  S*d£*  +  dJf*  +  dS* 

( for  cylindrical  waves  with  the  magnetic  vector  parallel  to  the  generators)  or  in 
the  form  (3)  i  , 

(iii)  ds 2  =  pi(^,^)(mJ’  -f  n)*(d£*  +  dJff2)  +  dJ * 

(for  spherical  and  plane  waves)  where  J  is  given  by  (17),  p  satisfies  (4),  and 
m  =  0  for  plane  waves,  or  m  ^  0  for  spherical  ones. 


As  we  shall  see,  this  last  theorem  suggests  a  convenient  method  for 
writing  down  or  generating  various  classes  of  solutions  of  Maxwell’s  equa¬ 
tions  for  spherical  and  plane  waves. 

§5.  The  foregoing  theorems  bring  out  several  points  of  physical  and 
mathematical  interest.  If,  for  simplicity  one  fixes  attention  on  the  usual 
plane  wave  solutions,  in  which  the  phase  difference  y  is  a  constant,  then  it 
can  be  seen  that  such  solutions  correspond  to  a  singular  solution  of  the  non¬ 
linear  integro-differential  equation  in  (20)  with  the  exponent  0,  namely 


(21) 


V  = 


c  -)-  d 
2 


where 

c  =  -  arg  [e  (to  +  *«„)  -  io], 
d  =  arg  [y(to  +  ia0)]. 

The  interesting  question  then  arises:  if  the  phase  difference  of  a  solution 
differs  slightly  from  (21)  for  some  particular  <p- value,  will  the  perturbation 
decay  or  grow  as  one  proceeds  in  the  direction  of  propagation  toward 
decreasing  y- values?  Put  in  other  terms,  is  the  phase  difference  of  the  usual 
plane  wave  solutions  in  stable  or  unstable  equilibrium? 

An  elementary  investigation  shows  that  the  latter  alternative  is  the  case 
in  ordinary  conducting  media.  For  simplicity  let  us  take  o,  e,  y  as  real  and 
<*„  =  0.  Then  d  —  0  and  c  =  tan-1  o/we.  Now  if  n  is  a  small  deviation  of  y 
from  the  value  c/2,  n  must  satisfy  the  differential  equation 


(23) 


n”  +  4tt'  tan  +  4rr  |l  tan2— j  =  0. 


Since  <p  increases  negatively  in  the  direction  of  propagation,  n  is  an  oscillating 
increasing  function  in  this  direction,  and  the  singular  solution  is  unstable. 
This  peculiar  behaviour  of  the  phase  difference  between  the  E  and  H  vectors 
is  to  be  investigated  in  considerably  more  detail  in  another  paper,  in  which 
the  physical  consequences  are  discussed.  We  merely  point  out  here:  1)  the 
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restriction  to  plane  waves  is  not  essential  to  the  above  argument,  2)  the 
form  of  the  vectors  E0  and  H0  indicates  that  the  amplitudes  of  these  waves 
also  depend  in  a  curious  way  on  the  phase,  and  3)  the  phase  difference 
behaves  differently  for  cylindrical  waves  as  compared  to  the  other  types.  In 
particular,  for  cylindrical  waves  of  our  type  no  solution  for  which  y>  —  con¬ 
stant  seems  to  exist. 

The  next  point  of  interest  is  that  in  the  case  of  spherical  or  plane  waves, 
the  curves  S  •=  const  and  JS  =  const  form  an  isometric  orthogonal  net  on 
surfaces  of  constant  phase  <p,  but  are  otherwise  rather  general.  Moreover  if 
one  pair  of  functions  (S,  3if),  and  a  solution  p  (S,  3^“)  of  4  are  known  all 
other  su.ch  pairs  and  solutions  may  be  determined  by  the  use  of  harmonic 
functions  [5],  For,  let  w  be  any  harmonic  function  of  f  =  S  -f  itff.  Then  a 
new  pair  (S',  2tf")  is  determined  by 


(24) 

and  the  function 

(25) 


-f-  i3^S'  =  w(S  -f- 


P'(S',JS')  =  p(S,3f)  \dwjdn 


is  also  a  solution  of  an  equation  of  the  form  (4),  written  in  terms  of  S'  and 
3%” .  This  is  easily  checked  by  direct  substitution.  For  any  such  solution 
p(S,Jf)  of  (4),  the  magnitudes  of  the  gradients  of  S  and  Jf?  are  of  course 
given  by 


(27) 


(26)  |  VS  |  =  |  VJf  |  =  p~'(S,  +  n)-1 

where  J(x,  y,  z)  =  constant  is  a  suitable  sphere  or,  if  m  —  0,  a  plane.  In 
addition,  if  1/  and  ijt  are  unit  vectors  in  the  directions  of  the  old  electric  and 
magnetic  vectors  respectively,  and  we  write  dwjdt,  —  \dw/dC\  exp  {«?}, 
then  for  the  new  solution,  the  corresponding  new  directions  are  given  by 

i/'  =  if  cos  T)  —  ijt  sin  rj, 
ijT'  =  if  sin  r\  -f  ijf  cos  rj. 

Thus  we  see  that  the  factors  VS  and  Vjf  appearing  in  the  solution  (13)  of 
Maxwell’s  equations  are  determinable  for  our  spherical  and  plane  wave  cases 
by  the  methods  of  conformal  mapping.  For  the  cylindrical  wave  case  we 
obtain  only  solutions  for  which  essentially 

either  |V^|  =  J~x,  j  V3tf>\  =  1 

or  |VJf|  =  |V<f|  =  1. 

In  order  to  calculate  explicitly  the  exponential  factors  in  (13),  a  solution 
rp(<p)  of  equations  (20)  must  be  found,  and  then  <p  must  be  determined  in 
terms  of  J  by  inversion  of  (17). 


(28) 
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§6.  We  conclude  with  several  elementary  illustrative  solutions  of 
Maxwell's  equations  in  the  form  (13).  In  these  examples  we  shall  not  con¬ 
cern  outselves  with  a  general  determination  of  <p  (z,  y,  z),  leaving  this  for  a 
later  paper,  but  will  consider  only  solutions  for  which  y>  —  constant.  This 
automatically  excludes  the  case  of  cylindrical  waves.  In  such  cases  we  may 
take,  with  q>0  =  <Pi , 

(29)  <f  =  +  9>0 

so  that 

(30)  E0  =  V<?  exp  {—  WiM,  H0  =  exp  {—hlKIKV'-^} 
where  lx  ,  l2  ,  k t,  k2  are  all  constant.  Hence  we  need  to  determine  only  the 
explicit  forms  of  J ,  Vg,  and 

As  a  first  example  consider  a  cylindrical  coordinate  system  (g,  6,  z), 
in  which 

(31)  ds*  =  d e2  +  Q2  dd*  +  dz*. 

Rewrite  this  as 

(32)  ds 2  =  g*[(d  In  e )2  +  d0*]  +  dz*. 

Apparently  we  may  take  m  —  0,  n  —  1,  p  =  g  and 

g  =  In  q, 

=  0, 

S  =z. 

Then  the  electric  and  magnetic  vectors,  by  formulas  (26)  and  (30),  are  pro¬ 
portional  respectively  to 

and  the  waves  travel  along  the  z-axis.  This  corresponds  to  the  well-known 
solution  for  plane  waves  travelling  along  a  wire  conductor  coincident  with 
the  z-axis. 

As  the  second  example,  consider  a  spherical  coordinate  system  ( r ,  0,  A) 
in  which  the  line  element  may  be  written  as 

ds 2  =  dr 2  -}-  ra  sin*  A  |ji0)2  +  J- 

Then  we  may  take  tn  —  1,  n  =  0,  p  =  sin  A  and 

.  ,  * 

g  —  In  tan—  , 

2 

sue  =  e. 

J  =r. 
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so  that  the  waves  are  travelling  radially  and  the  electric  and  magnetic  vec¬ 
tors  are  both  proportional  to  1/r  sin  X,  but  point  in  the  polar  and  azimuthal 
directions  respectively.  Imagine  a  perfectly  conducting  surface,  in  the  shape 
of  an  infinite  right  circular  cone,  to  be  placed  with  apex  at  the  origin,  and 
axis  along  the  2-axis.  Also  suppose  a  perfectly  conducting  wire  to  coincide 
with  the  cone  axis.  Then  our  solution  corresponds  to  electromagnetic  waves 
traveling  along  the  wire  in  the  space  between  it  and  the  cone  surface. 

In  this  example,  as  in  the  previous  one,  the  direction  of  either  of  the 
two  symmetric  isometric  coordinates  might  have  been  used  to  correspond 
to  the  direction  of  the  electric  vector,  the  other  to  the  magnetic  vector,  with 
a  corresponding  change  in  interpretation  of  boundary  conditions. 

§7.  The  authors  wish  to  express  their  appreciation  to  Professors 
W.  Magnus  and  F.  Reiche  for  carefully  reading  and  checking  the  results 
contained  herein  and  for  giving  them  many  useful  suggestions. 
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On  the  Non-Existence  of 
Continuous  Transonic  Flows  Past  Profiles  I* 

CATHLEEN  S.  MORAWETZ 


In  this  paper  we  shall  show  that  the  perturbation  problem  belonging  to 
a  two-dimensional  steady  transonic  flow  past  an  obstacle  is  not  correctly 
posed.  This  theorem  has  been  proposed  in  various  forms1  [1,  2,  3,  4]  as  an 
explanation  for  the  breakdown  of  continuous  transonic  flow.  No  rigorous 
proof  has  been  given  before.  The  statement  of  the  theorem  may  be  found 
as  Conjectured  Theorem  C  in  [1],  The  basic  plausibility  arguments  were 
developed  by  Busemann  [2],  Frankl  [3]  and  Guderley  [4], 

The  perturbation  problem  may  be  described  roughly  as  follows.  Sup¬ 
pose  for  some  Mach  number  at  infinity,  <  1,  there  is  a  steady  conti¬ 
nuous  symmetric  transonic  flow,  past  a  profile  given  by  y  =  4;  Y  (x),  with 
continuously  differentiable  potential  <p  and  stream  function  ip.  Analytic 
expressions  describing  such  flows  are  known,  see  for  example  the  work  of 
Lighthill,  Cragg  and  Goldstein,  or  Tomatika  and  Tamada.  <p  and  ip  form  a 
solution  of  a  boundary  value  problem  for  a  pair  of  nonlinear  elliptic-hyper¬ 
bolic  equations  of  first  order  and  ip  vanishes  on  y  =  ±  Y(a;). 

Consider  any  neighboring  profile  P-f-  6P given  by  y  =  ±  ( Y(x )  +  SY(z)). 
If  the  boundary  value  problem  described  above  is  properly  set  there  should 
exist,  for  general  variations  of  the  profile  <3Y (z),  a  flow  with  stream  function 
ip  +  dip  and  potential  <p  +  <595,  having  the  same  Mach  number  and  flow 
direction  at  infinity,  satisfying  ip  -f-  dip  =  0  on  y  =  ±  (Y (z)  -f  <5Y (*))  and 
for  which  Sip,  d<p  and  their  derivatives  are  small  of  the  order  of  6Y.  Such  is 
the  case  if  the  original  flow  is  everywhere  subsonic  provided  SY  is  sufficiently 
smooth.  For  mixed  flows  we  know  that  we  cannot  vary  the  profile  by  putting 
a  concave  piece  into  the  supersonic  region,  see  for  example  [5,  p.  294].  Nor 
can  we  have  a  flat  piece  [6]  or  a  piece  with  zero  curvature  [7]  in  the  super¬ 
sonic  profile.  These  conditions  limit  the  admissible  functions  SY  (z)  but  we 
shall  show  that  much  more  is  true. 

We  shall  prove  within  the  framework  of  the  perturbation  theory  that 


♦The  results  of  this  paper  were  obtained  under  the  sponsorship  ot  the  U.S.  Army, 
Office  of  Ordnance  Research,  Contract  DA-3O-O09-ORD-835  and  were  announced  in  Abstract 
433,  Bulletin  of  the  American  Mathematical  Society,  Volume  60,  No.  4,  July  1954. 

'See  [1]  for  a  more  complete  bibliography. 
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for  convex  symmetric  bodies®  the  variation  of  the  profile  can  he  prescribed  ad 
most  outside  some  finite  arc  containing  the  point  of  maximum  velocity.  Except 
for  an  arbitrary  constant  8Y  ( x )  on  the  arc  is  determined  by  the  variation 
6Y(x)  outside  the  arc. 

The  proof  is  based  on  a  uniqueness  theorem.  d<p  and  dip  satisfy  a  pair  of 
mixed  elliptic-hyperbolic  equations  of  first  order.  If  terms  of  order  greater 
than  one  in  d<pt ,  8<py ,  d ipm ,  8y>t  are  neglected,  these  equations  are  linear  with 
coefficients  depending  only  on  <px,  q>v.  Furthermore  d<px,  8<py  vanish  at 
infinity  and  the  value  of  dip  is  determined  on  the  profile  by  8Y  (x).  We  prove 
that,  except  for  two  constants,  dip  and  dip  are  uniquely  determined  every¬ 
where  by  the  value  of  <5F ( x )  on  part  of  the  profile,  so  that  dY (x)  is  deter¬ 
mined  on  the  rest  of  the  profile. 

We  would  have  a  rigorous  proof  that  continuous  transonic  flows  past 
profiles  do  not  exist  in  general  if  the  same  uniqueness  theorem  could  be  prov¬ 
ed  without  neglecting  the  higher  order  terms  in  dip  and  dip,  see  Conjectured 
Theorem  A  [1], 

In  what  follows  we  formulate  the  problem  not  in  terms  of  the  variations 
dip  and  dip  but,  completely  equivalently,  by  considering  a  set  of  flows,  with 
the  same  Mach  number  and  flow  direction  at  infinity,  past  profiles  depending 
differentiably  on  a  parameter  r  and  also  on  an  arbitrary  function.  The 
equivalent  non-existence  theorem  is  that  the  velocities  can  not  depend 
differentiably  on  r  except  for  a  particular  choice  of  the  arbitrary  function. 
This  demonstrates  that  there  do  not  exist  continuous  transonic  flows  with 
the  smoothness  properties  one  might  have  expected  from  a  study  of  the 
behavior  of  subsonic  flows. 

I  am  very  much  indebted  to  L.  Bers  for  his  valuable  criticism  and  many 
suggestions.  In  particular,  the  treatment  of  the  point  at  infinity  is  due  to  him. 

1.  Formulation  of  the  Problem 

Consider  a  given  smooth,  two-dimensional,  steady,  irrotational,  com¬ 
pressible  flow  JF0 ,  with  potential  <p(x,  y),  stream  function  ip(x,  y),  and  com¬ 
plex  velocity  u  —  iv,  past  a  smooth  convex  profile  P0 .  The  free  stream  Mach 
number  M0  is  less  than  1  but  there  is  a  finite  supersonic  region  next  to  the 
profile.  The  flow  «^0  will  be  referred  to  as  the  unperturbed  flow.  The  profile 
and  flow  are  symmetric  and  we  shall  consider  only  the  upper  half  of  the 
x,  y-plane.  The  flow  JF0  is  subsonic  except  in  some  finite  neighborhood  of  the 


•Manwell  [8]  has  shown  for  mixed  flow  around  a  cylinder  with  zero  velocity  at  infinity 
that  there  are  infinitely  many  variations  of  the  cylinder  which  admit  continuous  solutions  of 
the  perturbation  problem.  But  they  are  such  that  this  boundary  value  problem  is  also  badly 
set. 
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profile  where  it  is  supersonic.  The  curve  dividing  the  subsonic  and  the  super¬ 
sonic  region  is  called  the  sonic  curve.  It  intersects  the  profile  at  Xt  and  X2 . 
Let  S  be  the  segment  of  the  profile  cut  out  by  two  Mach  lines,  C1  and  C4 , 
issuing  from  a  point  X  on  the  sonic  curve  (see  Figure  1).  S  contains  the  point 
of  maximum  velocity  on  the  profile.  The  end  points  of  S'  are  denoted  by 
X3  and  Xt. 


Let  &(x,  y,  t)  and  W(x,  y,  r)  be  the  velocity  potential  and  stream 
function  of  a  family  of  irrotational  steady  symmetric  continuous  flows  , 
depending  differentiably  on  r,  with  the  same  free  stream  Mach  number  as  , 
around  a  set  of  profiles  P(r)  that  coincide  with  P0  except  on  S.  For  r  =  0 

the  flow  is  SFq.  The  functions  co(x,  y)  =  —  <P(z,  y,  r)|T_o  and  co*(x,  y)  = 
d  dr 

—  ^{x,  y,  x)!^  are  the  perturbation  potential  and  stream  function, 

OX 

respectively,  and  are  defined  in  the  domain  <?(P)  exterior  to  P0.  It  is 
assumed  that  &x ,  ,  'Px ,  Wy  are  differentiable  with  respect  to  r  and  that 

r-differentiation  commutes  with  differentiation  in  the  x  and  y  directions. 
The  potential  0  and  stream  function  W  are  related  by  the  equations 


(1) 


The  density  g  is  a  given  function  of  velocity,  g  =  p(|wj2),  |  wj2  =  0%  +  0%. 
The  functions  <o  and  oj*  satisfy  (1)  differentiated  with  respect  to  r  at 
r  =  0;  that  is,  in  ^(P), 


(2) 


—  OJ*  =  q(wx  +  2mpx(<px(0x  +  <pycoy)), 
o>*  =  +  2nupy(<pxmx  +  ?/»,)), 


where  m  =  -^-  (dg/d  |af)  and  o  =  e(M*)  |T-0. 

The  function  m  is  given  by  Bernoulli’s  equation,  qdq  +  g-h*  dg  =  0, 
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where  q  is  the  velocity  of  the  unperturbed  flow,  that  is, 


(3) 


1 


m  =  — 


2ca 


and  c  is  the  speed  of  sound.  We  assume  that 

<4> 

which  is  true  for  all  gases  whose  pressure-density  relation  satisfies 

>0. 

The  boundary  condition  is  W  —  0  on  P  and  for  y  —  0,  and  hence 


(5) 


to*  =  0,  =  0, 


on  P0  —  S  and  on  the  ®-axis. 

On  the  remainder  of  the  profile  we  set  y  —  Y(x,  x)  =  Y0{x)  J-  rY^a;) 
+  •  ■  •  where  Yr{x)  is  a  function  that  vanishes  at  the  end  points  of  S,  Xa 
and  X4  .  The  boundary  condition  on  S  is  W (a;,  y,  r)  =  0.  Differentiating 
this  condition  with  respect  to  r  at  r  =  0  yields 

(6)  to*  +  y)vYj  =  0  on  S. 

The  potential  <p  of  the  unperturbed  flow  is  assumed  to  be  smooth: 
First  the  third  derivatives  of  <p  all  exist  in  <o(P)  and  the  second  derivatives 
are  Holder  continuous  on  the  boundary  except  at  the  nose  and  trailing 
edge.  Secondly,  the  upper  half  of  the  domain  outside  P  is  mapped,  con¬ 
tinuously  differentiably  except  at  the  nose,  the  tail  and  at  oo,  in  a  one-to- 
one  way  into  the  u,  v-plane  where  u  =  <px(x,  y),  v  =  <py(x,  y)  and  the 
Jacobian 

(7)  xuy,  —  xuyv  ^  0  except  at  oo. 

We  also  assume  that  on  the  profile  in  the  supersonic  region 


(8) 


—  is  monotonic, .  , 
dd  ‘ 


q 2  ^  c2. 


The  functions  <p  corresponding  to  transonic  flows  computed  by  the 
hodograph  method  have  these  properties. 

At  infinity  we  have 

(9)  4>x  -  qn  >  0 

where  q9  is  the  speed  at  infinity. 

Since  qm  is  assumed  to  be  independent  of  r,  we  have 

(10)  to,  —  icoy  — >■  0  at  oo. 

Under  these  conditions,  we  can  prove  the  following  lemma. 
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lemma  1.  If  (o  and  <w*  satisfy  (2),  (5),  (10)  in  <?(P)  and  have  continuous 
derivatives  for  x2  +  y2  >  R2  then 

(11)  (1  —  Ml)-i(ou  —  iwv  =  b(u  -  qx  —  iv)~±(l  -f  o(l)) 

where  b  is  a  constant. 

If  we  assume  that  the  flow  is  incompressible  for  speeds  close  to  qx  , 
Lemma  1  is  easily  proved.  First  0X  —  i&y  —  qx  is  analytic  in  x  +  iy  for 
large  x2  +  y2.  Therefore,  by  the  condition  of  symmetry  and  the  single¬ 
valuedness  of  the  stream  function  we  find 


(12) 


-  loo 


,  _fiW_  , 

(a:  +  iy)2  \  x  +  iy 


where  a0  ,  a1 ,  •  ■  •  are  differentiable  with  respect  to  r.  Differentiating  (12) 
at  r  =  0,  which  is  easily  justified,  we  have 


(13) 


ft)x  —  icoy  = 


«Or(0) 

(x  +  iy)2 


(1  +  0((x  +  iy)-1) . 


Inverting  (12)  for  r  =  0  and  differentiating  with  respect  to«  =  it  —  qm, 
we  find 

d(x  +  iy)  =  «o;,(0)  ,,  , 

d(u  —  iv  —  qx)  2(w  —  iv  —  q)*1* 

Therefore 


ft>„ 


,tt>p  =  2(u  —  ~h>  —  it  1  F  °('?)  ■ 


2(w 


?•)* 


This  yields  (11)  since  =  0. 

For  general  compressible  flow  Lemma  1  is  established  in  Appendix  1, 
by  using  the  corresponding  knowledge  for  compressible  flows  given  by  Bers 

in  t11^ 

The  flows  are  assumed  to  have  continuous  velocities  and  piecewise 
continuous  accelerations,  except  at  X0  ,  Xx  and  X2 ,  which  are  differentiable 
with  respect  to  x  at  r  =  0;  that  is, 


(14)  to  has  continuous  derivatives  and  piecewise  continuous  second 
derivatives  in  <?(P)  —  A0  —  X1  —  X2  . 

It  follows  from  (14)  and  (2)  that  the  same  condition  holds  for  co*.  It 
is  natural  to  require  that  the  perturbation  velocities  a>x  ,  co,  remain  bounded 
at  X0 ,  Xx  and  Xt ,  but  actually  our  argument  holds  if  we  merely  assume 

(15)  cox  and  co,  are  o{d~*)  where  d  is  distance  from  A0  ,  Xx  or  A,  .  x  is 
1  for  X0  and  1/8  for  and  Xt . 
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It  has  been  conjectured  that  the  perturbation  problem  might  be  cor¬ 
rectly  posed  if  a  sufficiently  strong  singularity  were  permitted  at  Xx  or  X2  . 
Condition  (10)  indicates  what  least  order  of  singularity  is  necessary. 

Let  us  call  a  pair  of  functions  (co,  co*)  a  solution  of  Problem  I  if  they 
satisfy  (2),  (5),  (6),  (10),  (14),  (15).  From  Lemma  1,  co  also  satisfies  (11). 
The  statement  that  the  perturbation  problem  is  not  correctly  posed  means 
that  there  does  not  exist  a  solution  of  Problem  I  except  for  particular  choices 
of  the  functions  Yx(x).  We  shall  prove  in  fact  the  following  theorem. 

theorem  I.  There  is,  except  for  a  constant  factor,  at  most  one  junction 
Yx(x)  such  that  a  solution  of  Problem  I  exists. 

The  function  co  will  be  called  a  solution  of  Problem  II  if  it  satisfies 
(2),  (5),  (14),  (15)  and 

(IS)*  cou  —  icov  —  0((u  —  qx  —  iv)~Vt )  as  x2  +  y2  ->  oo. 

Note  that  no  boundary  condition  is  imposed  on  5. 

We  shall  prove: 

theorem3  II.  The  only  solution  of  Problem  II  is  co  =  constant,  co*  ==  0. 
Theorem  II  implies  Theorem  I. 

In  fact,  suppose  Problem  I  has  two  solutions  corresponding  to 
Yx  =  Yn{x),  Yx  =  yi2(a:).  Let  the  two  solutions  be  (a^ ,  co*)  and 
(co2,  co*)  respectively  and  let  the  constants  b  in  (11)  be  bx  and  b2  respec¬ 
tively.  We  may  assume  that  cox  =  co2  =  0  at  oo.  If  bx  or  b2  is  zero,  then 
(«! ,  co* )  or  (co2 ,  co*)  is  a  solution  of  Problem  II.  If  bx  and  b2  are  different 
from  zero,  consider  the  pair  of  functions  (b^  —  bxco2 ,  b2co*  —  bxw* ) .  They 
satisfy  (2),  (5),  (14)  and  (15)  and  also  (13)*  at  oo.  Hence  they  form  a 
solution  of  Problem  II.  By  Theorem  II,  b2u>*  =  bxco*  and  by  (6),  ^2^11  — 

MV 

The  main  theorem,  Theorem  II,  is  a  uniqueness  theorem  for  a  system 
of  equations  (2)  of  mixed  type  with  a  homogeneous  boundary  condition  (5). 


2.  Problem  II  in  a  Modified  Hodograph  Plane 

Theorem  II  is  proved  in  the  Q,  d-plane  where  6  and  a  are  certain  simple 
functions  of  the  unperturbed  velocities. 

(a)  Transformation  of  the  differential  equations. 

First  we  consider  Problem  II  in  the  hodograph  plane  of  .  We 
introduce  the  unperturbed  velocities 


«  =  y ). 

_  V  =  <p„(x,  y), 

•The  relation  between  Theorem  II  and  the  non-existence  statement  was  pointed  out 
by  L.  Bers. 
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(17) 


(18) 


as  independent  variables.  Equations  (2)  become 

—  —  0J* %y  =  Q0Ju(  (1  +  2mui)<Pxx  +  2 muv<pxy) 

+  e«D„((l  +  2  mu2)q>xy  +  2  muv<pvv), 

Kfxx  +  ™t<Pxv  =  Q“>u{(  1  +  2mv2)<Pzv  +  2muv<pxx) 

+  e«>»((l  +  2mv2)<pyy  +  2  muv<pxy) 

or,  by  solving  for  co*  and  co*  , 

(<Pxx(Pvy—<plv)a>Z  =  e^u[<Pxy{  ( 1  +  2mu2)<pxx+  2muv<pxy  +  (1  +  2  mv2)(pyy} 
+ 2muv<pxx<pyy] + Q(tiv[<pvv{4muv<pxv+  ( 1 + 2  mv2)<pyy} 

+  (1+2  mu2)(ply], 

{'PxxVvv-'Plv)^  =  -  &»«[?**{  (!+ 2mu2)<pxx+imuv<pxy}+  (1+ 2tnv2)<p\y\ 
—  QMv  [<Pxv{il  +  '2mu*)(pzx+  2muv<pxy+  (1  +  2mv2)<pyy} 
+2muv<pyy<pxx]. 

But  for  t  =  0  we  have  from  (1)  differentiated  with  respect  to  x  at  t  =  0, 
(1  +  2mu2)tpxx  +  \muv<pxy  +  (1  +  2  mv2)cpyy  =  0. 

Equation  (18)  reduces  to 

(<Pxx<Pvy  ~  u  =  QWxx<Pvv  ~  <Ply)[2mUV(Ou  -  (1  +  2mM>J, 

{<Pxx<Pyy  ~  9&K*  =  e(<Pxx<Pyy  ~  9&)[(!  +  2 ~  2mUV(Ov] . 

Since  the  Jacobian  <Pxx<py,  —  <p\y  =  uxvy  —  uyvx  does  not  vanish  except  at 
oo  by  (7)  we  find  using  (3) 

0_1C2CO*  =  —  UV(Du  —  (c2  —  U2)d)v, 
g_1c2to*  -  (c2  —  v2)cou  +  UV(Ov  . 

If  we  introduce  as  independent  variables  the  angle 

v 


(19) 


(20) 


6  =  tan-1- 


u 


and  the  speed  q  =  V u2  +  v2,  equations  (19)  can  be  reduced  to 


(21) 


c2q(o*  =  q(c 2  —  q2)coe , 

«>*  =  —  • 


These  equations  have  been  derived  by  Manwell  in  [8], 

Note  that  the  equations  satisfied  by  the  perturbation  stream  function 
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and  potential  are  very  similar  to  the  equations  for  the  stream  function  y> 
and  potential  cp  as  functions  of  6  and  q: 

c2m  =  -  9*)v»  - 

<Pe  =  -  e_1?w  ■ 


The  perturbation  stream  function  replaces  the  potential,  the  perturba¬ 
tion  potential  replaces  the  stream  function  and  o~x  replaces  g. 

Next  we  set 


(22) 


a=-fdl=-r 

(c*  -  q*) 


dq 


exp 


K(a)  —  g 


2  _  q2 


{/>}■ 

\-2\'}dq\ 


where  c*  is  the  value  of  q  for  which  q  =  0.  We  assume  without  loss  of 
generality  that  the  stagnation  density  is  1.  Note  that  for  small  q, 
a  ~  —  log  q. 

Equations  (21)  then  reduce  to 


(23) 


(o*  =  —  Km9  , 
CJ*  =  0>„  ■ 


Note  that  K(a )  has  the  properties 

oK(o)  S?  0 
dK 

(24)  —  >  0  for  a  ^  0. 

da 


which  we  shall  need  later.  For,  by  (22),  a  and  K(a)  vanish  only  for  q  =  c*. 
Also 


and 


and 


dK 

da 


dc 2  dc* 

dq 2  d{\ w\)2  I=0 


by  (4) 


(c2  -  q*)/K  >0  by  (22). 


From  (23)  differentiated  with  respect  to  0  we  obtain  an  equation  for  to, 

(25)  K(om  4-  (oaa  =  0. 

The  differentiation  can  be  easily  justified  from  our  assumptions. 

(b)  Transformation  of  the  boundary  condition. 


I 


* 
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From  (5)  and  (23)  we  have 

(26)  —  dot*  =  Katedo  —  coadd  =  0 

on  P0  —  S  and  the  x-axis  and  hence  also  on  77  —  Z  and  A,  the  images  of 
P0  —  S  and  the  x-axis  in  the  0,  a-plane. 

Note  that  the  function  w  satisfies  an  equation  of  mixed  type  similar  to 
the  Tricomi  equation,  with  a  homogeneous  boundary  condition  on  an  open 
arc. 

Alternatively  from  (23),  at*  is  a  solution  of  the  equation 

(27)  K<4e  +  <°*a  ~  It  it™*  =  ° 

and  the  boundary  condition  (5)  reduces  to 

(28)  at*  =  0  on  PI  —  S  and  A. 

The  uniqueness  of  solutions  of  (23)  which  are  constant  on  such  open 
arcs  as  II  —  £  has  been  studied  in  [9,  10]  but  the  uniqueness  theorem  for 
(23)  and  (26)  or  (27)  and  (28)  has  not  been  investigated. 

(c)  The  mapping  in  the  0,  a-plane. 

The  image  of  the  exterior  of  half  of  the  profile  P0  ,  r?(P0),  and  the 
x-axis  in  the  6,  a-plane  is  indicated  in  Figure  2.  It  consists  of  (a)  Q  which 


is  bounded  by  II,  the  image  of  P0  —  S  and  A  the  image  of  the  x-axis,  and 
the  images  and  y2  of  Cx  and  C2 .  (b)  Q*  which  is  bounded  by  E,  the  image 
of  5  and  y1  and  yt . 
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The  nose  and  the  trailing  edge  of  the  profile  P0  are  stagnation  points 
and  therefore  correspond  to  a  =  +  oo.  Suppose  the  half  angle  at  the  nose 
is  eN  and  at  the  trailing  edge  et.  Since  the  profile  is  convex  the  domain 
Q  -f-  Q*  lies  between  0  =  eTII  and  0  —  —  eNIJ.  We  denote  by  IJT  the 

point  <7  =  +  oo,  0  rgj  0  eT  and  by  IJN  the  point  a  =  +  oo,  —  sN  0  0. 

The  point  II x  ,  0  —  0,  a  —  ax  —  a{qai)  is  the  image  of  the  point  at  infinity 
in  the  physical  plane.  The  cr-axis  for  a  >  aK  corresponds  to  the  a-axis 

outside  P0  .  The  axis  a  =  0  is  the  sonic  curve.  The  points  X0 ,  X1 ,  X2  , 

X3  ,  Xt  are  mapped  into  770 , 77j  ,  /72  ,  773  ,  774 . 

The  arc  77  satisfies  the  important  condition 

(29)  Kda 2  +  d6 2  ^  0. 

This  expression  can  be  negative  only  for  a  <  0.  But  by  (5),  if  Kda 2  +  dO2 
vanishes  on  77,  then  Kyfy  -f  yl  also  vanishes,  hence  by  the  Chaplygin 
equations  and  equation  (22), 

_  _  d  (a,  y)  __  d(<p,  y>)  _  d(0,  a)  /  d{(p,  y) 

XuVv  X'>Vu  3{u,v )  d(6,  a)  d(u,v)/  d{x,  y) 

=  Q~3{K^e  +  wDifl  +  =  0 

d(u,  v) 

since  d(6,  o)/d(u,  v)  is  bounded.  But  this  violates  condition  (7). 

Since  and  y2  are  characteristics  the  slopes  of  y1  and  y2  are  given  by 
the  equations 

dd  -  V—  K  da  on  ylt 

(30)  - - 

dd  —  —  V  —  K  da  on  y2. 

Since  the  profile  is  convex  we  have  dd  >  0  on  77  —  27. 

(d)  Behavior  of  derivatives. 

By  (9)  the  derivatives  a>x  and  cov  are  continuous  in  <?(P)  except  at 
X0  ,  Xx ,  X2  and  at  infinity  and  therefore  as  functions  of  0,  a  they  are  con¬ 
tinuous  except  at  770  ,  77j ,  772 , 77v  ,  77r  and  77*  .  The  functions  xe ,  xa  , 
ye ,  ya  are  also  continuous  except  at  IIN  ,  77r  and  Ux  since  the  mappings 
(x,  y)  (u,  v)  ->  (0,  a)  are  continuously  differentiable  except  at  these 
points.  We  have  then 

(31)  cog  and  oja  are  continuous  except  at  770 , 77j ,  772  ,  II N  ,  77r  and  77^  . 

We  shall  need  a  condition  on  e oa  at  770 , 773  and  774  .  By  (15)  cox  and  coy 
are  o{dX()  where  d  is  the  distance  from  X0 ,  Xx  and  X?  respectively.  The 
functions  xu  ,  xv  ,  yu ,  yv  are  boundedat  these  points.  H  ence  a>u  and  wv  and 
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therefore,  by  (20)  and  (22),  coe  and  ma  are  o(d”).  Since  the  mappings  (x,  y) 
^  (u,  v)  ^  (0,  a)  are  differentiable  at  XQ ,  X±  ,  and  X2  it  follows  that 

(32)  are  °^ai  +  ~  0«)2)~“'/2)  where  **  =  1 > 

6i  =  0O  at  II0,  x{  =  1/8,  0,  =  01(  02  at  JJX  ,  II2. 

We  next  determine  the  behavior  of  the  derivatives  co9  and  co0  at  1JN 
and  IJT .  These  two  points  correspond  to  the  nose  and  the  trailing  edge  where 
the  velocity  is  always  zero.  (By  the  convexity  condition  on  P0  the  angle  at 
the  trailing  edge  is  finite).  It  can  be  shown  by  the  use  of  the  asymptotic 
expansions  at  the  nose  and  tail  for  subsonic  flow  given  in  [11]  (see  Appendix  1 
for  details)  that 

(33)  ojeea,  (oaea  are  bounded. 

We  verify  (33)  here,  assuming  for  simplicity  that  the  flow  is  incom¬ 
pressible  near  zero  velocity.  Let  z  =  x  iy  and  zT  be  the  value  of  z  at  the 
trailing  edge  and  tn  the  half-angle  there.  Map  a  neighborhood  of  the 
trailing  edge  onto  a  vertical  half  disc,  dte  £  >  0,  |  f  |  <  1  in  such  a  way  that 
the  trailing  edge  goes  into  the  origin,  and  a  part  of  the  profile  into  the 
imaginary  axis.  Then  £  ~  const,  (z  —  zT)ll(l~E) .  As  a  function  of  £  =  f  +  it}, 
0  +  iW  is  analytic  at  £  =  0  and  +  iW  —  aj(r)£(l  +  a](r)£  +•••)•  It 
may  be  easily  justified  that  we  may  differentiate  with  respect  to  r  and 
obtain  co  -f-  io>*  —  aOr(0)£(l  +  0(C)).  Differentiating  with  respect  to  z  we 
obtain 

(34)  a>x  —  ic »„  ~  const,  (z  —  zr)t/(1_f) 
and  from  the  expression  for  0  -f-  iW  at  r  =  0, 

(35)  <px  —  i<py  =  u  —  iv  ~  const,  (z  —  zr)e/(1_€). 

Furthermore  differentiating  (35)  we  find  that  xu ,  xv  ,  yu  ,  yv  are 
0((m  —  jv)1/<€-2)).  Thus  cou  and  a>v  are  0((u  —  fw)(1“e)/e).  By  (20)  and  (22)  it 
follows  that  as  o  go,  (33)  holds  for  Tlr  since  e  ^  Similarly  we  can 
establish  (33)  for  IJN. 

The  behavior  of  w9  and  coa  at  77„  is  derived  from  (13*)  and  (22), 

(36)  coe,  v)a  are  o(|0  +  i(a  —  OI-1,4)  at  . 

Problem  II  has  been  transformed  into  Problem  III  in  the  0,  0-plane; 
a  function  w(6,  a)  is  a  solution  of  Problem  III  if  it  satisfies  (25),  (26),  (31), 
(32),  (33),  (36)  in  the  domain  Q  +  Q*.  Theorem  II  is  then  equivalent  to 
showing  that  any  solution  of  Problem  III  is  constant  in  Q  +  Q*.  This 
follows  from  Theorems  III  and  IV. 
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theorem  III.  The  solution  of  Problem  III  is  constant  in  SI. 

theorem  IV.  Suppose  co  satisfies  (25),  (31)  and  (32)  in  a  domain  R 
bounded  by  four  characteristics  yx ,  y2  ,  y3  ,  yt  where  yx  and  yi  correspond  to 
greater  values  of  a.  If  co  is  prescribed  on  yx ,  y2  then  co  is  uniquely  determined 
in  R.  Or  equivalently  if  co  =  0  on  yx  +  y2 ,  co  =  0  in  R. 

Theorem  IV  which  is  a  slight  generalization  of  a  standard  theorem  is 
proved  in  Appendix  2.  The  domain  R  is  taken  to  be  the  domain  bounded 
by  ylt  y2  and  the  characteristics  through  their  end  points.  Then  co  =  0 
in  SI  -f-  R  I)  SI  -j-  SI*. 


3.  Proof  of  Theorem  III 

Theorem  III  is  proved  by  the  A  FC-method  of  Friedrichs.  Let  to  be  a 
solution  of  Problem  III  and  consider  the  integral 

(37)  Is  =  JJ  ( Bme  +  C(oa)  (KcoM  +  cooa)  dOdo  =  0 

Sid 

where  SIs  is  a  subdomain  of  Q  and  B  and  C  are  functions  of  0  and  a  to  be 
chosen  later.  The  domain  Qs  is  bounded  by  the  curve  36 6  which  consists  of 
the  two  characteristics  yx  and  y2  and  the  boundary  curves  77  —  E  and  A, 
except  in  the  neighborhood  of  770  ,  TIX  ,  772  and  77M  where  it  consists  of  four 
curves  7?w  ,  RXi ,  7?M  and  RM  which  exclude  770 , 771 , 772  and  IIX  and  whose 
length  shrinks  to  zero  as  <5  ->  0,  see  Figure  3.  Rw  ,  Rxs  ,  R^  and  7?M  will  be 


chosen  more  specifically  later.  We  denote  the  portion  of  36 d  which  coincides 
with  II  —  yx,  yt,  by  nd  ,  yu  ,  yu ,  respectively. 

We  shall  show  how  to  choose  B  and  C  in  such  a  way  that  unless  co  is 
constant 
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lim  7.  >  0. 


<s— >o 


This  contradiction  will  prove  the  theorem. 

First  B  and  C  are  to  be  chosen  so  that  we  may  integrate  (37)  by  parts 
and  obtain 

h  =  J  J  kb»  +  {KC)^~  +  KC°) 

+4«5(*»-c-)}‘ 


■  d6da 


(38) 


or 

(39) 


-f  (jl  —  [K(o\da  —  2a)e(oadd  —  u?„da) 
J  2 

Q 

-f - {Ko)\dO  -f-  2Ka>ea>ado  —  Q%dO) 

2 


Id  =  hs  +  I‘u+  7m  +  Iu 


where  Ilt  ,  /M  ,  /M  ,  Iu  are  the  area  and  line  integrals  over  the  elliptic  and 
hyperbolic  regions  respectively.  Thus, 


IU=  /  J  {j«5(  -KBe+  (KC)a)  -(0,)(oa(Ba+KC6)  +  ^{B.-C^dBda, 


(40) 


%  a^O 


C  B  C 

Im—  <p  -—(Kajldo—2cD9(oadQ—coldo)Jr  —  (K(oldOJr2K(oe(uado—cold6). 

j  z  z 

3»i, 

By  (26)  we  replace  wadQ  by  Koj6da  in  the  first  parenthesis  and  Ko>eda  by 
( oad0  in  the  second,  and  obtain 

7W  =  ^  —(K<ol+a>l)(Cdd  —  Bda)+  Y  ■--(Kco2edo~-2u)ea}ad0--oyLoda) 

(41 )  n6,a^0  '  0  RiS 

c 

+  —  {Km\dB  -j-  2Ko)$o>ado  —  a>%dd). 
z 

The  integrals  in  the  hyperbolic  region  are, 

4,=  J7  {~o>l(~KBe+  (KC)a)  —wewa(Ba+KC9)  + 


(42) 


C  B  C 

Ii&—  r  —{K<ti\d<j—2(x>e<oadQ—o?0do)-\ — (Ka%dQ-]-2K(Oga)ada—a%dQ) 
J  2  2 


* 

i 
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which,  by  (26)  and  (30),  may  be  written  as 

hi=  J  —  (—Bdo—CdO)  ((o„+ V—  Kcoe)2+  J  -i  (— Bda  —Cdd)  (coa—V - Kcoe )2 

Ylt  7  m 

(4s)  +  /  ^(,+K(^'hcdo-B^+  f 

R06+R16+R26 » 

c 

—  2ci)ga)ad6 — a%d8 )  -| —  {Kca\d6  -(-  2K(oea>„do — a%dd). 

2 


We  note  that  1  -f-  K 


ldoY‘ 
\  Id) 


> 


0  by  (29). 


By  (31),  the  integrals  will  exist  for  all  <5  ^  0  if  B  and  C  are  continuous 
in  the  closure  of  Qd  and  are  bounded  at  IJN  and  77r  by  (33). 

We  shall  show  that  functions  B  and  C  exist  such  that  Iu  and  I.M  are 
positive  and  bounded  waay  from  zero  uniformly  in  d  and  such  that  lim 

and  lim  Iu  are  non-negative.  It  follows  then  that  o>e  =  wa  =  0  in  Q.  This 

A— >0 

will  prove  Theorem  III. 

The  integrals  Iu  and  Iu  are  positive  definite  and  lim  Iu  ,  lim  I1S  are 


non-negative  if 

i->o  «-*o 

-  KBe  +  (KC)a)  >  0 

in  Q,  a  =4  0, 

B„  -f-  KCe  =  0, 

(44) 

Be  Ca  >  0 

in  Q,  a  ^  0, 

K(Cd6  -  Bda)  ^  0 

on  IId  +  A, 

—  Cdd  —  Bda  Y  0 

on  y1  and  yt 

and  if 

(45) 

]im  7J  =  lim  Ir=  0 

<>->•0  <S->0 

where  / J  ,  I RiS  ,  are  the  line  integrals  along  Ra  for  a  7>  0,  i  =  0,  1,  2,  3, 
and  a  0  respectively,  i  =  0,  1„2. 

In  the  lower  half  plane,  a  fYt  0,  we  set 


(46) 


B  =  B{9,  0),  a  <  0, 
C  =  0. 


Then  conditions  (44)  for  a  <  0  will  be  satisfied  if 
Bt{6,  0)  >  0  on  01  £  d  ^  02  , 
(47)  5(0,  0)^0  for  0O  ^  0  ^  02  , 

5(0,  0)^0  for  Bl  ^  Q  ^  0O  , 
where  0,  =  0(17,),  i  =  0,  1,  2. 
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For  the  purpose  of  constructing  B  and  C  for  a  S:  0  we  introduce  an 
auxiliary  plane4  into  which  we  map  Si,  a  ^  0.  We  set 


(48)  fi  =jWK(o)  do. 

Then  the  upper  half  of  the  d,  u-plane  is  mapped  smoothly  into  the 
upper  half  of  the  A-plane  where  A  =  6  +  */*•  Let  770 ,  TIX ,  II2 ,  IIT  ,  77v  , 
77w  be  mapped  into  17*  ,  II*  ,  II*  ,  11}  ,  11%  ,  11*  .  Let  Ra  ,  IIt ,  Sit ,  A 
and  77  —  Z  for  a  Sg;  0  be  mapped  into  7?J  ,  II t ,  Qs  ,  II*.  Note  that  since 
on  II  —  S,  dd  >  0  for  6  7S  0  we  have 


(49) 


dd  >  0  on  77*,  6  0. 


Using  Privaloff’s  theorem  on  conjugate  functions  it  is  easy  to  show 
that  there  exists  a  function  /(A)  =  /x(0,  ft)  +  if2(d,  fi)  which  is  analytic  in 
A  =  6  +  i/i  in  the  interior  of  Si*,  continuous  in  the  closure  of  Si*  except  at 
77*  ,  77j*  ,  77*  ,  77*  ,  77*  ,  and  satisfies  the  boundary  conditions: 


(50) 


f2  =  0  on  A*, 
h  =  jg  on  n* 


where  cos  j  —  dd/ds,  sin  j  =  du/'ds  on  77*  and  by  (49)  we  may  take 
—  jr/2  ^  j  n/2,  g  is  a  given  very  smooth,  say  infinitely  differentiable, 
function  satisfying 


(51) 


OggsS  1, 

g  =  0  near  77*  ,  77*  ,  77f  ,  77*  where  dd/i  >  0. 


From  the  assumption  that  the  second  derivatives  of  <p  are  Holder  continuous 
it  can  be  easily  shown  that  jg  is  Holder  continuous.  On  fi  =  0, 


(52) 


U  =  +  for  d2  ^  d  ^  d0  , 

z 

ft  =  for  0O  ^  d  £  d2  . 

z 


At  77*  and  77*  where  A  =  0„  i  =  1,  2,  it  can  be  shown  that 

(52)  /  =  —  \  log  (A  —  df)  +  continuous  function. 

At  77*  where  A  =  0O  , 

(53)  /  =  log(A  —  do)  -j-  continuous  function. 

At  77*  ,  77*  were  fi  00, 

(54)  /  is  bounded. 


4In  the  variables  (9,  ft),  equation  (25)  takes  on  the  normal  form. 
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The  functions  B  and  C  are  chosen  for  a  ^  0  as 

(55)  VKC  -iB  =  exp  {/(A(6,  /,))}. 

We  must  show  now  (a)  that  B  and  C  are  continuous  at  a  —  0  and  all  the 
integrals  Iu  ,  1 M  ,  I M  ,  I u  exist  and  (b)  that  the  conditions  (44)  and  (45)  are 
satisfied. 

From  the  properties  of  /(A)  we  see  that  e*  is  a  regular  analytic  function 
of  the  complex  variable  A  in  Q*  with  a  square  root  singularity  at  II*  and  77*  , 
and  vanishes  to  first  order  at  77*  by  (52)  and  (53).  As  p  ->  oo  at  77?  and 
m,  e*  is  bounded  by  (54).  Also  !%e  ef  =  0  on  Jm  A  =  0. 

By  (50),  (51)  and  (52)  we  have  |/2|  <  max  (nj2,  |/|)  on  the  boundary 
and  by  (49),  |/|  <  nj 2.  Therefore  |/2|  <  n/2  and  VKC  >  0  almost 
everywhere  in  £*j5ince  by  (50),  (51)  and  (52)  it  is  not  a  constant.  Further¬ 
more  by  (52),  VKC  =  0  on  //  =  0  and  therefore  since  VKC  —  iB  is  analytic 
up  to  fi  =  0,  by  the  reflection  principle,  we  have 

VKC  =  0  .  on  fi  —  0, 

(56)  VKC  >0  in  Q*, 

(VKC)m  >0  on  ,u  =  0. 

We  next  check  (a)  and  (b). 

(a)  By  (56)  VKC  vanishes  at  least  like  /i  on  ju  =  0.  Therefore  C 
vanishes  like  a  on  a  =  0  by  (48)  and  (24).  By  (47)  then  B  and  C  are  conti¬ 
nuous  on  a  =  0  in  . 

Plainly  the  integrals  Iu  ,  Iy, ,  ,  Iu  converge  absolutely  for  all  <5^0 

since  B  and  C  are  bounded  and  continuous  in  Qt ,  d  yt  0. 

(b)  To  check  (45)  we  must  fix  Rid .  For  a  ^  0,  Rid ,  i  =  0,  1,  2,  is  the 

image  in  the  d,  <r-plane  of  a  semi-circle  7?,*  of  radius  5  in  the  0,  //-plane.  For 
o  0,  consists  of  the  segments  cut  out  by  the  lines  a  =  —  d,  |fl  —  0,.| 
=  <5,  see  Figure  3.  7?M  is  a  circle  of  radius  <5. 

First  we  consider  the  line  integral  at  77*,  , 

c  ^  c 

~  —  {Ko)\da  —  2co8cy adB  —  colder)  +  —  (KcoldO  2Kcoecoada  —  coldO). 

By  (36),  cje  and  coa  are  o(<5-*).  The  coefficients  B  and  C  are  bounded. 
Therefore  lim  IR  =  0. 
a->o 

For  7Jw  ,  i  =  0,  1,  2,  we  use  the  d,  //-plane: 
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n<‘“h 

J  RiS> 


B 


(Kafydct  —  2  coe(o„dd  —  co2do) 


+  —  ( Kco2ed6  -f  2(oe(oada  —  co2ad6) 
=  y  {VK(o2ed/i,  -  2(oeco„dd  = 


+  ■ 


Va:c 


V Kw\dQ  -f-  2coe(oadf,i - -  i0 

VX 


)' 


Now  from  (52),  (53)  and  (55)  we  have 

VKC  —  iB  =  0(A  —  0,)=* 

where  x,  =  1,  0,  =  0O  at  /7a  and  0t  =  0j ,  02  at  TIX  and  772 . 

From  (24),  dK/da  >  0,  we  see  that  .K  ~  iC'(0)er  and  by  (48),  /t  ~ 
f  (if'(0))1/«<//»,  K  ~  (f)s/3(i^'(0))J/a/r2/a.  Hence  since  w8  ,  &>a  are,  by  (32), 
o(((0  —  0,)2  +  o,2)x,/2)  they  are  o(<52*i/3)  on  I?*  .  We  set  tan  v=  ///(0  —  0.) 
on  1?,*  ,  then  VK  =  (f)l/»(/C'( 0))1/j<51/a  sin1/a  r  on  I?*  and  thus 


Ktj = »(#",+'+l)  [c™'s  >*+/;*]  - 

Substituting  the  values  of  from  above  and  xt-  from  (32) 

^  =  «U>- 

Hence 


For  the  integrals  in  the  hyperbolic  region  we  have,  using  (46), 


_  «V*V*0  —  WodQ)- 

By  (52),  (53)  and  (55),  |£(0,  0)|  =  0{(6  -  04)s‘)  and  |£(0,  0)|~ 
const.  (5*‘  on  |0  —  0f|  =  <5.  By  (32),  |<w9|  and  |cae|  are  o(d~Hi)  on  all  of 
Rm  .  Hence 

i~t)  =  0{i r2*^1)  +  o(d~2x‘)  fjOi  -  o\x<  do  =  o(r2x'+i'+1), 

and  by  (32),  lim  R ^  =  0.  Condition  (45)  has  been  proved. 

It  remains  only  to  establish  the  inequalities  (44)  and  (47).  We  trans¬ 
form  (44)  and  (47)  to  the  0,  //-plane  and  show  that 
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(57)  -  K(Be  -  (VKC)„)  +  1^VXC>0 

(58)  B„  +  (VKC)e  =  0, 

(59)  (VKCl+i 

(60)  VKC  dQ  -  Bdfi  ^  0 

(61)  Be  >  0 

(62)  B  ^  0 

R  ^  0 


in  72*. 


in  13*, 

on  77,  /l, 
on  p  —  0, 

on  p  =  0,  0t  5S  0  ^  0O, 
on  p  =  0,  0O  £=i  0  02. 


Since  VKC  —  iB  is  analytic  in  X  —  0  -)-  ip,  (57),  (58)  and  (59)  hold  by 
virtue  of  the  Cauchy-Riemann  equations,  (56)  and  (24).  Using  (55)  and 
(50),  on  77, 


VKC  dd  —  Bdp  —  exp  {Me  /}  cos  (j  —  ./m  f) 

=  exp  {Me  /}  cos  (/(  1  -  g)). 

Now  |;’(1—  g)|Si|/|  by  (51)  and  by  (49),  —  n/2  gi  /  ^  n/2,  so  that 
cos  /(I  —  g)  >  0.  On  A,  dd  =  0,  IS  =  0.  Hence  (60)  holds. 

(61)  follows  from  (56)  since  (VKC)/I  =  3?e .  Finally  by  (55),  B  — 
—  exp  {Me  /}  sin  Jm  f  and  condition  (62)  follows  from  (52). 

The  proof  of  Theorem  III  is  now  complete. 


Appendix  1 

Behavior  at  Infinity  and  at  the  Nose  and  Trailing  Edge 

In  this  appendix  the  references  labelled  B  are  to  Bers  [11], 

At  infinity.  Proof  of  Lemma  1.  First  we  shall  need  the  asymptotic 
behavior  of  the  unperturbed  velocity  at  infinity.  We  use  the  notations  of 
B  Lemma  8a.  The  distorted  velocity  w*  for  the  flow  is  an  analytic  func¬ 
tion  of  £  =  £  -j-  it]  where  t  is  a  mapping  of  a  neighborhood  of  z  =  oo  onto 
a  neighborhood  of  f  =  oo,  conformal  to  the  flow  metric.  The  distorted  velo¬ 
city  w *  is  defined  in  the  subsonic  region  in  terms  of  the  velocity  w  —  <px  — 
i<pv,  see  B  (7.7).  For  £  ->  oo  it  follows  that 

(Al)  w*  =  w*  +  —  +  -p- d - ■ 

We  may  assume  the  mappings  z  ^  £  are  symmetric  with  respect  to  the 
horizontal  axes.  Then  ,  w*a  ,  •  •  •  are  all  real. 
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As  w  qK  =  wx  >  0  in  the  notation  of  Bers  (1 1),  we  find  by  expand¬ 
ing  B  (7.7),  for  j  z  |  oo 


(A2) 


w 


*1  =  A(M  -  KIM1  +  0(1)) 


where  A  is  a  positive  constant.  Hence  from  (Al) 


(A3)  \w\  -  Kl  = 


(l  +  o(l)) 


W- 


wz 


loo 


w\ 


2® 


+  0(f-») 


Kl)- 


Expanding  arg  w  =  arg  id*  in  powers  of  £  from  (Al)  and  combining 
with  (A3)  we  obtain 

(A4)  w-wm  =  w*a  ((1  +  0(1)l&e  £-i  +  ^  Jm  £~A  +  0 (£-»). 

\  A  WZ.  / 


The  mapping  from  the  physical  plane  into  the  £-plane  is  conformal  with 
respect  to  the  flow  metric: 

(A5)  d£2  -)-  drf  =  (c2  —  u2)  dx2  +  2 uv  dx  dy  +  (c2  —  v2)  dy2. 


Therefore  as  £  ->  oo,  u  ->  wx  ,  »->0we  may  choose  the  mapping  so  that 
(A6)  2/f  ->  0  as  z  ->  co 

vl  -  Ml 

and  thus 

(A7)  £  =  (*  +  iV  1  -  Mly)  (1  +  o(l)) . 

Substituting  (A7)  in  (A4)  we  have 

(]  _  ^  _  _  ... 

—  @e(x+iV  1 —Mly)  Jm(x+iV  1 —Mly) 

+o((x+iVTZZ Mly)-2)  j  +0((x+iVl— Mly)~2) 

Now 

ip  =  f  qvdx  -r-  gudy  =  —  Jm  f*  gwdz 

J  tj'  J  X'f 

is  single-valued  since  ip  =  0  on  the  profile  P0  .  Thus  Jm  jy  owdz  =  0 

where  the  path  of  integration  may  be  taken  as  the  ellipse  x2jr  (1— Ml)y2  = 
R2.  For  R  ->  oo,  we  find  by  direct  computation  that 

Iw  1  \ 

w*  — -  -) - —  —  1  =  0  and  therefore  wJ  =  0.  Thus5  from  (A8), 

K  aVi  -Ml’ 

•cf.  [12]. 
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(A9)  w  -  ww  =  0((x  +  iV  1  — 

and  in  (Al), 

(A10)  w*  =  w*  +  ~%r  +0(f-3). 


Secondly,  we  must  derive  estimates  for  xu  ,  xv  ,  yu  ,  yv  at  co  which  we 
shall  need  to  estimate  cou  —  icov  ,  we  differentiate  (A10)  with  respect  to  £ 
and  find  d£ldw*;  differentiate  B(7.7)  with  respect  to  u  and  v  and  use  (A6). 
This  shows  that  xu  ,  xv  ,  yu  ,  yv  are  0{\w  -  q„  |-s/*)  or  by  (20),  (22)  and  (48), 
xe ,  *„  ,  y, ,  y„  are  0( |A  -  XJ~1')  where 

A  =  0  +  iy,  A«,  =  *>(?«)• 

Since  a>x  —  ioov  0  it  follows  that 

(All)  coe  —  *'gv  =  0(W-,/*)- 

Finally,  to  derive  (11),  we  use  the  fact  that  coe  —  is  a  pseudo-analy¬ 
tic  function  of  the  first  kind.  For,  since  co  and  co*  satisfy  (2),  they  satisfy 
(23)  and  therefore  co  satisfies 

dK 

(A12)  %  +  %  —  0 

where  y  =  jWK(a)  da,  (48).  It  follows  by  the  similarity  principle,  see  [13], 
that  in  the  neighborhood  of  (0,  yx) 

(Al 3)  coa  -  ico „  =  exp  {s (0,  p)}W{X  —  XJ. 

Here  s  is  Holder  continuous,  real  on  0  =  0  and  vanishes  at  (0,  ym),  and  W 
is  analytic.  Since  co*  satisfies  (5)  ony  =  0  we  have  by  (26),  co„  =  0  on  0  =  0, 
y^yx.  Introducing  A*  =  (-  t(A  -  AJ)V*  =  0*  +  iy*  we  see  that 
3te  W  vanishes  on  Jm  A*  =  0  and  by  reflecting  about  Jm  A*  =  0  that  W  is 
a  singlevalued  function  of  A*,  0  <  jA*|  <  jAj|  .  Therefore  by  (All) 

CO 

(A14)  <o«  —  =  e*  2  an?*n- 

-2 


Since  s  is  continuous  and  vanishes  at  (0 ,yx), 

(A15)  co6  -  iw„  =  «_,A*-J(1  +  o(l))  +  a-, A*-X(l  +  o(l)). 
From  (23)  and  (26)  we  have 


§ 


doo*  =  —  Jm 


ioop)  (dB  +  idy)  =  0 
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where  the  path  of  integration  is  a  circle  of  radius  (5  about  (0,  fix). 

Letting  8  ->  0  we  obtain  by  (A15),  Sm  =  0.  By  (26)  and 

Jm  s  =  0  on  0  =  0  we  have  /m  a_2  =  0.  Hence  a_2  =  0. 

Therefore 

(A16)  (oe  —  ic a,,  —  +  o(l)) . 

By  (48),  (22),  (20)  and  (9), 

**»  =  %^1(1  —  M2«o)1/’(1  +  oi1))  +  "o0!1)- 
«>»  =  (oeq~1(  i  +  o(l))  +  eo„(o(l)). 

Equations  (A17)  yield  (11). 

At  the  Nose  and  Trailing  Edge.  Derivation  of  (33).  For  a  compressible 
flow  \<px  —  i<pv\  =  \w\  ~  const  \z  —  zT\el{1~e)  where  e/7  is  the  half  angle 
of  the  trailing  edge,  z  =  x  -f  iy  and  zT  is  the  value  of  z  at  the  trailing 
edge.  This  relation  is  determined  in  B  p.  478. 

To  estimate  xu  ,  xv  ,  yu  ,  yv  we  use  the  method  of  B  §8.  We  can  derive 
a  representation  of  the  form  B(8.7)  for  the  distorted  velocity 


■(f)  =  ^(i-I)'. 


Here  f  is  again  a  mapping  conformal  to  the  flow  metric,  say,  of  the  whole 
subsonic  domain  onto  the  exterior  of  the  unit  circle.  The  point  f  =  1  cor¬ 
responds  to  the  trailing  edge  and  the  function  G(C)  is  analytic  and  Holder 
continuous  at  f  =  1.  The  last  remark  implies  that  the  derivatives  of  G(C) 
are  0(|£  —  l)1-®)  for  some  a  >  0,  Hence  by  differentiating  z»*(£)  we 
estimate  dw*fd£  or  d^jdw*.  From  B  (7.7)  differentiated  we  find  dw*/dw. 
From  the  flow  metric  we  determine  ,  xn  ,  y( ,  yv .  We  find  that  xu  ,  xv  , 
yu,  yv  are  Oflwf^-*'). 

Since  wx  —  iw¥  -*■  0  at  the  stagnation  points  we  have  <o„  —  ia>v  — 

o{w’~ 2)  or,  by  (13)  and  (16),  to,  and  co,  are  o  ^exp  j*  | - 1  j  (0  +  ia)  j  j. 

Since  e  <  |  it  follows  that  cojC®,  to are  bounded,  (33). 

The  leading  edge  may  be  treated  in  the  same  way. 


Appendix  2.  Proof  of  Theorem  IV 

Let  to  vanish  on  yx  +  yx  and  to*  vanish  at  some  point  on  yx  4-  y2 .  We 

j»8,  a 

again  apply  the  A  BC-method  but  to  the  function  W  =  J  <o*do  +  to  dd. 
It  is  not  difficult  to  show  by  (23)  that  W  is  single-valued,  vanishes  on  yx 
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and  ya  and  satisfies  (25).  Let  the  intersection  of  y1  and  yt  be  at  (d0  ,  a0)  and 
that  of  ys  and  y4  be  at  (60  ,  ax)  where  a0  ,  at  ^  0.  From  (25),  as  in  (38),  we 
have  for  any  subdomain  Ra  of  R,  bounded  by  Ga 


(A18) 


0  =  JJ  CWJKWw  +  WJddda 


J 


-KW*edO  +  2We  W„da  -  W2a  dd) 
g»  2 


+  Jf{iW*(KC)„  -  We  Wa  KCd  -  W!  CJ  ddda. 


We  take  Ra  as  the  domain  bounded  by  y1 ,  y2  ,  y3  ,  yt  and  the  segment  Sd 
cut  out  by  a  =  a0  —  <3.  Then  substituting  W„  —  co*,  We  =  u>,  yields 


r  CdO  , -  C  CdO  , - 

0  =  (F - (V —  Koj  to*)2  ct - —  ( V  —  Kw  —  co*)2 

J  Vl+y,  2  J  y,+yt  2 

(A19)  +  f  -  (Kco2dd  +  2Kcoco*da  -  w*2dO) 

J  sa  2 


+  jj{W(KC)a  -  co(o*KCe  -  \<o*2C0)dQdo 

where  yx  and  y3  have  positive  slope  V—  K  and  yt  and  yt  have  negative 
slope  —  V—  K,  see  (30). 

We  take  C  =  1/K.  The  integrals  along  Sa  are  arbitrarily  small  as  6  -*  0. 
For,  da  —  0,  co  and  co*  are  easily  shown  to  be  bounded  by  (32),  and  thus  the 
integral  is  bounded, 


const. 

<^) 


f  de. 

J  si 


But  the  length  of  5,  is 

jsdO  =  0(2  jdV-K(o)  ■  da)  =  0(A) 

since  dK\da  ^  0.  Hence  the  integral  along  S{  is  O(A)  and  has  limit  zero. 

On  yx  and  yt  the  line  integrals  vanish  since  W  =  0.  On  y,  and  y4  the  line 
integrals  are  non-negative  since  C  <  0,  dd  >  0.  The  area  integral  reduces  to 

Ha 

But  it  is  positive  unless  co*  =  0  in  Rt  or  in  the  limit  co*  ==  0,  co  ==  0  in  R. 
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Appendix  3.  Flow  with  Weak  Shocks 


We  shall  show  that  we  cannot  lift  the  continuity  restrictions  (14) 
on  cox  or  a>„  in  the  supersonic  region  of  a  continuous  flow  and  obtain  a  per¬ 
turbation  solution  representing  a  flow  with  a  weak  shock. 

Let  us  suppose  that  the  perturbation  velocities  cox ,  <av  have  some 
discontinuities  in  the  supersonic  region,  that  is  for  a  <  0.  For  a  weak  shock 
we  may  neglect  changes  of  entropy  and  vorticity,  see  [5],  §117,  and  derive 
the  shock  conditions  by  integrating  the  equations  of  motion.  This  yields 


(A22) 


dy  dx 

Agu  •  - - Aqv  •  —  =  0, 

ds  ds 


dx  .  dy 

Au  • - 1 -  Av  •  —  =  0 

ds  ds 


where  A  indicates  the  jump  across  the  shock  and  the  derivatives  are  along 
the  shock  path.  Equations  (A22)  may  be  written  as 


(A23) 


A 


dW  _  A  d0 
ds  ds 


Since  the  unperturbed  flow  is  shockless  and  thus  <p  and  y>  have  con- 

30 1 

tmuous  derivatives  it  follows  that  &»  —  — 

dr 


dW 

and  co*  =  — 

t-o  dr 


r-0 


satisfy 


(A24) 


.  dco  ,  dm* 

A  —  —  A - =  0. 

ds  as 


In  terms  of  0  and  a,  these  equations  may  be  rewritten  as 
(A25) 

and  using  (23)  we  obtain 
(A26) 


dd  do 

0  =  A(oe—  +  Acoa  — , 
ds  ds 

_  .  *  dd  *  do 

0  —  Aa>9  — — [-  A(Oa  — 
ds  ds 


dd  do 

0  =  Aco9—+  Acog—, 
ds  ds 

.  .  dd  do 

0  =  Aa)g- - KAcOf-—. 

ds  ds 

Since  Acoe  and  Aco,  are  not  both  zero  it  follows  that 
(A27)  dd1  +  Kdo*  =  0, 

that  is,  the  shock  path  is  characteristic  and  furthermore  the  derivative  of 
w  along  the  characteristic  is  continuous. 

Suppose  there  exist  several  weak  shock  in  the  flows  &  and  that  has 


08 


CATHLEEN  S.  MORAWETZ 


no  shocks,  see  §1.  Then  co  will  be  a  solution  of  Problem  III  except  that 
either  (14)  holds  or  (A26)  and  (A27). 

Let  Qa  be  the  domain  in  the  Q,  cr-plane  described  in  §3  but  slit  along  the 
shock  paths.  Consider  the  integral 

JJ (Ba>t  +  C(og){K(ow  +  o)m)  d6da 

Os 

and  integrate  by  parts.  Equations  (38)  and  (39)  will  be  changed  only  by  the 
addition  of  the  line  integrals  I,  along  both  sides  of  the  shock.  By  the 
characteristic  condition  (A27)  this  integral  may  be  written  as 


/.= 


(-  B  da-  C  d6 )  (a>„  ±  V-  Ka >a)2 


-*<-.*-<»>£)  */(*)*. 

where  the  integral  is  taken  around  the  slit.  Since  by  (A24),  dco/ds  is  conti¬ 
nuous  across  the  shock  it  follows  that  Is  =  0.  Therefore  equations  (38)  and 
(39)  hold  and  again  w0  =s  a>a  =  0  or  to*  =  0  in  Q.  The  same  argument  also 
yields  the  same  result  for  the  domain  Q*.  This  means  that  the  perturbation 
problem  is  still  incorrectly  posed  even  if  weak  shocks  are  admitted. 
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On  the  Motion  of  Small  Particles  in  a  Potential 

Field  of  Flow 

ABRAHAM  ROBINSON 

University  of  Toronto 


1.  Introduction 

The  analysis  of  the  motion  of  small  particles  which  are  suspended  in  a 
moving  fluid  is  of  importance  for  the  problems  of  aircraft  icing,  and  it  may  also 
assist  in  the  understanding  of  the  silting  up  of  rivers  and  estuaries.  There 
exists  a  considerable  body  of  work  on  this  subject,  both  published  and 
unpublished  (see  the  list  of  references  at  the  end  of  the  paper).  There  are 
individual  differences  as  to  the  equations  by  which  the  motion  of  the  par¬ 
ticles  is  supposed  to  be  governed,  but  all  the  papers  on  this  problem  which 
are  known  to  the  writer  involve  the  computation,  for  the  most  part  step-by- 
step,  of  the  trajectories  of  the  individual  particles. 

In  the  present  paper,  we  assume  that  the  equation  of  motion  of  a  particle 
is  given  by  Stokes’  law,  or  more  generally,  that  the  fluid  exerts  on  the  par¬ 
ticle  a  drag  force  which  is  proportional  to  the  vector  difference  between  the 
velocity  of  the  particle  and  the  fluid  velocity  (compare  [1,2]).  However, 
instead  of  confining  ourselves  to  the  trajectory  of  an  individual  particle, 
we  consider  the  (virtual)  field  of  flow  which  would  be  produced  by  a  conti¬ 
nuous  distribution  of  small  particles  of  given  size  each  governed  by  the  speci¬ 
fied  equation  of  motion.  We  derive  various  general  properties  of  this  virtual 
field  of  flow,  including  a  counterpart  of  Kelvin’s  theorem  on  circulation. 
Moreover,  we  obtain  a  compact  formula  for  the  total  mass  of  the  particles 
deposited  on  a  slender  obstacle  or  on  a  low  elevation  in  a  plane  wall. 

2.  General  Analysis 

Let 

<t>0  =  y>  *>  t) 

be  the  velocity  potential  of  the  field  of  flow  F0  of  a  compressible  or  incom¬ 
pressible  fluid,  where  x,  y,  z,  denote  the  space  coordinates  and  t  denotes 
the  time.  Let 

q0(*.  y.  z>  0  =  K  »  vo  .  wo)  =  “  grad  to 

be  the  corresponding  velocity  vector. 
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Let  m  be  the  mass  and  r  =  (x,  y,  z)  the  position  vector  of  a  small 
particle  which  moves  in  the  given  field  of  flow,  F0.  We  assume  that,  owing 
to  its  smallness,  the  particle  does  not  interfere  with  the  field  of  flow  of  the 
fluid,  and  that  its  equation  of  motion  is  given  by 

da  r 

(2.1)  m—  =  -k(q-q0) 

where 

.  .  dr 

q  =  [u,  v,  w)  =  — 


is  the  velocity  of  the  particle  and  k  is  a  constant.  In  particular,  if  the  particle 
has  the  shape  of  a  sphere  of  radius  R,  then  according  to  Stokes’  law 

(2.2)  k  =  6  nfiR 


where  /i  is  the  coefficient  of  viscosity  of  the  fluid. 

(2.1)  represents  a  system  of  three  ordinary  differential  equations  of  the 
second  order  for  the  quantities  x,  y,  z  regarded  as  functions  of  t.  Thus,  the 
solution  of  (2.1)  is  made  determinate  by  the  specification  of  the  six  quan¬ 
tities 


(2.3)  x,  y,  z,u  —  — ,  v 

at 

for  a  given  time  t  =  t1. 

Now  suppose  that  the  functions 


(2.4)  =  iq(x,  y,  z),  vx  =  vx{x,  y,  z)  wx  =  w^x,  y,  z) 


are  specified  in  some  region  of  space.  We  may  then  regard  x,  y,  z,  iq  ,v1,w1 
as  a  set  of  initial  values  (2.3)  for  arbitrary  but  fixed  t  =  tx  ,  so  that  (2.1) 
determines  x,  y,  z,  u,  v,  w  for  subsequent  times  t.  Hence  q  =  (u,  v,  w)  may 
be  regarded  as  the  velocity  vector  of  a  (virtual)  fluid,  where  in  particular 

q  =  qi  =  K . vi , ^i)  for  t  =  h- 

Also,  in  this  sense,  we  may  replace  (2.1)  by 


(2.5) 


Dq 

=  -  *(q  - 


where  DqjDt  denotes  differentiation  following  the  motion  of  the  virtual 
fluid.  In  scalar  notation 

.  ,  Du  idu  du  du  du\ 

(2.6)  =  m  +  u-  +  v-  +  w  =-*(«-  »„).  etc. 

We  shall  denote  the  field  of  flow  of  the  virtual  fluid  by  F. 
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Now  let  C  be  a  closed  curve  in  the  field  F,  and  let  r  and  7^  be  the  cir¬ 
culation  round  C  of  the  velocity  vectors  q  and  q0 , 

r  =  J  qis  =  J  udx  +  vdy  +  wdz, 

F0  =  j  q0is  =  u0dx  +  v0dy  +  w^Lz. 

J  c  J  c 

We  consider  the  rate  of  change  of  r  in  time  as  C  moves  with  the  field 
F  (i.e.,  with  the  virtual  fluid).  Then 


Now,  by  (2.5) 

/.if  *  "  -  i  J>  “  q°)da  “  “  l  (r~  r»> 

while 

icqi^‘  =  icqiq  =  iSc“m  =  ,>- 

Hence 

dr  k 

<27) 

Moreover,  we  shall  assume  that  ro  is  a  constant,  and  hence  we  can  replace 

(2.7)  by 

5  (r“ 

(2.8)  r=r0  +  ye-*‘/m 

where  y  is  a  constant  which  may  depend  on  C.  By  the  theorems  of  Stokes 
and  Kelvin,  (2.8)  holds  provided  only  the  field  F0  is  irrotational  in  the  region 
swept  by  C  during  the  time  under  consideration.  More  particularly,  if  the 
field  F0  possesses  a  one-valued  velocity  potential  <j>0  ,  then  ro  =  0  and  there¬ 
fore 

(2.9)  r  =  ye~ulm. 

It  follows  that  if  the  circulation  r  round  the  given  curve  C  equals  zero  at 
any  time  then  71  =  0  always. 

If  we  are  given  initially  only  the  conditions  (2.3)  for  a  single  particle 
at  time  t  =  tl,  then  we  can  always  assume  F  =  0  at  time  t  =  t1 ,  and  hence 
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at  all  subsequent  times  by  defining  that  the  velocities  be  equal  to  the  con¬ 
stants  u,  v,  w  specified  by  (2.3)  everywhere  in  the  field  F  at  time  t  —  tr . 
More  generally,  we  have  shown  that  if  the  virtual  field  of  flow  F  is  irrota- 
tional  at  any  one  time,  then  it  is  always  irrotational,  and  hence,  possesses  a 
velocity  potential  <f>.  It  will  be  assumed  from  now  on,  that  this  condition  is 
indeed  satisfied.  Thus 

(2.10)  curl  q  =  0,  q  =  —  grad  <f>. 

It  follows  that  we  may  replace  (2.6)  by 


with  two  similar  equations.  Integrating,  we  obtain  the  following  counter¬ 
part  of  Bernoulli’s  equation: 


(2.12)  m 


d<t>  llUV  /d<j>\2  (dd>\2\) 

-4  +  iUJ  +(4)+(a0)f  -*»-«+*«• 


Conversely,  it  is  not  difficult  to  verify  that  if  <f>  is  a  solution  of  (2.12), 
for  specified  <f>0  ,  then  q  =  —  grad  tf>  is  a  solution  of  (2.6). 

For  steady  conditions,  we  obtain 


(2J3)  ((S)  +  (§)  +  (S)  )  ~  k[(f>  ~(f>o)+H 


where  H  may  now  be  taken  to  be  independent  of  the  time. 


3.  Flow  Past  an  Obstacle 

Equation  (2.13)  may  also  be  written  as 
(3.1)  \m(u2  +  v2  +  w2)  =  k(<f>  —  <j}0)  +  H. 

Put 

d<f>' 

<f>' =  <f>  —  <j>0 ,  u'  =  u  —  u0=  —  —,  etc., 
then  (3.1)  yields 


(3.2) 

where 


W(K  +  +  uou'  +  vov'  +  u>ow')  =  H'  +  H 


q0  =  Vu%  +  v2  +  q'  =  V  u'2  +  v'2  +  w'2. 


We  shall  now  suppose  that  q'  is  small  compared  with  q0  .  This  amounts 
to  the  assumption  that  the  difference  between  the  velocity  of  the  particle 
and  the  velocity  of  the  fluid  is  of  the  nature  of  a  (first  order)  correction  to  the 


MOTION  OF  PARTICLES  IN  A  FIELD  OF  FLOW 


73 


fluid  velocity.  For  example,  this  assumption  will  be  satisfied,  generally 
speaking,  in  the  case  of  uniform  fluid  flow  past  a  slender  obstacle,  provided 
the  particle  is  at  rest  relative  to  the  fluid  initially.  The  assumption  breaks 
down  in  the  close  neighborhood  of  a  stagnation  point  but,  as  in  thin  aerofoil 
theory,  the  overall  results  may  still  be  expected  to  be  valid. 

Omitting  q '2  in  (3.2)  we  obtain  for  the  steady  case, 

(3.3)  m(u0u'  +  v0v'  +  w0w’)  =  k<f>'  +  H  —  . 

Now  the  equations  of  a  streamline  in  the  field  F0  are 


dx  u, 
ds  q0 


o  W 


dz  Wn 


ds  q0  ds  f/o 


where  s  measures  the  arc  of  the  streamline.  Also,  differentiating  <f>{)  with 
respect  to  s,  we  obtain 

d<f>0 

Hence,  along  the  given  streamline 

dy 


wn 


dx  u0  dy  v0  dz 
Wo~  Wo~  fo  Wo  ~ 

Accordingly,  we  may  modify  the  left  hand  side  of  (3.3)  in  the  following 


«8‘ 


way: 


'  m(u0u'  +  v0v'  +  w0w')  =  nuf0  ^ 


dx  d<f>'  dy  d<j>'  dz  d<f>' 


\d<f>0  dx  d<f>0  dy 
Inserting  this  in  (3.3),  we  obtain 


dz_df\_ 

d<f>0  dz)  m  0  d<f>0 


(3.4) 


nufojr  =  ki'+H- 
“<Po 


which  is  an  ordinary  linear  differential  equation  of  the  first  order  for  </>' 
as  a  function  of  (f>  along  any  given  streamline.  The  general  integral  of  (3.4)  is 

w  IJSI/HS-IM-JSIK 

Now  for  the  case  of  uniform  fluid  flow  past  a  slender  obstacle,  with  main 
stream  velocity  U  in  the  direction  of  the  a;-axis,  we  may  write 


<f>o  =  —  Ux  +  <l>0  <  «o  =  U  — 


K 

dx  ’ 


Vn=  — 


dy' 


wR  = 


dz 
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The  induced  velocities  —  ^2,  —  ^2,  —  ^2  are  small  compared 

ox  dy  dz 

with  U  (except  near  a  stagnation  point)  and  therefore 

where  we  have  omitted  terms  of  the  second  order  of  smallness  compared 
with  UK  Also,  if  the  particle  is  at  rest  relative  to  the  fluid  far  upstream 
of  the  obstacle,  then 

u  =  U,  v  =  w  =  0,  </>'  =  <f>  —  <f>0  —  0 

in  that  region,  and  so,  applying  (3.1), 

ImU*  =  H. 


Hence 


-*  =  * 


U 2 


H 


*»*•  *  ‘  *1  U‘-2U< 

'  dx 


■  u  dx' 


Accordingly,  (3.5)  becomes  approximately 


d<}>0- 


The  constant  in  the  integral  J  ( k/mU 2)  d<f>0  is  at  our  disposal,  and  so  we 
may  choose 

f*#o 


f  kd<f>  o 

J  mU‘ * 


—  Ho 


where  A  =  kjmU 2.  Then 

(3.7)  + 


since 


-^2  =  C7  4- 

^  a*  ’ 

To  insert  the  appropriate  limits  of  integration  in  (3.7),  we  bear  in  mind  that 
as  x  tends  to  —  oo,  <f>0  oo  and  <f>'  ->  0.  Hence 


(3.8) 


where  we  have  used  $0  to  distinguish  the  variable  of  integration,  or 


U 

4 
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(3-9)  + 

Once  again,  we  may  simplify  the  analysis  by  taking  into  account  that 
the  streamlines  of  F0  are  nearly  parallel  to  the  ar-axis  and  that,  in  the  first 
approximation, 

-  Ux,  d-p  =  -U. 

OX 

Hence 

(3.10)  *-*•-[_  {U+d^)e~UHX~X)dX’ 

where  we  use  X  to  indicate  the  variable  of  integration.  However,  this  addi¬ 
tional  simplification  will  not  be  required  for  the  analysis  of  the  next  section. 


4.  Incompressible  Flow  in  Two  Dimensions 


Let  us  now  assume  that  the  fluid  flow  is  two-dimensional  and  incom¬ 
pressible,  and  takes  place  in  planes  parallel  to  the  x ,  y-plane.  We  may  then 
use  z  to  denote  the  complex  variable  x  -j-  iy.  Also,  <f>0  now  is  the  real  part  of 
a  complex  potential  w0 , 

wo  —  fa  +  *Vo 

and  we  may  write 


Wo  =  0 o  +  iW0 


where  &0  is  the  variable  of  integration  in  (3.9).  But  the  integration  in  (3.9) 
is  performed  along  a  streamline,  and  so  0O  =  const  =  rp0,  throughout. 
Also 


where  dwjdz  is  an  analytic  function  of  z,  and  hence  of  w0  .  Thus 
j"  fu  +  ^  d0o  =  die  (u  +  dWo 

and  therefore 

(4.1)  j,  =  ®e  |w0  -  i  J”  dW 0}  . 

If  the  upper  limit  of  the  integral  in  (4.1)  were  a  fixed  finite  complex 
number,  it  would  follow  immediately  that  the  expression  within  the  curly 
brackets  is  an  analytic  function  of  w0,  and  some  reflection  shows  that  this 
is  true  even  in  the  present  case.  Hence  <f>  is  the  real  part  of  a  complex  po¬ 
tential  w,  where 
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(4.2)  w  =  w0  —  j  e  A!W,°  W||)  dWo 

or 

l  /•*o  /  \  ,fry 

(4.3)  «*(*„)  =  t*0(z0)  +  i  J  (U  +  p  *. 

The  lower  limit  of  the  integral  in  (4.3),  —  oo,  indicates  that  the  integral 
is  to  be  taken  from  a  point  far  upstream  of  the  obstacle.  Apart  from  this, 
the  complex  representation  has  the  advantage  that  the  path  of  integration 
may  now  be  deformed  at  will  in  any  finite  region  of  the  field  F0 . 


5.  The  Mass  Flow  of  Particles 


Assume  now  that  particles  of  equal  size  (i.e.  of  equal  m  and  k)  are 
distributed  uniformly  and  densely  in  the  region  far  upstream  of  the  obstacle. 
We  may  then  define  the  density  q  of  our  virtual  fluid  as  the  average  total 
mass  of  particles  per  unit  volume.  The  variation  of  q  following  the  motion 
of  the  fluid  is  given,  as  usual,  by  the  equation  of  continuity 

1  Dp 

<51>  -eu~~ divq- 

But  under  the  conditions  of  the  preceding  section,  <f>  is  the  real  part  of 
a  complex  analytic  function  and,  therefore,  satisfies  Laplace’s  equation 

A(p  =  ~  div  q  =  0. 

Hence 


(5.2) 


Dq 

Dt 


—  0, 


showing  that  the  virtual  fluid  is  incompressible. 

Equation  (5.2)  is  only  approximately  true.  It  can  be  proved  more 
rigorously  that  if  the  real  fluid  is  irrotational  and  incompressible,  and  if  the 
virtual  fluid  is  at  rest  initially  (at  time  t  —  0,  say)  then 

(5.3)  div  q  0. 

This  shows  that  <f>  is  a  subharmonic  function, 

(5.4)  A<f>  =  —  div  q  ^  0. 

In  order  to  establish  (5.3),  we  apply  the  operator  div  to  (2.5).  This 
yields,  on  the  right  hand  side 

(5.5)  —  k  div  (q  —  q0)  =  —  k  div  q, 
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while  on  the  left  hand  side 


(5.6) 

where 


m 


/Dq\  ID  \ 

div  —  =  m  I  —  div  q  +  A ) 
\Dt!  \Dt  M  / 


/du\2  dv  du  dw  du  du  dv  /dv\2  dw  dv  du 

IcW  dx  dy  dx  dz  +  dy  dx~^~  \3y/  +  d y  dz^  dz 

/ du\ 2  /dv\2  /dw\2  / du\ 2  ldv\2  Idw\2  t 

\dx)  \dx)  \3a:/  ^  \dy)  \di //  1 3 yj  \ 


dw  dv  3 w 
dx  ^"3^3^^” 


du\2 

rJ  + 


/dw\2 

\dzj 

/ 3v\ 2  /  3te^\ 2 
\3  z)  3z/ 


Thus,  A  can  be  represented  by  a  sum  of  squares  and  therefore 

A  S:  0. 


In  deriving  this  relation,  we  have  made  use  of  the  fact  that  curl  q  =  0, 
by  (2.9).  Combining  (5.5)  and  (5.6),  we  obtain 

D  k 

(5.7)  —  div  q  -) —  div  q  =  —  mA. 

Dt  m 

This  is  an  ordinary  differential  equation  for  div  q,  following  the  motion 
of  the  virtual  fluid,  and  since  the  virtual  fluid  is  supposed  to  be  at  rest  for 
t  —  0,  we  have  the  initial  condition 

div  q  —  0  for  t  =  0. 

With  this  condition,  the  solution  of  (5.7)  is 

(5.8)  div  q  =  —  me~ktlm  [* A'eht,mdt. 

J  o  ■ 

Thus, 

div  q  0 


as  asserted  by  (5.3). 
(5.9) 


It  then  follows  from  (5.1)  that 


Dp 

—  >  0 
Dt 


which  shows  that  the  density  of  the  virtual  fluid  is  a  non-decreasing  function 
of  the  time,  following  the  motion  of  the  fluid.  However,  in  the  sequel  we 
shall  adopt  the  approximate  conclusion  (5.2)  according  to  which  q  is  actually 
constant.  Then  the  mass  flow  of  particles  across  a  curve  connecting  two 
points  A  and  B,  in  unit  time,  is  given  by  the  product  of  the  density  q  and 
of  the  difference  yA  —  fB,  where  y>  —  Sm(w)  is  the  stream  function  of 
the  field  F, 


(5.10) 


M  —  —  q  Sm(wB  —  wA). 
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However,  this  formula  becomes  inapplicable  if  any  of  the  streamlines 
of  F  which  cross  the  curve  connecting  A  and  B  have  already  crossed  the 
solid  boundary  of  the  obstacle  further  upstream.  Excluding  this  case,  we 
assume  more  particularly,  that  A  and  B  are  situated  on  the  same  streamline 
of  the  field  F0  .  Then 

Sm(w0(zB)  —  w0(zA))  =  0 


where  zA,  zB  are  the  complex  numbers  corresponding  to  A,  B.  Inserting 
(4.3)  in  (5.10),  we  have 


(5.11) 

1 

II 

* 

[  e^B)  P  P  +  d™°\  dW°  e-m*dz 
l  J  _»  V  dz  1  dz 

and  modifying  (5.11)  slightly 

=u|b 

1 

II 

|  e-wodz 

(5.12) 

We  now  consider  flow  without  circulation  round  a  symmetric  obstacle 
at  zero  incidence;  we  take  A  as  the  front  stagnation  point  and  B  as  a  point 
on  the  top  surface  of  the  obstacle.  We  may  choose  A  as  origin  of  coordinates, 
then  the  negative  x-axis  is  a  streamline;  we  may  assume  Jm{w)  =  y>0  =  0 
along  that  line.  It  follows  that  y>0  =  0  also  on  the  top  surface  of  the  obstacle, 
more  particularly  between  A  and  B.  Also  dWJdz  is  real  on  the  negative 
ar-axis  and  so  the  first  term  within  the  curly  brackets  of  (5.12)  is  real. 
Accordingly,  the  expression  for  M  reduces  to 


M  —  - 


dW„\dW,  _ 


(5.13) 

v  l  J  0 

since 

Wo  e-ZW, 

U  dz  j  dz 


e~AW°  dz 


dz 


is  real  along  the  top  surface  of  the  obstacle.  But 
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and,  since  the  obstacle  is  slender,  we  have  the  linearized  boundary  condition 

®o  =  Uf'(x) 

where  y  —  f(x)  is  the  equation  of  the  top  surface.  Hence  (5.13)  may  be 
replaced  by 

r*B  dW 

(5.14)  M  —  —  q  \  f'(x)e~^Wo~w o(‘b)) - -dz. 

'  Jo  dz 

Finally,  replacing  W0  —  w0 (zB)  and  (dWJdz)dz  by  the  first  approximations 
—  U  (x  —  xB )  and  —  Udx,  and  bearing  in  mind  that  X  =  k/mU2  we  obtain 

(5.15)  M  =  qU  f'(x)e-k(*B-*)l™u  dx 

where  xB  =  Ste  (zB).  This  is  a  compact  formula  for  the  total  mass  of  the 
particles  which  strike  the  top  surface  of  the  obstacle  between  A  and  B,  in 
unit  time.  This  formula  applies  equally  to  an  unstaggered  cascade  of  geo¬ 
metrically  similar  symmetric  obstacles  at  zero  incidence,  and  to  an  elevation 
in  a  straight  wall  (e.g.  river  bottom)  which  is  parallel  to  the  main  stream. 

Using  (3.10)  we  might  have  arrived  at  the  final  formula  (5.15)  also  by 
a  less  careful  analysis. 


6.  Examples 


Write 

£  =  X/Xg  ,  T  =  f(XB)/Xg  ,  g(£)  =  f(XB^)lf(XB)- 

Then 

/»  =  /'(f)/(*a)/*s  =  *?'(£)• 

Substituting  in  (5.15),  we  obtain 

M  =  qUxxb  fVd)*-0'1^  dt 
(6-1)  J°,a 

=  QUt(xB)jog'(£)e-°Md! 

where  a  is  the  non-dimensional  parameter  kxBlmU.  As  [cr  tends  to  zero, 
M  tends  to 

(6.2)  M0  =  eUf(xB)  fo  g’(£)  d£  =  qUf[xB). 

This  is  in  agreement  with  the  fact  that  for  vanishingly  small  viscosity 
the  particles  retain  the  main  stream  velocity,  so  that  M  becomes  equal  to 
the  product  of  density,  main  stream  velocity,  and  frontal  area.  Accordingly, 
it  is  appropriate  to  restate  the  general  result  in  terms  of  the  coefficient 
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Suppose  first  that  the  front  portion  of  the  obstacle  has  the  shape  of  a  wedge, 

f(x)  =  cx,  c  >  0. 

The  shape  of  the  rear  portion  of  the  obstacle  is  irrelevant,  but  we  may 
assume,  for  example,  that  we  are  dealing  with  a  double-wedge  aerofoil.  Then 


and  hence 
(6.4) 


f(xB)  =  CXB  ,  g(f)  =  cx£B/cxB  =  f 


rj  =  f  g-oh-f)  d£  =  e-°  eofj  =  —  (1  —  e~a). 
J  0  l<7  Jo  a 


Next,  consider  the  case  of  an  obstacle  of  elliptic  shape  which  is  given 
by  the  equation 


(6.6) 

Then 


y 


(«  ~  a)2  ,  t  =  i 

a2  +  b2 


and  the  front  stagnation  point  is  located  at  the  origin,  as  required.  We 
choose  xB  at  the  midchord  of  the  obstacle,  xB  =  a,  so  that  f{xB)  =  b.  Then 


g(|)=&j/l- 

g'(f)  = 


(a|  —  a)2 


Vl  -  (1  -  f)». 


1  -  f 


(6.6) 


Vi  -  (i-l)2 

V=  f  7  C 

J  oVl-(l-f)2  Jo 


VI  -  Cs 


r«. 


It  is  not  difficult  to  verify  that  along  the  entire  segment  AB  of  the  con¬ 
tour  of  the  obstacle,  the  streamlines  of  F  cross  into  the  obstacle  from  up¬ 
stream  so  that  we  are  dealing  with  a  real  case.  The  integral  in  (6.6)  can  be 
evaluated  in  terms  of  standard  transcendental  functions  (see  [6]  p.  136, 
No.  28) 


n 


V  =  1  —  J  (A(«)  -  Li(*))- 

In  this  formula,  I^a)  is  a  modified  Bessel  function  of  the  first  kind, 


(6.7) 
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+  i)W+V 

and  L^ct)  is  a  modified  Struve  function, 

T  /  \  £  (a/2)2n+2 

Ll(<r)“i0r{n+f)r(w  +  f)- 

Table  1  shows  values  of  rj  for  various  values  of  a,  for  wedge  and  ellipse. 
For  a  spherical  particle,  k  is  given  by  (2.2),  according  to  Stokes’  law,  while 
tn  =  (4:j3)nQ'R3,  where  q'  is  the  density  of  the  particle.  Hence 

k  xB  $,UXB 

°  —  ~rn~TJ  ~  2g'R2U  ' 

If  the  particle  is  a  water  droplet,  and  the  surrounding  medium  is  air, 
then  we  may  take  q'  =  1,  /j,  =  1.7  x  10-*  in  c.g.s.  units.  Typical  values 
for  xB  ,  R,  U  are,  in  the  same  units, 

xB  =  3  X  102,  R  =  3  x  10~3,  U  —  1.5  X  10* 

for  the  aircraft  icing  problem.  This  yields  a  =  1.7. 


TABLE 

1 

r)  (wedge) 

7]  (ellipse) 

0.0 

1.000 

1.000 

0.2 

0.906 

0.865 

0.4 

0.824 

0.734 

0.8 

0.752 

0.630 

0.8 

0.688 

0.543 

1.0 

0.632 

0.468 

2.0 

0.432 

0.310 

Although  no  direct  comparison  is  possible,  our  results  appear  to  be  con¬ 
sistent  with  [2],  The  fact  that  our  approximations  break  down  in  the 
neighborhood  of  the  front  stagnation  point  does  not  invalidate  these  results. 
To  see  this,  we  shift  the  point  A  a  short  way  upstream,  and  remove  the  path 
of  integration  AB  from  the  neighborhood  of  the  front  stagnation  point. 
Then  our  approximations  may  be  assumed  to  hold  uniformly  on,  and  up¬ 
stream  of,  the  modifed  AB.  However,  the  actual  flow  across  AB  is  the 
same  as  before,  by  continuity,  and  the  mathematical  result,  given  e.g.  by 
(5.13),  has  not  changed,  by  Cauchy’s  theorem.  This  argument  shows  that 
our  overall  results  may  be  expected  to  hold  in  spite  of  the  breakdown  of  our 
assumptions  near  the  front  stagnation  point. 
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7.  The  Influence  of  Gravity 


The  influence  of  gravity  has  been  excluded  from  our  considerations  so 
far.  Suppose  now  that  a  conservative  body  force,  which  is  given  by  the 
potential  £3,  acts  on  the  particle.  Then  (2.1)  and  (2.5)  have  to  be  replaced  by 


(7.1) 
and 

(7.2) 


—  k(q  —  q0)  —  m  grad  £3 


Dq 

w—  -  -  k(q-  q0) 


—  m  grad  £3, 


respectively.  Now 

Jc  grad  £3  d s  =  0 

when  the  integral  is  taken  round  a  closed  curve  C,  for  one-valued  £3.  Hence 
(2.8)  and  its  general  consequences  still  hold.  In  place  of  (2.12)  and  (2.13) 
we  obtain 


=  k(<f>  —  <j>0)  —  m£3  +  H(t) 


and 


{7,4)  (s) +  (ft)  +  (i) )  =  ^  ~  ^  ~  mQ  +  H- 


\dy, 

Moreover,  if  Q  is  a  harmonic  potential, 


A£2  =  div  grad  £3  =  0, 

then,  by  applying  the  operation  div  to  (7.2),  we  may  still  deduce  (5.3)  under 
the  conditions  stated  in  Section  5. 

Let  us  now  consider  the  particular  case  of  flow  past  an  obstacle  under 
two-dimensional  conditions,  as  in  Section  4.  We  suppose  that  the  conser¬ 
vative  body  force  is  the  force  of  gravity  acting  in  a  direction  perpendicular 
to  the  free  stream.  Then 


(7.5)  Q  =  gy 

provided  we  can  neglect  the  buoyancy  of  the  particle.  If  we  wish  to  take  the 
buoyancy  into  account,  we  only  have  to  replace  g  in  (7.5)  by 


r- 


where  m  is  the  mass  of  the  particle,  as  before,  and  m'  is  the  mass  of  the  fluid 
displaced  by  it.  We  then  obtain  from  (7.4) 
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(,-6)  im  KS) +(§)*] = ~  m  -  ms'y + h 

for  steady  conditions.  If  no  obstacle  is  present,  <f>0  =  —  Ux,  and  (7.6)  is 
satisfied  by 


(7.7) 


m 


6  =  —  Ux  +  —  g'y 


which  corresponds  to  the  steady  fall  of  the  particles  with ,, terminal  velocity” 
(7.8) 

The  corresponding  value  of  H  is 


m  ,  m  —  m 

V=~JS= - ~e- 


(7.9)  H  -  im(U‘  +  Pa) 

and  this  must  also  be  the  value  of  H  for  steady  flow  past  an  obstacle,  as  can 
be  seen  by  considering  conditions  far  upstream.  For  that  case,  we  write 


m 


<£  =  —  Ux  +  —  g'y  +  <f,’  =  —  Ux  +  Vy  +  <f>' 


so  that 


d<b' 

u  =  U  —  ~- 
dx 


Then  (7.6)  becomes 

T-[(sr+d) 


v  =  —  V  — 


36'.  361 

2  U-^~  +  2Vf- 

3x  3yl 


3£ 


k(6’  -  6o) 


where  6o  is  the  induced  velocity  potential  of  the  fluid  flow,  6o  =  <^o  +  Ux. 

We  now  assume  that  the  terminal  velocity  V  is  small  compared  with  the 
free  stream  velocity  U  and  we  may  suppose  that  the  same  applies  to  the 

36'  36' 

induced  velocities  in  the  field  F,  —  —  and  —  — . 

3x  3y 

Then  (7.10)  becomes 

(7.11)  mU^  +  k6'  =  k6i. 


The  solution  of  this  differential  equation  is,  with  the  appropriate  limits  of 
integration, 

(7.12)  6'  =  m  e~U**  jX  „  o  eUXX  dX- 

In  this  formula,  X  =  k/mU *,  as  before,  and  we  have  used  in  order  to 
indicate  that  the  function  depends  on  the  variable  of  integration  X.  The 
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lower  limit  of  integration  has  been  chosen  so  as  to  ensure  that  <f>'  vanishes  far 
upstream,  as  required. 

Integrating  by  parts  in  (7.12), 

f = j [0; ««“] '  -  j‘‘  ‘iaix 

But 

<f>=-Uz+Vy+'{,\  +  =  U  +-—■ 

ox  ox 

and  thus 

(7.13)  <f>  =  ^  -  J*  (u  +  ^f)  e~UM*~X)  dX  +  Vy. 

Comparing  (7.13)  with  (3.1),  we  see  that,  within  the  stated  approxi¬ 
mations,  the  velocity  field  due  to  gravity  is  additive  to  the  field  calculated 
previously.  The  result  still  represents  the  field  of  flow  of  an  incompressible 
fluid.  Thus,  in  order  to  obtain  the  total  mass  of  the  particles  which  strike  a 
segment  AB  of  the  contour  of  the  obstacle,  we  only  have  to  add  to  (5.15) 
the  mass  which  would  strike  the  obstacle  in  unit  time  if  it  were  at  rest.  For 
the  examples  calculated  in  Section  6,  this  quantity  is  given  by  qVxb  . 
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On  a  Generalization  of  the  Normal  Basis  in 
Abelian  Algebraic  Number  Fields* 

M.  NEWMAN  AND  OLGA  TAUSSKY 

National  Bureau  of  Standards 


It  is  known  [1,  2]  that  in  certain  normal  algebraic  number  fields  a 
so-called  ,, normal”  basis  can  be  chosen  for  the  integers  which  consists  of 
the  conjugates  of  one  number.  Let  F  be  the  field  in  question,  and  assume 
that  it  is  of  degree  n  and  that  a(S<>,  i  =  1,  2,  •••,«,  is  a  normal  basis  for  F, 
where  the  S{  are  the  elements  of  the  Galois  group  of  F. 

With  any  integral  basis  oq  ,  a*  ,  •  •  • ,  a,  one  considers  the  matrix 
D  =  (a)s*’ ) ,  i,  k  —  1,  2,  •  •  • ,  n.  Since  the  square  of  the  determinant  of  D 
is  the  discriminant  of  F  we  call  D  a  discriminant  matrix.  Consider  this 
matrix  for  a  normal  basis.  We  then  have 

D  =  ( * <s<s*>),  i,  k  =  1,  2,  •  •  • ,  n. 

This  special  D  is  symmetric  if  and  only  if  the  field  is  abelian.  For  complex 
matrices  symmetry  is  here  used  in  the  ordinary  sense,  not  in  the  hermitian. 

It  will  be  shown  in  §1  that  for  an  abelian  field  this  special  D  is  always 
normal. 1 

The  main  question  (which  is  considered  in  §3)  is  whether  there  is  any 
other  possibility  for  the  discriminant  matrix  to  be  normal.2  It  is  shown 
that  for  cyclic  fields  there  is  no  other  possibility. 

In  §2  the  uniqueness  of  the  normal  basis  is  studied.  It  is  shown  that 
for  n  —  2,  3,  4,  6  the  normal  basis  is  unique  apart  from  permutations,  and 
that  for  all  other  «  the  normal  basis  is  not  unique. 

In  the  general  case  the  problem  studied  here  leads  to  the  study  of  group 
matrices;  for,  the  discriminant  matrix  derived  from  a  normal  basis  is  a 
permutation  of  a  group  matrix.  In  the  case  of  cyclic  fields  we  are  led  to 
circulants. 

§1.  theorem  1.  Let  F  be  an  abelian  field  of  degree  n  with  a  normal 
basis  a(S<)  and  D  the  discriminant  matrix  corresponding  to  this  basis.  Then 
D  is  normal. 


•This  work  was  supported  (in  part)  by  the  Office  of  Naval  Research. 

*A  matrix  A  is  called  normal  if  A  A*  =  A*  A  where  A*  =  A',  the  transposed  conjugate 
of  A.  A  complex  matrix  A  with  A  «=  A'  is  in  general  not  normal. 

•See  [3]. 
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Proof :  As  mentioned  in  the  introduction,  it  is  clear  that  D  is  symmetric 
for  abelian  fields.  It  remains  to  show  that  D  is  normal  even  if  F  is  a  complex 
field.  The  (i,  k)  element  of  DD*  is 

J  a(s,sr)  a<s*sr>. 

f"l 

The  corresponding  element  of  D*D  is  the  conjugate  complex  value.  It  will 
now  be  shown  that  the  element  is  real.  Let  a(SiSr)a(S*Sr)  be  complex  for  some 
value  of  r.  We  will  see  that  the  sum  also  contains  its  conjugate.  Among 
the  S{  there  is  a  substitution  Sm  with  S£  =  1  such  that  Q{Sm)  =  q  for  all  q 
in  F.  We  then  have  by  virtue  of  the  commutativity  of  the  Galois  group  of  F 

a<sisr>  a(s*s'>  =  a<su>sisr)  a(sis*sr) 

=  a«s‘Ws‘V 

where  Sr  =  Sa  Sr .  This  proves  the  theorem. 

§2.  Any  integral  basis  is  obtained  from  the  given  normal  basis  by 
means  of  a  substitution  ( aik )  where  aik  are  rational  integers  and  \aik\  =  ±  1. 
When  searching  for  an  integral  basis  with  normal  discriminant  matrix 
we  may  simply  find  a  permutation  of  the  given  normal  basis  or  another 
normal  basis. 


theorem  2.  Let  F  be  a  cyclic  field  of  degree  n  with  normal  basis  a(Sf). 
If  n  —  2,  3,  4,  6,  the  normal  basis  is  unique.  For  all  other  n  the  normal  basis 
is  not  unique. 

Proof:  Assume  that  =  S'  where  S  is  a  generator  of  the  cyclic  Galois 
group  of  F.  Consider  a  unimodular  matrix  (aik)  of  rational  integers  such  that 
(aik)(tx.{S,))  is  again  a  normal  basis. 

We  prove  first  the  following  lemma: 


lemma  1.  Let  (aik)  be  a  unimodular  matrix  of  rational  integers  and  alS() 
a  normal  basis  in  a  cyclic  field.  In  order  for  (aik)  (a(S<))  to  be  a  normal  basis  also 
it  is  necessary  and  sufficient  that  (aik)  be  the  product  of  a  circulant  matrix  and 
a  permutation  matrix. 

Proof:  A  circulant3  matrix  is  usually  assumed  of  the  form 


(#1  »  »  1  ”  i  ®n) 


ax  a%  •  • '  an 

1 


In  the  cyclic  case  the  discriminant  matrix  of  a  normal  basis  multiplied 
on  the  left  by  a  suitable  permutation  matrix  becomes  a  circulant.  Converse- 


“For  a  detailed  study  of  circulants  see  e.g.  [4]. 
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ly,  if  the  discriminant  matrix  is  a  permutation  of  a  circulant  then  the 
basis  from  which  it  is  derived  is  normal. 

Assume  then  that  the  discriminant  matrix,  multiplied  by  a  suitable 
permutation  matrix  is  the  circulant  (a<s<~*>).  Use  now  the  fact  that  the 
product  of  circulants  and  the  inverse  of  a  circulant  is  again  a  circulant.  It 
follows  that  (aik)  (a(S<))  will  be  again  a  normal  basis  if  and  only  if  (aik)  is  a 
permutation  of  a  circulant. 

From  Lemma  1  it  follows  that  the  non-uniqueness  of  the  normal  basis 
depends  on  the  existence  of  a  non-trivial  unimodular  circulant  of  rational 
integers.  (A  circulant  is  called  trivial  if  it  has  only  one  non-zero  entry  in 
each  row).  Theorem  2  is  therefore  a  consequence  of 

theorem  2'.  An  nxn  circulant  of  rational  integers  cannot  be  unimodular 
for  n  =  2,  3,  4,  6  unless  it  is  trivial.  For  all  other  values  of  n  there  are  non¬ 
trivial  unimodular  circulant  sA 

Proof:  A  circulant  {a^  ,  az  ,  •  •  • ,  an)  is  unimodular  if  and  only  if  ax  -f- 
a2£  +  •  •  •  -f-  aJln~1  is  a  unit  in  the  field  generated  by  the  w-th  roots  of  1, 
for  each  of  the  n  values  of  £  such  that  £"  =  1.  Since  for  n  —  2,  3,  4,  6  the 
field  of  the  «-th  roots  of  unity  does  not  contain  any  units  of  infinite  order  it 
follows  that  for  these  degrees  all  a{  but  one  are  zero. 

Suppose  now  that  n  —  2l‘n1 ,  where  %  is  odd  and  5:  5.  Then  there  is 
an  odd  k  with  ( k ,  nx)  —  1  such  that  3  sS  k  gi  %  .  For  this  k,  it  is  easily 
shown  that  e  =  (1  +  ) /  ( 1  +  C2")  is  a  unit  in  the  field  of  the  w-th  roots 

of  1,  for  each  of  the  n  values  of  £  such  that  £*  =  1.  Thus  the  circulant  based 
on  e  is  unimodular  and  non-trivial.  This  leaves  only  the  n’s  of  the  form 

2'',  3  •  2"  for  consideration.  In  view  of  the  fact  that  if  a1A-  a.t£  -\ - -f  an£n-1 

is  a  unit  in  the  field  of  the  «-th  roots  of  1  for  each  of  the  n  values  of  £  snch 
that  £"  =  1  then  «i+  a2(C1,k)k  -)-•••  +  an(^llk)k{n~1)  is  a  unit  in  the  field  of 
the  nk- th  roots  of  1  for  each  of  the  nk  values  of  £1,lc  such  that  (£i/*)»*  =  l, 
it  is  only  necessary  to  exhibit  a  non-trivial  unimodular  circulant  for  n  =  8 
and  n  =  12.  For  n  —  8,  the  circulant  (1,  1,  0,  —  1,  —  2,  —  1,  0,  1)  is  uni¬ 
modular  and  non-trivial.  It  is  derived  from  the  unit  (1  -(-  V 2)2.  For  n  =  12, 
the  circulant  (3,  2,  1,  0,  —  1,  —  2,  —  2,  —  2,  —  1,  0,  1,  2)  is  unimodular  and 
non-trivial.  It  is  derived  from  the  unit  (2  -f  V 3)2. 

§3.  We  now  deal  with  the  main  question.  Denote  the  symmetric 
discriminant  matrix  of  the  fixed  normal  basis  of  an  abelian  field  by  D  and  put 
(aik)D  =  A.  The  problem  is  to  find  a  unimodular  matrix  (ailc)  of  rational 
integers  such  that  A  A*  —  A*  A.  We  prove 


4We  are  indebted  to  Dr.  S.  Gorn  for  turning  our  attention  to  the  problem  of  the  unique¬ 
ness  of  the  normal  basis.  Recently  H.  P.  F.  Swinnerton-Dyer  found  another  proof  of  Theorem  2'. 
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theorem  3.  Let  F  be  an  abelian  field  of  degree  n  with  a  normal  basis. 
In  order  for  the  discriminant  matrix  A  =  ( aik)D  to  be  normal  it  is  necessary  and 
sufficient  that  (  be  a  group  matrix  corresponding  to  the  Galois  group 

of  F. 

Proof5:  The  equation  A  A*  =  A*  A  leads  to 

(aik)DD*(aikY  =  />*(**)'(**)£>. 

We  next  note  that  DD*  is  a  matrix  consisting  of  rational  numbers.  For, 
let  be  the  normal  basis  from  which  D  is  derived.  Then  DD*  = 

(a<s<s*>)(a(s‘s*))',  i.e.  the  (i,  £)-element  of  DD*  is  2  a(S(S,)a(S‘Sf).  If  we  now 

f— 1 

replace  a  by  a(S)  where  S  is  any  one  of  the  Sr  we  obtain  in  virtue  of  the 
commutativity  of  the  Sr 

'  J  a(s<ss,)  a  (S*SSr) 

r-l 

which  coincides  with  the  previous  expression  since  SSr  runs  through  all 
elements  of  the  Galois  group  with  Sr  .  Thus  DD *  is  rational. 

Hence  it  follows  that  also  the  right  side  D*(aik)’(aik)D  is  unaltered  if 
a  in  D  is  replaced  by  a(S)  where  5  is  any  of  the  elements  of  the  Galois  group  of 
F.  However,  replacing  a  in  D  by  a(S)  effects  a  permutation  of  the  rows  of 
D,  i.e.  can  be  written  in  the  form  PSD,  where  Ps  is  a  permutation  matrix. 
We  thus  obtain 


C) 


Psiai*)’ iaik)  Ps  —  iaik)'(aik) 


where  the  Ps  are  n  x  n  permutation  matrices  which  form  a  representation 
of  the  Galois  group  of  F.  In  the  case  that  F  is  cyclic  the  Ps  can  be  taken  as 
the  powers  of 


0 

1 

0  • 

•  0 

0 

0 

1  • 

•  0 

0 

0 

0  • 

•  1 

1 

0 

0  • 

•  0 

P  = 


It  is  then  known  that  the  previous  relation  implies  that  ( aik)'(aik )  is  a  circu- 
lant.  In  the  general  abelian  case  it  can  be  concluded  that  (aik)'(aik)  is  a 
group  matrix  with  respect  to  the  Galois  group  of  F.  For,  let  Ps  be  the 

(1  2  •  •  •  ft  \ 

).  Then  (*)  implies  that  the  matrix  obtained  from 

V*  '  ‘  ’  sn/ 

(«•)'(«*)  by  permuting  the  rows  according  to  Ps  is  equal  to  the  matrix 
obtained  from  [aik)'(aik)  by  permuting  the  columns  according  to  Pj1 .  It 


•K.  Goldberg  had  previously  found  a  different  proof  in  the  case  of  cyclic  fields. 
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follows  that  the  st-th  row  of  (att)  '(aa)  is  obtained  by  applying  the  per¬ 
mutation  Pj1  to  the  first  row.  Since  the  Galois  group  is  transitive,  every 
row  of  {aik)'(aik)  is  obtained  from  its  first  row  in  this  manner. 

We  will  now  discuss  the  case  of  cyclic  fields  only. 

It  was  shown  in  §2  that  for  n  =  2,  3,  4,  6  only  trivial  circulants  exist, 
hence  for  these  values  of  n  the  circulant  (aik)'(aik)  being  symmetric  is  the 
unit  matrix.  This  implies  that  ( ai]c )  is  an  orthogonal  matrix;  thus  each  row 
contains  only  one  element  (which  is  ±  1)  different  from  0  in  our  case  of 
rational  integral  elements.  Such  a  matrix  will  be  called  a  generalized 
permutation  matrix.  We  now  show 

theorem  4.  Let  F  be  a  cyclic  field  with  normal  basis.  There  exists  no 
basis  with  normal  discriminant  matrix  unless  it  is  a  normal  basis  or  a 
generalized  permutation  of  a  normal  basis. 

Theorem  4  follows  from 

theorem  4'.  A  unimodular  circulant  of  the  form  A  A'  where  A  is  a  matrix 
of  rational  integers,  is  equal  to  CC'  where  C  is  again  a  unimodular  circulant  of 
rational  integers. 

Proof:  We  find  it  convenient  to  prove  the  following  lemma  first: 

lemma  2.  Let  nbe  a  positive  integer,  and  suppose  that  M  is  a  generalized 
n-cycle  with  det  M  —  det  P.  Then  M  and  P  are  similar  elements  of  the  group 
of  generalized  permutation  matrices. 

Proof:  We  show  first  that  the  generalized  w-cycle 


0 

o, 

0  • 

•  •  0 

0 

0 

■ 

•  •  0 

0 

0 

0  • 

•  •  On- 

On 

0 

0  • 

■  ■  0 

is  similar  to  P,  if  alol  ••■<;„=  1  (this  is  the  necessary  and  sufficient  condition 
that  det  K  =  det  P).  Put  E  =  diag  (ej  ,  Ej  ,  •  •  •  ,  en)  where  e,  ±  1  and  the 
exact  values  of  e(  are  to  be  determined  later.  Then 
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We  seek  such  that 
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£1^2  —  al  i 
£>j£3  =  <7j  , 


en-l£n  —  an~l  > 

£ne  1  =  °n  • 

Since  •  •  •  an  —  1,  these  equations  have  the  solution 

£i  =  ±  1, 

s2  —  Ct-fa  , 
e3  =  OjCfiSi , 


en  —  OjOj  •  •  •  a^xh  • 

We  conclude  therefore  that  E~XPE  —  K,  so  that  K  and  P  are  similar 
elements  of  the  group  of  generalized  permutation  matrices. 

Suppose  now  that  M  is  a  generalized  tt-cycle  with  det  M  —  det  P. 
We  can  write  M  —  HM0  ,  where  H  is  a  diagonal  matrix  of  ±  l’s  and  M0 
is  an  n- cycle  in  the  usual  sense.  Hence  there  is  a  permutation  matrix  Q  such 
that  M0  =  Q'PQ,  and  so  M  =  HQ'PQ.  Since  Q'  is  a  permutation  matrix 
and  H  is  a  diagonal  matrix  of  ±  l’s  there  is  some  other  diagonal  matrix 
H0  of  ±  l’s  such  that  HQ'  =  Q'H0 .  Thus  M  —  Q'H0PQ,  which  implies 
det  H0  —  1,  since  det  M  =  det  P.  Hence  det  H0P  =  det  P,  and  so  H0P 
and  P  are  similar  elements  of  the  group  of  generalized  permutation  matrices. 
The  same  is  therefore  true  for  M  and  P,  and  Lemma  2  is  proved. 

We  return  now  to  Theorem  4'.  Let  A  be  a  unimodular  matrix  of  rational 
integers  of  order  n  such  that  A  A'  is  a  circulant.  Then  A  A'  commutes  with 
P,  so  that  ( PA){PA )'  =  AA'.  Thus  the  matrix  A~lPA  is  orthogonal,  and 
since  A  and  P  have  rational  integral  elements  and  are  unimodular,  A~XPA 
is  a  generalized  permutation  matrix.  Put  A~XPA  =  M,  so  that  PA  =  AM. 
We  remark  first  that  M  is  of  period  n,  since  M  is  similar  to  P.  Furthermore, 
write  A  as  the  matrix  of  its  row  vectors 

f 


Then  PA  =  AM  implies  that 

a2  =  a-xM, 

ag  =  OgM , 


<h  =  arM- 
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These  in  turn  imply  that 

<h 

axM 
A  =  a±M2 

a± M n_1 

We  now  show  that  M  is  a  generalized  w-cycle.  For,  assume  that  M  is 
the  product  of  at  least  two  disjoint  cycles.  We  shall  show  that  this  implies 
det  A  is  even,  which  is  a  contradiction  since  A  is  unimodular.  In  order 
to  show  this,  it  suffices  to  consider  A  modulo  2.  We  can  then  regard  M  as  a 
permutation  matrix.  Number  the  elements  so  that  the  cycles  are  (1  2  •  •  •  %) 
and  +1  ix  +  2  •  •  •  i2).  The  first  of  the  cycles  affects  only  the  first  ix 
columns  of  A ,  the  elements  in  each  row  of  this  submatrix  being  a  permuta¬ 
tion  of  ,  an  ,  •  •  •  ,  «!,■  .  The  determinant  of  A  therefore  is  divisible  by 
ax  —  -f-  Ojjj  -j-  •  •  •  -f-  au  ,  modulo  2.  Applying  the  same  argument  to  the 

second  cycle  shows  that  det  A  is  divisible  by  <r2  =  -f  «1<i+2  +  *  ••  +  #!<, , 
modulo  2,  and  applying  the  argument  to  the  product  of  the  cycles,  we  find 
that  det  A  is  also  divisible  by  ax  +  <r2 ,  modulo  2.  At  least  one  of  the  numbers 
a x ,  a2 ,  ax  +  a2  is  even.  This  implies  that  A  is  singular  modulo  2,  so  that  M 
can  not  be  the  product  of  more  than  one  disjoint  cycle. 

Since  M  and  P  are  of  equal  determinant,  Lemma  2  tells  us  that  there  is 
a  generalized  permutation  matrix  Q  such  that  M  —  Q'PQ.  Hence 

<hQ'Q  ci 

A==  “iQ'PQ  = 

a1Q'Pn~1Q  CjP”-1 

with  =  OiQ’. 

This  concludes  the  proof  of  Theorem  4'. 
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of  the  Second  Order  * 

WOLFGANG  WASOW 

University  of  California  at  Los  Angeles 


1.  Introduction 

This  article  is  concerned  with  the  asymptotic  calculation  of  the  solutions 
of  boundary  problems  for  nonlinear  differential  equations  of  the  second 
order  depending  in  a  singular  manner  on  a  small  parameter  e.  More  precisely, 
we  shall  study  real  differential  equations  of  the  type 

(1.1)  ey"  =  F(x,  y,y',  e) 
together  with  boundary  conditions  of  the  form 

(1.2)  y(x)  =  4  ,  ytf)  =  4 , 

where  a,  fi,  lx ,  l2  are  constants  independent  of  s.  The  precise  assumptions 
on  F  ( x ,  y,  y' ,  e)  will  be  found  below.  The  most  restrictive  one  is  that 
F{x,y,  y',  e)  must  be  linear  in  y’ ,  i.e.  that 

(1.3)  F(x,  y,  y',  e)  =  Fx{x,  y,  e)y'  +  F2{x,  y,  e). 

If  the  ’’reduced”  differential  equation  F(x,  y,  y',  0)  =  0  possesses  a 
solution  y0[x)  satisfying  one  of  the  two  boundary  conditions  (1.2)  it  is  natural 
to  ask  whether  the  full  problem  has  a  solution  y( x,  e )  that  tends  to  y0(x)  as 
e  ->  0  and  to  attempt  a  construction  of  y(x,  e)  by  means  of  some  perturba¬ 
tion  scheme. 

The  questions  concerning  the  existence,  uniqueness  and  convergence, 
as  e  -*  0,  of  solutions  of  (1.1)  and  (1.2)  have  been  dealt  with  by  E.  Codding- 
ton  and  N.  Levinson  [1]  and,  in  an  earlier  but  less  complete  paper,  by 
R.  von  Mises  [2],  N.  I.  Brish  [3]  using  a  result  of  M.  Nagumo  [4]  has  proved 
somewhat  more  general  results  when  F  is  of  the  form  F  (x,  y,  y' ,  e)  = 
Fx(x,  y)y'  +  F%{x,  y,  y')  with  F3(x,  y,  y')  bounded  for  all  y' .  As  pointed  out 
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in  [1]  the  boundary  value  problem  cannot  be  expected  to  possess  a  solution 
when  F  is  more  strongly  nonlinear  in  y' . 

The  present  paper  differs  from  these  investigations  in  that  it  contains  a 
constructional  scheme  for  the  representation  of  the  solution  by  means  of  converg¬ 
ent  and  asymptotic  series  which  reveal  its  analytic  structure.  The  calculation 
of  the  terms  of  these  series  can  be  effected  by  quadratures.  No  use  is  made 
of  the  results  of  [1],  so  that  our  results  furnish  as  a  side  product  a  new 
existence  proof.  The  construction  of  our  series  requires  analyticity  of 
F(x,  y,  y' ,  e)  with  respect  to  y,  which  is  more  restrictive  than  the  hypotheses 
of  [1]  and  [3],  as  far  as  y  is  concerned.  On  the  other  hand,  it  is  assumed  in 
[1]  and  [3]  that  F  is  independent  of  e.  It  is  not  difficult  to  extend  our  proof 
to  the  case  that  a,,  ft,  lx ,  l%  are  analytic  functions  of  e.  Generalizations  to 
certain  higher  order  differential  equations  are  very  likely  possible. 

As  in  most  perturbation  problems  in  which  a  small  parameter  multiplies 
the  highest  derivative,  there  occurs,  near  one  endpoint  of  the  interval  a 
phenomenon  of  non-uniform  convergence,  as  e  ->  0,  which  is  frequently 
called  a  boundary  layer  phenomenon,  because  the  boundary  layers  of  viscous 
flows  are  the  best  known  physical  applications  of  such  mathematical 
occurrences.  Usually,  two  analytic  representations  are  needed  in  such  prob¬ 
lems,  one  valid  in  the  boundary  layer,  the  other  one  outside,  and  to  calculate 
the  correct  transition  from  one  region  to  the  other  is  then  an  involved  and 
difficult  task.  For  linear  ordinary  differential  equations  asymptotic  expres¬ 
sions  valid  umformly  in  the  whole  closed  interval  under  consideration  have 
been  calculated  before  [5],  But  the  author  is  not  aware  of  such  representa¬ 
tions  for  nonlinear  equations.  There  exists,  however,  a  literature  on  the 
asymptotic  representation  of  solutions  of  initial  value  problems  for  non¬ 
linear  differential  equations  without  a  parameter,  which  uses  methods  re¬ 
lated  to  those  of  the  present  paper.  See,  for  instance,  the  articles  [6],  [7], 
[8]  by  Trjitzinsky  and  by  Malmquist.  There  are  also  points  of  contact  with 
the  investigations  of  I.  M.  Volk  [9]. 

2.  Statement  of  the  Problem 

To  fix  the  ideas  let  us  assume  that  e  >  0.  The  modifications  for  e  <  0 
are  obvious.  Our  perturbation  method  is  then  based  on  the  following  as¬ 
sumptions. 

assumption  A.  The  “reduced  differential  equation” 

(2.1)  F^x,  y,  0)y'  +  Ft{x,  y,  0)  =  0 
possesses  a  solution 

(2.2) 


y  =  y0  (*) 
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for  which 

(2.3)  y0(P)  —  4  and  F^x,  y0(x),  0)  <  0,  a  ^  x  g  /9. 

If  the  solution  y0(x)  satisfies  the  conditions 

y0(x)  =  4  and  F^x,  y0{x),  0)  >  0,  a  ^  x  ^  /9, 

instead  of  (2.3),  the  transformation  x*  =  —  x,  y*  =  —  y  changes  the  given 
problem  into  one  for  which  Assumption  A  is  true. 

assumption  B.  The  functions  Fj(x,  y,  e),  j  =  1,  2,  are  regular  analytic 
with  respect  to  y  and  e  and  of  class  C(2)  with  respect  to  x  in  a  region  R  of  the 
x,  y,  e-space  that  contains  in  its  interior  all  points  y  —  y0(x),  x  x  iS  /?,  e  =  0. 

The  assumption  of  analyticity  with  respect  to  e  can  be  replaced,  without 
additional  complications  by  the  milder  one  that  the  F^x,  y,  e )  admit  asymp¬ 
totic  expansions  in  powers  of  e. 

Without  loss  of  generality  we  may  assume 

(2.4)  a  =  0,  p=l, 

(2.5)  y0(x)  =  0, 
and 

(2.6)  F(x,  0,  y',  0)  =  -  y'. 

For,  the  transformation  of  the  variables  x,  y,  e  into  x*,  y*,  e*  by  means 
of  the  equations 

(2.7)  x*  =  y0(t),  0)  dt /  y0{t),  0)  dt, 

(2-8)  y*  =  y-  y0(»). 

(2.9)  e*  =  -e/  jyi(t,y0(t),0)dt, 

reduces  the  given  problem  to  one  having  these  properties,  as  can  be  seen 
by  direct  verification. 

With  these  assumptions  the  differential  equation  (1.1)  can  be  written 
in  the  form 

(2.10)  ey"  +  p(x,  e)y'  +  q{x,  e)  y 

=  ea(x,  e)  f(x,  y,  e)yy'  +  g(x,  y,  e)yi, 

where  all  coefficients  in  R  are  regular  analytic  with  respect  to  y  and  e  and 
in  class  C(2)  with  respect  to  x.  Moreover, 

(2.11)  p(x,  0)  =  1. 

Our  aim  is  to  construct  a  solution  of  the  differential  equation  (2.10)  that 
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assumes  the  prescribed  boundary  values 

(2.12)  y( 0,  e)  =  y°,  y(l,  e)  =  0 

and  tends  to  zero  in  0  <  x  5S  1,  as  e  -*■  0. 

3.  A  Preliminary  Transformation 

Our  construction  of  the  required  solution  consists  of  two  parts.  In  the  - 
first  part  we  calculate  a  particular  solution  yl(x,  e)  of  the  differential  equa¬ 
tion  (2.10)  that  passes  through  the  right  endpoint  x  =  1,  y  —  0,  but  not 
generally  through  the  left  endpoint  x  =  0,  y  —  y°,  and  which  has  the  pro¬ 
perty  of  being  in  the  whole  closed  interval  0  :£  x  ^  1  asymptotically  re¬ 
presentable  by  an  ordinary  power  series  in  e.  This  solution  has  no  ’’boundary 
layer”  for  small  e  at  x  =  0.  In  the  second  part  u  is  shown  that  the  required 
solution,  which  passes  also  through  x  —  0,  y  =  y°,  differs  from  y(x,  e)  by  a 
convergent  series  that  is  asymptotically  zero  except  in  an  infinitesimal  inter¬ 
val  near  x  =  0. 

In  the  present  section  the  differential  equation  will  be  subjected  to  a 
preparatory  transformation  consisting  in  the  subtraction  of  a  certain 
function  y*  (x,  e)  which  will  later  be  shown  to  be  asymptotically  equal  to 
&(*.  «)• 

We  start  from  the  observation  that  (2.10)  can  be  formally  satisfied  by 
a  series  of  the  form 

(3-1)  2>r(*)er, 

r— 1 

since  insertion  of  this  series  into  (2. 10)  followed  by  expansion  and  identifica¬ 
tion  of  like  powers  of  e  leads  to  recursion  formulas  of  the  form 

(3.2)  v[  +  q{x,  0)vx  =  a(x,  0), 

(3.3)  v'r  +  q(x,  0)vr  =  —  v”_x  -f  K(x,  vx ,  ■  •  • ,  vr_x ,  v[ ,  •  •  •,  v'^),  r>  1, 

which  can  be  solved  successively.  We  can— and  will— even  impose  the 
boundary  condition 

(3.4)  »,(!)  =  0,  r  —  1,2,-*-. 

The  series  so  determined  is  in  general  divergent  but  by  virtue  of  a 
theorem,  a  proof  of  which  can  e.g.  be  found  in  [10],  there  exists  a  function 
y*(x,  s)— in  fact,  infinitely  many  such  functions— analytic  in  e  and  of  class 
C(,)  with  respect  to  x,  in  0^z^l,0<£^e1  (ex  a  constant),  that  satisfies 
the  asymptotic  relation 

(3.5)  y?(*>  e)  vr(x)er, 

r-1 


I 


i 
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Also,  this  asymptotic  relation  may  be  twice  term-wise  differentiated 
with  respect  to  x.  (The  proof  in  [10]  deals  only  with  asymptotic  series  with 
constant  coefficients.  But  the  additional  facts  concerning  the  differentiable 
dependence  on  x  require  only  a  simple  corollary.)  Because  of  (3.4)  and  the 
construction  in  [10]  we  have  also 

(3.6)  yf(l,e)  =  0. 

The  function  yf  (x,  e)  is  in  general  not  a  solution  of  the  differential 
equation  (2.10),  but  with  its  help  we  shall  be  able  to  show  the  existence  of 
a  true  solution  yx{x,  e)  for  which 

(3.7)  y1{x,e)~yf(x,e)l  0  ^  x  ^  1,  0^e^elt 
and 

(3.8)  ffi(l,  s)  =  0. 

In  order  to  achieve  this  we  introduce 

(3.9)  z  =  y  —  yf{x,  e) 

as  new  dependent  variable  into  (2.10)  and  show  that  the  resulting  differential 
equation  for  z  possesses  a  solution  with 

(3.10)  z(x,  e)  ~  0,  Ogxgl,  0  <£  e  ex , 

(3.11)  2(1,  e)  =  0. 

The  transformation  results  in  a  differential  equation  of  the  form 

(3.12)  ez"  +  p(x,  e)z'  +  q[x,  e)z  =  a (x,  e)  +  h(x,  z,  z',  e) 
where 

(3.13)  «(*,  e)  =  f(x,  y f ,  e)  yf  yf  +  g[x,  yf ,  e)  yf  +  ea(x,  c) 

-  +  Pi*.  «)»*'  +  ?(*,  e)yf } 

and 

i 

(3.14)  h(x,  z,  z',  e)=f(x,  yf+z,  e){yf+z){yf'+z')  +g{x,  yf+z,  e){yf+z)* 

-  /(*.  y f,  e)yfyf  -  g(x,  yf,  «) yf. 

00 

Since  yf  (x,  e)  has  the  asymptotic  expansion  2  vrix)e'  which  satisfies 

r-1 

(2.10)  formally,  it  follows  that 

(3.16)  <x(x,  e)  t**  0,  0  e  . 

Expanding  the  right  member  of  (3.14)  with  respect  to  powers  of  z  and 
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considering  that  e  can  be  factored  from  the  right  hand  side  of  (3.5),  we  see 
that 

(3.16)  h(x,  z,  z\  e)  —  e(n(x,  e)z'  +  q(x,  e)z)  +  h*(x,  z,  z’,  e) 
with 

(3.17)  h*(x,  z,  z',  e)  =  f*(x,  z,  e)zz'  +  g*(x,  z,  e)z2. 

Here  n,  q,  f*  and  g*  have  asymptotic  power  series  expansions  in  e  and 
are  in  class  C<2>  with  respect  to  x.  Moreover,  f*  and  g*  are  analytic  in  z  at 
2  =  0.  Combination  of  (3.16)  and  (3.17)  with  (3.12)  leads  to  the  form 

(3.18)  ez"  +  P*{%>  £)z'  +  ?*(*,  e)z  =  »l(x,  e )  +  h*(x,  z,  z',  e) 

for  our  differential  equation,  which  is  the  basis  of  our  further  analysis.  We 
observe  that 

(3.19)  p*(x,  0)  =  p(x,  0)  =  1,  q*{x,  0)  =  q(x,  0), 

because  of  the  factor  e  in  the  right  member  of  (3.16). 

Our  subsequent  arguments  are  based  on  the  asymptotic  properties  of 
linear  differential  equations  whose  homogeneous  part  is  the  left  member  of 
(3.18).  We  therefore  collect  in  the  next  section  some  lemmas  concerning 
such  equations. 


4.  Lemmas  on  Linear  Equations 

In  the  following  the  symbol  [/(*)]  will  be  used  to  designate  functions 
of  the  form 

[/(*)]  =  f(x)  +  eE(x,  e), 

where  E  (x,  s)  and  its  first  derivative  with  respect  to  x  are  bounded  in 
0  £  *  £  1,  0  ^  e  ^  . 

The  first  lemma  below  is  a  direct  application  of  the  standard  asymptotic 
theory  for  linear  differential  equations  (see  [11]  and  [5]). 

LEMMA  4.4.  The  “variational  differential  equation” 

(4.1)  eV”  +  p*(x,  e)V’  +  q*(x,  e)V  =  0 
possesses  two  solutions  of  the  form 

(4.2)  V^x,  e)  =  e~«*[v*(x)l  Vt(x,  e)  =  [»(*)] 
where 

(4.3)  ct  =  e-1, 

(4.4)  v(x)  =  exp  ||*  q*  [t,  0)  <ft}, 
and 


v*(x)  —  v-^x)  expj  — J*w(^,  0)<ft}. 
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corollary.  In  terms  of  its  boundary  values  V  (0,  s),  V  (1,  e)  every  solution 
V  (x,  s)  of  (4.1)  can  be  written  in  the  form 


V(x,  e)  =  Ax(a)Vx{x,  e)  +  A2{e)V2(x,  »>, 


where 


A1(e)  =  V (0,  e)[l]  —  V (1,  «)[®(1)-*], 

A2(e)  =  -  V (0,  e)e-°[v*(l)/v(l)]  +  F(l,  *)[»(1)"‘]. 

lemma  4.2.  The  solution  WxJx,  e)  of  the  differential  equation 

(4.5) 
with 

(4.6) 

is  of  the  form 

Wx  (x,  e)  =  f*  K(x,  t,  e)<p(t,  e)  dt 
■I  x0 


eW"  +  p*(x,  e)W'  +  q*{x,  e)W  =  <p(x,  e) 
WXt(x0  ,  e)  =  ,  e)  =  0 


(4.7) 

with 


(4.8) 

that 


~v(x)- 

U(/)J 

— i 

'♦si 

* 

Jj 

» 

K(x,  t,  e) 

Proof:  An  application  of  the  method  of  variation  of  parameters  shows 

,  s  rw  e)Vi(x,  ej-VAx,  e)V,(t,  e)  /t  ^ 
w-{x' e>  =  *> * 

From  this  formula  the  statement  of  the  lemma  is  obtained  by  substitution 
of  (4.2)  followed  by  some  elementary  manipulation. 

lemma  4.3.  The  solution  of  the  differential  equation 

(4.9)  eW"  +  p*(x,  e)W'  +  q*(x,  c)W  =  e-*°x<p{x,  e) 

for  which 

W(0,  e)  =  W(  1,  «)  =  0 


(4.10) 

is  of  the  form 

(4.11) 
with 

(4.12) 


W(x,  e)  —  ee  °xco( x,  e) 


a>w(x,  e)  |  sup  \<p(x,  e)|, 

OSfSl 

0£eSe, 


v  =  0,  1. 


The  constant  c  is  independent  of  <p(x,  e)  and  e. 


* 

a 


100 


WOLFGANG  WASOW 


Proof:  The  solution  V(x,  e)  of  (4.1)  that  assumes  the  boundary  values 

(4.13)  F(0,  e)  =  1,  7(1,  e)  =  0 

is  of  the  form 


(4.14)  V(x,  s)  =  [»*(1)» (*)/»(  1)]}. 

This  follows  directly  from  the  corollary  to  Lemma  4.1.  Now, 


(4.15)  W(x,  e)  =  Wx(x,  e)  -  Wx( 0,  e)V(x,  e), 

and,  by  Lemma  4.2, 


(4..6)  W, (*,  .)  =  .-|[{|}|]  [£f 


<p((,  e)  dt. 


Lemma  4.3  is  a  simple  consequence  of  the  last  three  formulas. 


5.  A  Particular  Solution  without  Boundary  Layer 

From  Lemmas  4.1  and  4.2  it  follows  that,  if  the  non-linear  integral 
equation 

z(x,  e)  =  [* K(x,  t,  e)  (<x(f,  e)  +  h*(t,  z(t,  e) ,  z'(t,  e),  e)}dt 
(5-1)  J° 

-  Jotf(l,  t,  e)  {«(/.  e)  +  h*(t,  z(t,  e),  z'(t,  e),  e)}dt  ■  7^(1,  e)  Vt(x,  e) 

possesses  a  solution,  then  it  is  also  a  solution  of  the  differential  equation 
(3.18).  Moreover, 

(5.2)  *(1,  e)  =  0. 

We  shall  now  construct  such  a  solution  by  Picard’s  iteration  method 
and  show  that  it  is  asymptotically  equal  to  zero. 

Let 


(5.3)  zQ(x,  e)  =  0, 

and  define  zn+1(x,  e),  n  =  0,  1,  —  ,  as  the  value  of  the  right  member  of 
(5.1)  when  z  and  z'  are  replaced  by  zn  and  z'n ,  respectively.  Because  of 

(3.15)  and  (3.17)  we  have  from  (5.1)  that 

z^{x,  e)  ~  0, 


(5.4) 


v  =  0,  1,  0^a;^l, 
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Generally, 

4+1  (*.  4  “  Zr](X *  £) 

=  t,  e){h*(t,  zr(t,  e),  z'r(t,  e),  e)  —  h*(t,  zr_x{t,  e ),  4-i  (t,  e),  e)}dt 

(5.5)  J° 

—  \X K{  1,  t,  e)  {h*(t,  zr(t,  e),  z'T{l,  e),  s) 

—  h*(t,  z^t,  e),  Zr_x(t,  fi),  e)}dt-  V^{\,  e)  V£} (x,  e),  v  =  0,  1. 

The  inequalities 

(5.6)  \z?(x,  e)|  ^  Mi? ,  v  —  0,  1,  0  iS  x  1,  0  <  s  5S  ex 

where  N  is  an  arbitrary  positive  integer  and  M  a  constant  depending  on  N 
and  gj ,  are  certainly  true  for  r  —  1,  in  view  of  (5.4).  Let  us  assume,  for  a 
proof  by  induction,  that  they  have  been  proved  for  r  ^Ln.  It  follows  then 
from  (3.17)  that— at  least  for  £j  so  small  that  zr(x,  s)  and  z'r(x,  e)  lie  in  the 
domain  of  regularity  of  h*  (t,  z,  z' ,  e)  — 

zr(t,  e),  z'r(t,  e),  e)  -  h*(t,  z^t,  e),  4_i(*,  4*  4  | 

(5.7)  ^  Ml£( \zr{t,  e)  -  zr_i (t,  e)  |  +  j z'r(t,  e)  -  4_i(<>  e)  \ ) , 

r  <Ln,  Ogigl,  0  ^  e  ^  . 

(The  factor  e  might  even  be  replaced  by  eN.)  M1  depends  on  M.  For  abbre¬ 
viation  we  introduce  the  notation 


(5.8)  Cr  =  sup  ( I zr(t,  e)  -  zr_i(f,  e)|  +  |4(*>  4  ~  4-i(<*  41)- 

0S«S1 

OSeitx 

Insertion  of  (5.7)  and  (4.8)  into  (5.5)  followed  by  summation  over  v leads  to 

(5.9)  Cr+i  ^  ,  r  ^  n,  0  ^  e  ex , 

with  another  constant  M, ,  and  hence  to 

-(») 


1 4+i  (**  41  = 


(5.10) 


fl+1 

2  4’)(*.  4  -  4-i  (*. 4 

r-1 


«+l  w+1 

f"l  r— 1 


eMjCi 

1  -  ’ 


v  =  0,  1,  O^x^l,  0  ^  e  ^  . 


But 


&  ^  4Me*, 

by  virtue  of  (5.3)  and  (5.5);  and  (5.10)  becomes,  therefore, 

I4’ii(*<  41  ^  1 - TF-  M*N>  *  =  0,  1,  0^*:gl, 

1  —  MfE 
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The  constants  M,  Ma  depend  on  ,  but  they  are  bounded  if  sx  is  small. 
For  et  sufficiently  small,  therefore,  4eMs/(l  —  Mae)  is  less  than  unity,  and 
(5.6)  is  proved  for  r  =  n  -f  1  also,  i.e.  it  must  be  generally  true. 

For  this  value  of  the  factor  M^s  in  (5.9)  is  less  than  unity  as  well. 
From  this  fact  the  convergence  of  the  sequence  zn(x,  s)  to  a  solution  z(x,  e)  of 
the  integral  equation  follows  by  the  usual  argument  of  Picard’s  iteration 
method.  At  the  same  time  one  proves,  that  z'n(x,  e)  tends  to  z'(x,  e).  Finally, 
letting  r  tend  to  infinity  in  (5.6)  we  see  that,  since  N  is  arbitrarily  large, 

(5.11)  zw(x,  e)~0,  r=0,  1,  0  ^  ^  1,  0  ^  e  ^  1. 

The  second  derivative  z"(x,  e)  is  then  also  asymptotically  zero,  in  con¬ 
sequence  of  the  differential  equation  itself.  By  repeated  differentiation  of  the 
differential  equation  the  formula  (5.11)  can  be  extended  to  as  many  deriva¬ 
tives  as  exist  and  are  continuous. 

The  theorem  below  summarizes  these  results. 

theorem  1.  The  differential  equation 

ey"  +  p{x,  e)y'  +  q(x,  e)y  =  ea(x,  e)  +  /( x,  y,  e)yy'  -f  g{x,  y,  e.)y\ 

with  f)(x,  0)  =  1,  possesses,  for  0  <i  e  ^  ,  solutions  that  admit  in  0  ^  x  ^  1 

a  uniformly  valid  asymptotic  representation  by  a  series  in  powers  of  e,  as 
0.  Among  these  solutions  there  is  one  y^x,  e),  with  yx(  1,  e)  —  0.  The 
asymptotic  series  may  be  twice  termwise  differentiated  with  respect  to  x. 

We  shall  see  later  that  yx(x,  e)  is  not  the  only  such  solution,  but  it  is 
clear  that  all  solutions  that  are  asymptotically  represented  by  a  power  series 
in  0  ig  x  iS  1  and  vanish  at  x  =  1  can  at  most  differ  by  functions  that  are 
asymptotically  zero,  since  the  coefficients  of  the  power  series  must  satisfy 
the  differential  equations  (3.2),  (3.3). 

6.  Construction  of  the  Boundary  Layer  Correction 

In  Section  3  the  function  y*  (x,  e)  could  have  been  any  function  satis¬ 
fying  the  asymptotic  relation  (3.5).  We  may  now  take  as  yf  (x,  e)  in  the 
transformation  (3.9)  the  particular  solution  yx(x,  e)  constructed  in  the  pre¬ 
ceding  section.  The  transformation 

z  =  y  —  «) 

leads  then,  as  before,  to  the  differential  equation  (3.18).  Since  y^x,  e)  is  a 
true  solution  of  the  differential  equation  (2.10)  we  see  from  (3.13)  that  we 
now  have 

(6.1)  a(x,  e)  s=  0 

instead  of  the  weaker  asymptotic  relation  (3.15).  To  the  required  solution  of 
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(2.10)  with  boundary  values  (0,  i/0)  and  (1,  0)  there  corresponds  a  solution 
of  (3.18)  with 

(6.2)  2(0,  e)  =  y°  —  yt{ 0,  e),  z{  1,  e)  =  0. 

The  existence  and  analytic  structure  of  z(x,  e)  could  probably  be  in¬ 
vestigated  by  Picard’s  method  making  use  of  Green’s  function.  This  method 
would  have  the  advantage  of  being  applicable  even  when  the  differential 
equation  is  not  analytic  in  y.  But  since  the  necessary  calculations  seem  to  be 
rather  involved,  and  since  we  are  aiming  at  an  analytic  expression  in  series 
form  we  shall  find  the  solution  instead  in  the  form  of  a  power  series  in  its 
initial  value  2(0,  e)  and  prove  the  convergence  by  the  method  of  dominating 
series. 

Let  us  set,  for  abbreviation. 


(6.3)  y°  —  yt{ 0,  e)  =  y 
and  insert  for  2  in  (3.18)  formally  the  series 

(6.4)  *  =  2 

f— 1 

where  ur  =  ur(x,  e)  are  functions  to  be  determined.  By  identification  of  like 
powers  of  /i  and  by  using  (6.1),  we  obtain  the  recursive  differential  equations 

(6.5)  eu"  +  p*(x,  e )u[  -f-  q*( x,  e)%  =  0, 

(6.6)  eu"  +  p*(x,  e)u'T  -f  q*(x,  e)ur  =  hr{ x,  ,  •  •  •,  Vu ,  u[  ,  •  •  •,  ur_x ,  e), 

r  >  1, 

where  the  functions  hr  have  the  following  properties: 

1)  hr  is  a  polynomial  in  the  and  u'} ,  j  —  1,  •  •  • ,  r  —  1,  linear  in  the 
u'j  combined  and  without  constant  or  linear  terms  in  the  u}  and  u}  combined. 

2)  The  coefficients  of  this  polynomial  possess  asymptotic  power  series 
in  e  each  term  of  which  is  a  function  of  x  of  class  C(2). 

If  the  ur  are  determined  from  these  differential  equations,  and  if  the 
resulting  series  (6.4),  as  well  as  its  termwise  derivative  with  respect  to  x 
converge  uniformly  and  absolutely  in  0  gi  x  ^  1,  then  it  represents  there  a 
solution  of  the  differential  equation  (3.18)  with  a(x,  e)  =  0.  In  fact,  if  we 
multiply  the  equations  (6.5)  and  (6.6)  by  /a,  and  /ur,  respectively,  and  sum 
over  r,  the  right  member  is,  in  view  of  (6.4),  precisely  the  function 
h*(x,  z(x,  e),  2'(x,  e),  e),  and  the  last  two  terms  in  the  left  member  become 

p*{x,  e)z'  -j-  q*{x,  e)z.  The  first  term  of  the  left  member,  i.e.  e  2  ur'ur,  must 

f“l 

therefore  also  be  uniformly  convergent  and,  hence,  equals  ez" . 
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The  boundary  conditions  (6.2)  are  satisfied  if  we  require  that 

(6.7)  Wl(0,  e)  =  1,  %( 1,  e)  =  0, 

(6.8)  ur( 0,  e)  =  ur(  1,  e)  =  0,  r  >  1. 

7.  Convergence  Proof 

The  asymptotic  form  of  u^x,  e)  is  immediately  determined  from  (6.5) 
and  (6.7)  with  the  help  of  (4.13)  and  (4.14).  We  find 

u^x.e)  =  *—{[#*(»)]  -  e^[v*(l)v{x)lv(l)]} 

=  x,  e),  0  ^  e  ^  e, 


(7.1) 

where 

(7.2) 


Wi{x,  «)| 


v  =  0,  1, 


%  being  a  constant  independent  of  e. 

The  convergence  of  the  series  (6.4)  will  be  proved  by  constructing  a 
00 

dominating  series  J  ur  with  constant  positive  coefficients.  As  a  first  step 

r— 1 

we  construct  a  function  h(z,z')  which  dominates  h*(x,  z,  z' ,  e).  Let  the 
constants  %  >  0  and  s,  be  so  small  that  the  two  series 

(7.3)  /*  (x,  z,  e)  =  2  y,(x,  e)z\  g*  (x,  z,  e)  =  2  <*•(*»  £)** 

i—O  •— 0 

for  the  functions  in  (3.17)  converge  uniformly  and  absolutely  for 

(7.4)  O^z^l,  |z|  ^  zx  ,  0  ^  e  ^  Sl  . 

If  C  is  a  constant  such  that 

(7.5)  \f*{*,z,e)  \&C,  \g*{x,z,e)  |  gC, 

in  that  domain,  then  both  functions  are  dominated  by  the  series 

ci(-V 

.tbVzi/ 


which  represents  the  function 


Hence 


Kr 


(7.6)  %*')  =  c(l--^)  V  +  z*) 

dominates  h*(x,  z,  z' ,  e),  in  the  sense  that  every  term  in  the  series  for  h* 


5 
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in  powers  of  z  and  z'  has  a  coefficient  whose  absolute  value  in  0  ^  x  ^  1, 
0  g  e  ^  has  the  corresponding  coefficient  in  the  expansion  of  h  as  upper 
bound. 

Next,  functions  hr  dominating  the  functions  hr  introduced  in  (6.6)  are 

oo  oo 

defined  as  follows:  Insert  into  h(z,  z')  the  formal  series  2  uTf*r,  2  vr/zT  for 

r— 1  r— 1 

z  and  z' ,  respectively,  and  rearrange  according  to  powers  of  fi.  The  coefficient 
of  fiT,  which  will  be  denoted  by  >  ^r-i :°i  > ' ' ' ,  vr-i)  is  a  poly¬ 

nomial  in  the  indicated  arguments,  linear  in  the  v,  with  positive  coefficients 
that  are  upper  bounds,  in  0  sS  x  ^  1,  0  sS  £  ^ ,  for  the  moduli  of  the 

corresponding  coefficients  of  the  polynomial  hr(x,  z1 ,  •  •  • ,  zr_x 

Cx  .  *)• 

We  now  proceed  to  define  a  sequence  4,  that  dominates  uT .  It  has 
already  been  proved— cf.  (7.1),  (7.2)— that 

(7.7)  | {x,  fi) |  ;£  t ~axux  ,  \u[  ( x ,  e)|  ^  oe~axu1  ,  0  fg  x  ^  1,  0  ^  e  5S  c1. 

We  shall  now  show  generally  that 

(7.8)  |  ur(x,  fi)  |  ^  e-°*4r  ,  \ur(x,  fi)  |  ^  ae~axur  ,  O^x^l,  0^e^£l 
where  the  ur  are  determined  from  the  recursive  relation 

(7.9)  ur  =  chr(uy ,  •  •  •  ,  ur_i  j  «| ,  •  •  • ,  wr_x).  r  —  2,  3,  •  •  •  . 

Consider  first  h2(x,  ullu1,  e).  This  function  is  a  polynomial  in  ur  and 
tq  without  constant  or  linear  terms  and  linear  in  u[ .  Hence,  we  obtain  from 
(7.7)  and  the  definition  of  h2( u 1 ,  tq)  the  inequality 

I K(x,  %  ,  u[,e)\  <;  ae~2axJi2(u1 ,  jq) 

and  therefore,  by  virtue  of  Lemma  4.3  and  equations  (7.9), 

(7.10)  | «£>  (*,  fi)  |  ^  c<fe~a%{ul ,  *q)  =  (fe~°xu2  . 

The  completion  of  the  inductive  proof  of  the  formulas  (7.8)  is  straight¬ 
forward  and  is  therefore  omitted. 

00 

If  we  can  show  the  convergence  of  the  series  J  urfir,  the  uniform  and 

r— 1 

absolute  convergence  of  the  series  (6.4)  will  follow.  Now,  the  equations 

00 

(7.9)  are  also  obtained,  if  the  series  4r(«r  is  inserted  for  |  into  the  equation 

r—  1 

(7.11)  £  -  ch(l  f)  -  ^iq  =  0 

and  the  coefficients  of  /S  are  set  equal  to  zero.  The  equation  (7.11)  is 
satisfied  for  £  =  fi  =  0,  and  the  partial  derivative  of  the  left  member  with 
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respect  to  f  is  different  from  zero  at  £  =  p  =  0.  Therefore  (7.11)  defines  { 
as  a  regular  analxtic  function  of  p  in  a  neighborhood  of  p  —  0.  The  series 


2  wryitr  represents  this  function,  hence  it  converges  for  sufficiently  small 

r-l 


H  say  H  ^  /*o  • 

00  /  00 

The  derived  series  ^  wr(a;,  s)//  is  dominated  by  J  «raJMr.  The  series 


oo  f*l  r— 1 

^eu'r(z,  e)pr  is  therefore  also  uniformly  convergent  in  the  same  interval. 


T—l 

It  then  follows,  from  the  differential  equation  (1.1)  itself  and  from  the 


formal  construction  of  the  uT(x,  e)  that  the  series  2  e2u'r'(x,  e)pr  converges 

r—  2 

uniformly  as  well.  (The  factors  e  and  ea  in  these  series  are  necessary,  if  we 
wish  our  statements  to  hold  for  x  =  0  in  the  closed  interval  0  <;  e  5S  . 

The  value  of  eu'T(x,  e)  at  s  =  0  is  of  course  to  be  defined  by  lim  eu'T(x,  e).) 

£~*0 

co 

Because  of  the  factor  e~ax  in  (7.8)  the  values  of  the  series  2  ur]  (x>  £)<“f> 

r— 1 

v  =  0,  1,  tend  exponentially  to  zero,  in  every  interval  0  <  x1  ^  x  ^  1. 

In  the  theorem  below  these  results  are  summarized  in  terms  of  the 
original  differential  equation  (1.1).  For  convenience  the  notation  has  been 
changed  in  a  few  points. 


theorem  2.  If  the  differential  equation 


ey"  =  Fi{x,  y,  e)y'  +  Fz(x,  y,  e),  e  >  0 

satisfies  Assumptions  A  and  B  of  Section  2,  then  there  exist  constants  ex  >  0, 
px  >  0  such  that  the  differential  equation  possesses,  for  0  <  e  ^  e1 , 

14  —  2/o(0)l  sS/q,  in  a.  r^x  ^  fi  a  solution  satisfying  the  boundary  condition 

y(<t,  e)  =  4  ,  y(P,  e)  =  l2 . 


This  solution  is  representable  by  a  uniformly  and  absolutely  convergent 
series  of  the  form 

(7.12)  y(x,  «)  =  2  «.-(*-  e)t*T- 

r-0 


The  convergence  is  uniform  with  respect  to  e,  p  and  x  for  0  e  2S  % , 

l4-yo(*)l  ^th>  a^z^p. 

The  terms  of  this  series  have  the  following  properties: 

1)  w0  (z,  0)  is  the  given  solution  y0(x)  of  the  reduced  equation  which  satis¬ 
fies  the  right  hand  boundary  condition  u0(p,  0 )  =  l2  ■ 

2)  u0(x,  e)  has  an  asymptotic  series  expansion  in  powers  of  e.  This  series 
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is  a  formal  solution  of  the  differential  equation  and  satisfies  the  right  hand  bound¬ 
ary  condition. 

3)  fi  =  l1  —  m0(«,  e). 

4)  The  functions  uT(x,  e),  r  >  0  are  of  the  form 

(7.13)  ur(x,  e)  =  exp  |  —  a  J*  Fx(t,  y0(t),  0)  dt  J  <or(x,  e) 

where  wr(x,  e)  is  bounded. 

The  series  (7.12)  may  be  twice  termwise  differentiated  with  respect  to  x, 
for  0  <  e  5S  Cj.  The  series  ^e'u^  (x,  e)[T ,  v  =  0,  1,  2  converge  uniformly  in 
the  full  range  0  5S  e  —  t/0(a)|  ^  fzx  ,  «  5S  x  ^  /?. 

REMARKS. 

1)  The  boundary  layer.  When  /z  =  0  then  y^(x,  e)  tends  to  y^  (x) 
uniformly  in  the  whole  interval  0  x  :g  1,  for  v  =  0,  1,  2, .  •  • ,  n,  where 
n  —  2  is  the  order  of  differentiability  of  F(x,  y,  y',  e)  with  respect  to  x. 
In  other  words,  there  is  then  no  boundary  layer  phenomenon.  This  is  a  little 
surprising,  for  one  often  hears  the  remark  that  the  boundary  layer  occurs, 
because  the  solutions  of  the  reduced  first  order  equation  cannot  satisfy  both 
prescribed  boundary  conditions.  From  this  view  point  one  would  expect  the 
boundary  layer  to  disappear  when  the  solution  y0(x)  —  u0(x,  0)  of  the  re¬ 
duced  equation  happens  to  pass  through  the  left  boundary  point  as  well, 
i.e.  when  lx  —  u0( a,  0)  =  0.  However,  y  —  0  means  that  lx  —  u0 (a,  e)  =  0 
or  lx  —  w0(a,  0)  =  0(e).  Thus,  when  both  boundary  conditions  are  satisfied 
by  y0(x),  Theorem  2  shows  that 

y{x,  e)  =y0(x)  +  0(e)  +  exp  {  -  o  J*  Fx{t,  y0(t),  0)  dt}  0(e). 

Therefore,  there  is  no  boundary  layer  in  y(x,  e)  itself,  but— unless,  excep¬ 
tionally,  the  0(e)  in  the  third  term  is  0(e2)— it  is  still  true  that 

limy' (a,  e)  =£y'0  (<*■). 

E—>0 

2)  The  leading  perturbation  terms.  The  calculation  of  the  series  (7.12)  re¬ 
quires  the  successive  solution  of  sequences  of  linear  differential  equations 
provided  y0(x)  is  given.  The  terms  of  the  asymptotic  series  for  u0(x,  e )  can 
be  determined  by  quadratures.  In  view  of  formula  (7.1)  the  same  is  true  of 
ux(x,  0)eax.  We  give  the  leading  terms  explicitly:  Let 

Mx)  =  F*(x>  yo(«).  «»(*).  °). 

p2(x)  =  Fx(x,  y0(x),  0), 

«(*)  =  Fe(x>yo(x)>y'o(x),Q). 

vx(x)  =  V  exp  j  P  pxp2l  dt  |  a  dt. 
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(The  subscripts  y,  s  indicate  partial  differentiations.)  Then  the  solution  of 
(1.1)  satisfying  (1.2)  is  approximately  equal  to 

y0(x)  +  ev^x)  +  (4  -  y0( a)  -  eMa))  expj  J*  (apa  +  ft#4)  *}  . 

We  omit  the  straightforward  calculation  that  leads  to  this  formula. 
The  error  is  0(e2)  in  every  fixed  subinterval  0  <  aq  ^  a:  ^  1.  In  the  whole 
interval  0  5S  x  ^  1  the  error  is  0((4  —  t/0(a)) 2  +  e2). 

3)  The  higher  perturbation  terms.  If  the  asymptotic  series  for  uQ(x,  e)  has 
been  calculated  to  within  terms  of  order  0(eN)— which,  by  virtue  of  the  for¬ 
mulas  (3.2),  (3.3)  and  (3.4),  requires  only  quadratures,  no  matter  how  large 
N— then  the  coefficients  of  the  differential  equations  (6.5),  (6.6)  can  also 
be  determined  with  the  same  accuracy.  Solution  of  these  linear  differential 
equations  yields  the  functions  cor(x,  e)  of  (7.13)  to  within  terms  of  the  same 
order.  If  F(x,  y,  y',  e)  is  analytic  in  x,  then  the  functions  [y*(x)]  and  (y(x)] 
of  (4.2)  possess  asymptotic  power  series  in  e,  [12],  and  the  approximate 
solution  of  (6.5)  and  (6.6)  to  any  desired  accuracy  can  also  be  performed  by 
quadratures. 


8.  Reordering  of  the  Series 


If  F ( x ,  y,  y' ,  e)  is  analytic  in  x  also,  our  results  can  be  made  more  explicit 
on  the  basis  of  the  following  lemma. 

LEMMA  8.1.  Let  tp{z,  e)  be  a  function  for  which,  as  e  ->  +  0  an  asymptotic 
expansion  of  the  form 

00 

<p(z,  e)  ~  ^  e>r (z) 

r*»0 

is  uniformly  valid  and  indefinitely  termwise  differentiable  with  respect  to  z,  in 
0  ^  z  ^  1;  then 

(8.1)  \\ma«-*cp(t,  e)  dt  =  exp{z,  e)  -  ee-™x(e), 

where  y>(z,  e)  also  possesses  an  asymptotic  expansion,  whose  leading  term  is 
(p0{z).  Moreover,  X(e)  ~  y>(0,  e). 

Proof:  By  means  of  »  +  1  integrations  by  parts  one  finds  that 


(8.2) 


J*  erao(‘-‘V(*,  e)dt 

-si  <-»•(=)>*■■  i  -  s'-  e) 

+  (-  e)dt. 
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One  more  integration  by  parts  shows  that  the  last  term  of  the  right 
member  is  0(en+2),  uniformly  in  0  ^  z  zS.  1.  Now,  by  assumption, 

(8.3)  ?><•>  (z,  e)  =  £  £><*>  (z)  +  £"+!£'•>  (z,  e), 

r-0 

where  each  (z,  e)  is  uniformly  bounded  in  0  e  ^  Cj ,  0  ^  ^  1.  In¬ 

serting  (8.3)  into  (8.2)  and  rearranging  according  to  powers  of  e  we  obtain 
a  relation  of  the  form 

(8.4)  f*  emo<‘-*V(<,  e)dt  =  e  £  £>,(*)  -  ee -mo*  J  e>,(0)  +  0(e"+1) , 

0  1—0  v«0 

with  y0 (z)  =  <p0(z).  Let  y>(z,  e)  denote  some  function  for  which 

(8.5)  y>(z,  e)  ~  f  £>»(«). 

i—0 

uniformly  in  0  ^  z  ^  1.  The  existence  of  such  a  function  is  assured  by  the 
theorem  already  used  in  formula  (3.5).  Since  the  sums  in  the  right  member 
of  (8.4)  differ  from  ip(z,  e)  or  y>(0,  e)  by  terms  of  the  order  0(e"+1)  and  n  is 
arbitrary,  the  proof  of  the  lemma  is  at  hand. 

With  the  help  of  this  lemma  and  those  of  Section  4  the  analytic  form 
of  the  functions  ur{x,  e)  can  be  studied  more  explicitly.  The  reasoning  con¬ 
sists  in  a  rather  simple  elaboration  of  the  previous  discussion  of  these 
functions.  We  content  ourselves  with  the  statement  of  the  resulting  formula 

(8.6)  ur  {x,  e)  =  2  <Pr,  {*’  E)e~’°x  > 

f— 0 

(provided  j*  Fx(t,  y0{t),  0 )dt  =  x)  in  which  all  <prt(x,  e)  possess  asymptotic 
power  series  in  e.  In  particular  it  turns  out  that 

(8.7)  tp^x,  e)~0,  r>  0. 

00 

The  formula  (8.6)  suggests  a  rearrangement  of  the  series  ^  ur(x,  e)fiT 

oo  r— 0 

according  to  powers  of  e~ax.  This  is  permissible,  if  the  series  (x,  e)fir  with 

r—0 

*?(*»  «)  =  £  \<Pr,{x>  ®)l  e~KX 

converges.  That  this  is  so  requires  again  only  a  somewhat  detailed  retracing 

00 

of  the  convergence  proof  for  2  urix>  e)pT>  which  will  be  omitted. 

r— 0 

The  rearrangement  leads  to  a  series  of  the  form 

m—0 
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Because  of  (8.7)  we  have 

a0(x,  s,  4)  ~  u0(x,  e). 

The  functions  am(x,  e,  4)  are  of  the  form 

00 

4)  =  2  (4  -  «o(«.  «)) '?«.(*.  «) 

m 

with  coefficients  gm  (x,  e)  analytic  in  x  and  possessing  asymptotic  series  in  e. 

oo 

Analogous  arguments  can  be  applied  to  the  series  2 e'ut) (x>  £)/-<r. 

r=l 

(r  =  1,  2).  The  theorem  below  states  these  results  in  the  notation  of  the 
general  case. 

theorem  3.  If  Fx(x,  y,  e)y  +  Fa(x,  y,  e)  is  analytic  in  x,  in  0  5S  x  ^  1, 
in  addition  to  satisfying  Assumptions  A  and  B,  then  the  solution  (7.12)  of  the 
boundary  value  problem  formulated  in  Theorem  2  is  of  the  form 

(8.8)  y{x,  e)  =  '£cm(x,  4  ,  e)  exp  {mo  Fx{t,  y0{t),  0)  dt\. 

m-0  l  J  a  J 

The  series  converges  uniformly  and  absolutely  for  0  ^  x  5S  1,  0  :£  s  ^  ex , 
\l1  —  y0{x)\<fi1.  The  functions  cm(x,  4 .  e)  are  there  regular  analytic  in  x  and 
in  4  and  possess  asymptotic  series  in  powers  of  e.  Moreover 

c0(x,  4  ,  e)  ~  u0(x,  e), 

and 

cm(x>  k>e)~  (4  ~  «o(*.  e))mcZ{x>  4  >  «). 

where  c*(x,  4  ,  e)  has  the  same  regularity  properties  as  cm(x,  lt,  s). 

The  series 

(8.9)  f  < b' ^cm(a:,  4  ,  «)  expjmo J*  *i(f,  y0{t),  0)  v  =  1,  2, 

converge  uniformly  in  the  same  range. 

9.  Stretching  of  the  Boundary  Layer 

Although  the  series  representations  in  Theorems  2  and  3  are  valid  in  the 
whole  interval  a  5S  a;  SjJ  /?,  they  have  certain  disadvantages  for  the  approxi¬ 
mate  calculation  of  the  solution  inside  the  boundary  layer.  The  approximate 
formula  derived  in  Remark  2  of  Section  7,  for  instance,  is  satisfactory  only, 
if  4  —  y0(«)  i-s  small,  more  precisely,  if  it  is  of  the  order  0(e).  In  the  hydro- 
dynamic  applications,  where  such  a  condition  is  usually  not  satisfied,  approx¬ 
imate  solutions  in  the  boundary  layer  are  in  general  obtained  by  a  „stretch- 
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ing”  of  the  independent  variable.  In  our  example  a  natural  stretching  trans 
formation  consists  in  setting 


(9.1)  x  =  a  +  ze,  y(x,  e)  =  y(a  +  ze,  e)  =  w(z,  e). 


This  changes  the  given  differential  equation  into 

dha  ^  ,  dw  /  dw  \ 

<9-2>  = 

(dw  \ 
z,  w,  — ,  ej  is  also.  The 

boundary  conditions  (1.2)  are  changed  into 

w( 0,  e)  =  lx,  w(((i  —  x)a,  e)  =  /2 . 


It  is  plausible,  but  not  obvious,  that  w(z,  e)  can  be  approximated,  at 
least  for  bounded  z  intervals,  by  some  solution  w0(z)  of  the  differential  equa¬ 
tion 


(9.3) 


dhv  o 
dz 2 


=  Fx{ «,  ®0,  0) 


dz 


Another,  even  less  self-evident,  plausibility  consideration  leads  to  the  bound¬ 
ary  conditions 

(9.4)  ^0(0)  =  4  ,  k/0(oo)  =  y0(a) 

for  this  approximate  solution. 

We  shall  now  confirm  these  conjectures  with  the  help  of  Theorems  2 
and  3.  The  series  (8.8)  for  the  full  solution  y(x,  e)  becomes,  after  the  trans¬ 
formation  (9.1), 

(9.5)  w(z,  e)  =  2  cra(«  +  ez,  4  ,  e )  expj  F^a.  +  ex,  y0(  cl  +  ex),  0)  izl. 

m-0  1  J  0  J 

According  to  Theorem  3  this  series  is  convergent,  uniformly  with  respect  to 
z  and  to  e  in  the  range 

O^z  £(P-  «)ct,  0  ^  e  £  ^ . 

For  e  =  0  it  defines,  therefore,  in  particular,  the  function 


(9.6)  w(z,  0)  =  J  cm(«,  4 , 0)  exp  {mzF^x,  y0(a),  0)}. 

m—0 

Since 

w(0,  0)  =  J  cm(«,  4,0)  =  lim  y(0,  e)  =  4  , 

oi-O  »-H> 


w(z,  0)  has  the  initial  value 

(9.7)  w( 0,  0)  =  4  . 
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As  z  ->  oo,  we  have,  thanks  to  the  condition  (2.3), 

w(co,  0)  =  c0( a,  4 , 0). 

But  by  Theorem  3, 

c0(«,  4.0)  =  w0(a,  0), 

and  by  Theorem  2, 

Mo(a>  0)  =  Vo (a)- 


Hence, 

(9.8)  w(co,  0)  =  y0(ct). 

Thus  w(z,  0)  satisfies  the  two  boundary  conditions  for  w0(z)  in  (9.4). 

The  series  (8.9)  becomes,  after  setting  x  =  a  -f-  ze, 

m  d”  (  r*  1 

2  +  a-  4 .  e)  exP  {  0  Fi(x  +  cr,  y0(a  +  er),  0)  rfr  j  ,  v  =  1,  2. 


This  is  the  termwise  derivative  of  the  series  (9.5);  it  is  uniformly  convergent, 
with  respect  to  z,  in  any  fixed  interval  0  ^  z  ^  zx  .  It  therefore  represents 
d’w  (z,  e)/dz'.  Since  the  convergence  is  uniform  with  respect  to  e  also,  for 
and  as  each  term  is  continuous  in  e  at  e  =  0,  we  have 


(9.9) 


lim 

e-*o 


d'w(z,  e) 
dz’ 


d?w(z,  0) 
dz' 


for  0  5S  z  2S  zt . 

Observe  that  the  formally  analogous  relation 


v  =  1,  2, 


(9.10)  lim  tf(x,  e)  =  y(,)  ( x ,  0) 

e-*-o 

is  not  true  at  *  =  a.  The  reason  is,  of  course,  that  each  term  of  the  series 
(8.9)  is  discontinuous  in  e  and  x  combined  at  e  —  0,  x  =  «.  In  (9.10)  the 
value  of  x  is  held  fixed  independently  of  e,  whereas  in  (9.9)  the  variable  x 
approaches  a  in  accordance  with  the  relation  x  =  a.  +  ez.  Thus  we  are 
dealing  with  two  entirely  different  passages  to  the  limit. 

From  (9.2)  and  (9.9)  we  conclude  that  w{z,  0)  satisfies  the  differential 
equation  (9.3).  We  have  therefore  proved 


theorem  4.  If  y(x,  e)  is  the j  solution  (7.12)  of  the  given  differential 
equation  with  the  boundary  values  (1.2),  then 

lim  y(a  +  ez,  e)  =  w0(z), 

«—►+<> 

uniformly  for  0  ^  z  ^  zx  (zx  arbitrary  but  independent  of  e),  where  w0(z) 
satisfies  the  differential  equation  (9.3)  and  the  boundary  conditions  (9.4). 
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